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I. ABSTRACT 


This report concerns several research tasks related to the origin and evolution of 
planetary atmospheres and the large-scale distribution of volatile elements in the Solar System. 

These tasks and their present status are as follows: 

1. We have conducted an analysis of the volatility and condensation behavior of compounds 
of iron, aluminum, and phosphorus in the atmosphere of Venus in response to published inter- 
pretations of the Soviet Venera probe XRF experiment data, to investigate the chemistry of 
volcanic gases, injection of volatiles by cometary and asteroidal impactors, and reactions in the 
troposphere. 

2. We have completed and are now writing up our research on condensation-accretion 
modeling of the terrestrial planets. 

3. We have laid the groundwork for a detailed study of the effects of water transport in the 
solar nebula on the bulk composition, oxidation state, and volatile content of preplanetary 
solids. 

4. We have completed an extensive laboratory study of cryovolcanic materials in the outer 
solar system. 

5. We have begun to study the impact erosion and shock alteration of the atmosphere of 
Mars resulting from cometary and asteroidal bombardment. 

6. We have developed a new Monte Carlo model of the cometary and asteroidal bombard- 
ment flux on the terrestrial planets, including all relevant chemical and physical processes asso- 
ciated with atmospheric entry and impact, to assess both the hazards posed by this bombard- 
ment to life on Earth and the degree of cross-correlation between the various phenomena (NO x 
production, explosive yield, crater production, iridium signature, etc.) that characterize this 
bombardment. 

Tasks 2 and 3 are principally aimed at understanding the distribution of volatile elements 
and the oxidation state of planet-forming solids in the inner Solar System. Tasks 1, 5, and 6 are 
all targeted on very important geochemical evolutionary processes, largely (but not entirely) 
aimed at integrating the latest emerging understanding of the chemical, physical, and dynamical 
properties and population statistics of the Earth-crossing asteroid and comet populations to as- 
sess the role of impacts of these bodies in volatile injection, shock-processing, and impact ero- 
sion of the atmospheres of the terrestrial planets. 

The purpose of these investigations has been to contribute to the developing understand- 
ing of both the dynamics of long-term planetary atmosphere evolution and the short-term 
stability of planetary surface environments. 
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II. INTRODUCTION 


Within the last few years the importance of impacts in shaping the composition and 
origin of planetary atmospheres has become widely accepted. At the time of accretion, as 
shown by Ahrens and coworkers, impact devolatilization of infalling volatile-rich asteroidal and 
cometary bodies must have been extremely efficient. Very early development of atmospheres 
must have occurred on the terrestrial planets, beginning at a time when only a few percent of 
the masses of Earth and Venus had accreted. 

The presence of an atmosphere during most of the accretionary era can lead, especially 
on the smaller terrestrial planets, to extensive atmospheric loss by explosive (hydrodynamic) 
blowoff. As first pointed out by ourselves, and also argued by Cameron and by Ahrens, the 
blowoff process may very severely deplete the total inventory of volatiles on a planet. G. H. 
Watkins, using a crude blowoff model, argued that more than 99% of the atmospheric mass of 
Mars could be lost during the accretion of the last 1% of the mass of Mars. Substantial loss can 
occur over billions of years even at present impact rates. Thus not only the accretionary era, 
but also the continuing impact history of a small terrestrial planet, may influence its present 
abundances of the rare gases and nitrogen. Recent work by Vickery and Melosh provides a 
strong quantitative basis for modeling of the explosive blowoff process and its application to the 
atmospheric evolution of Mars. 

We first proposed in 1974 that the late infall of cometary carbonaceous asteroidal 
material on Venus could, even in a conservative estimate, replenish water on Venus at a rate at 
least equal to the present-day escape rate. The present inventories of hydrogen (water), sulfur, 
chlorine and fluorine could be entirely provided by later infall, even at the present rate, con- 
tinuing for billions of years. This suggests that the water abundance on Venus may be in a near 
steady-state, not declining rapidly from some enormous initial amount. The discovery of a 
highly enhanced deuterium abundance on Venus by Donahue and coworkers during the 
Pioneer Venus program has been interpreted by several groups in terms of the unidirectional 
escape of a massive original supply of water. Strict logic dictates only that the mass of escaped 
water be at least 100-200 times the present inventory, on the order of 0.1% of a terrestrial 
ocean, but Donahue and several others have concluded that the original Venus water inventory 
was at least equivalent to Earth’s oceans. All such studies have ignored the very important in- 
fluence of water added to Venus by cometary and asteroidal infall. 

Cameron and Benz, Melosh, and others have argued that giant-impact models (Mars- 
sized impactors) for the origin of the Moon would have profound consequences for the early at- 
mosphere: almost instantaneous loss of most of the early atmosphere might have occurred. 
Also, Ahrens has pointed out that the blast wave from much later, smaller surface impacts 
(kilometer-sized bodies) may remove traces of atmosphere from Earth, and Walker has investi- 
gated the effects of the near-cylindrical atmospheric shock propagating from large entering 
bodies before striking the surface (which is ineffectual at removing atmosphere from bodies 
with escape velocities of 10 km/s or higher). Whether Earth impacts in the last 4 billion years 
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can remove significant amounts of atmosphere from Earth or Venus is not yet clear. Ahrens is 
addressing this problem, and Vickery and Melosh have recently presented criteria for the onset 
of explosive blowoff by the blast waves from km-size impactors, showing that they can indeed 
remove atmosphere from Mars, and providing a basis for extension of the theory to Earth. 

We suggested in 1982 that, over the last 4 billion years, the most important consequence 
of cometary and asteroidal impacts on Earth is severe shock heating of large volumes of atmo- 
sphere, which produces extremely high concentrations of nitrogen oxides. We have proposed 
an explicit model to explain not only the scope but also the selectivity of the Cretaceous extinc- 
tion event as a consequence of the toxicity of NO x and the destabilization of carbonates in the 
upper ocean by acid rain. Alvarez has accepted this suggestion as the most probable single me- 
chanism for mass extinction events. 

We have, in response to the interesting problems raised by the original accretion of the 
terrestrial planets and the later impact modification of their atmospheres, concentrated much 
of our attention on combined condensation-accretion modeling of the terrestrial planets (begin- 
ning with Mercury) and on the specific major effects of later impacts on the atmospheric evolu- 
tion on Venus (volatile addition), Earth (shock chemistry and the subsequent atmospheric and 
marine geochemistry), Mars (explosive blowoff), and Titan (shock chemistry). 


III. RESEARCH ACHIEVEMENTS 

We report here, organized according to task, summaries of our recently completed re- 
search tasks: 


1. Water and Cloud-Forming Condensates on Venus 

Achievements: Since our last peer review, we have modeled the time variation of the H 
and D escape fluxes from Venus in the presence of a stochastic cometary impacat flux. We find 
that the D enhancement observed on Venus by Pioneer Venus mass spectrometry is found even 
for model histories in which Venus is formed with no water content. The inference of early 
oceans on Venus suggested by others as the source of a the D:H enhancement is logically in- 
valid. A very wide range of early conditions, some wet and some dry, can be reconciled with 
the observations (4, 7, 1 1). 

We have also found that the strange and apparently unphysical altitude profiles found for 
water vapor by the Venera 11 infrared spectrometer experiment can be explained in part by 
considering the formation of gaseous carbonic acid near the planetary surface and in part by as- 
cribing much of the 0.94 nm opacity to gaseous sulfuric acid. Further, this model shows hope of 
reconciling almost all the conflicting evidence regarding the water abundance and its vertical 
distribution by considering the speciation of water, and by noJ assuming that the free water 
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abundance and the total water abundance are the same. A true total abundance of a few tens 
of ppm for water seems most likely (20). Reprints of this and other recent papers are attached 
as Appendices. 

We have also explored the evidence presented by Krasnopolskii (Icarus .80, 202 (1989)) 
for the presence of compounds of phosphorus and chlorides of iron and aluminum in the Venus 
clouds. So far we find no thermodynamic basis for the reported elements and compounds, and 
propose an explanation of the Venera cloud composition data in terms of an instability in the 
algorithm used to distinguish chlorine, sulfur and phosphorus in the Soviet XRF data (29). 


2. Condensation-Accretion Modeling of the Terrestrial Planets 

Achievements: We have carried out a suite of accretion models for the terrestrial planets 
using a run of mineral chemistry vs. heliocentric distance that is appropriate to conventional 
Solar Nebula models (28). The accretion probability distributions for the planets are modeled 
as Gaussians of variable width. We have reported the application of this model to Mercury (5), 
which presents the most extreme conditions of temperature and proximity to the Sun. These 
models all underestimate the density (metallic iron content) of Mercury, making it essential to 
postulate some fractionation mechanism such as heating by the superluminous early Sun, selec- 
tive accretion due to differences between the physical properties of iron and silicate mineral 
grains, or large-scale erosion of the crust and upper mantle by giant impact events. For the rest 
of the terrestrial plant region, generally reasonable densities, oxidation states, and volatile con- 
tents are found for intermediate values of the Gaussian half-width (HW). Very broad sampling 
functions (HW > 0.5 AU) erase the density and composition differences between the terrestrial 
planets, whereas very narrow sampling functions (HW< 0.1 AU) fail to provide enough 
volatiles. The paper reporting our results is in preparation (28). Several other papers and 
reviews relevant to the molecular, mineralogical, and isotopic composition of the Solar System 
have also been completed (1, 15, 23, 26, 30). 

We have recently compiled compositional data on the Earth and interpreted it in terms 
of condensation-related (volatility-dependent) and differentiation-related processes (especially 
core formation). This model provides strong constraints on the composition of the core, the 
bulk composition of the Earth, and the accretion sampling process that formed the Earth (17). 


3. Water Vapor Transport in the Solar Nebula 

Background: In 1981 JSL had the idea that radial turbulent transport of water in the 
Solar Nebula would efficiently pump water vapor out past the ice condensation threshold. This 
would lead to a buildup in the surface density of solids near Jupiter’s orbit and to progressive 
depletion of water across the asteroid belt, with a consequent dramatic enhancement of the 
stability of reduced materials such as hydrocarbons and enstatite chondrite minerals across the 
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asteroid belt. The actual work was begun seven years ago as a potential dissertation topic by 
William Drew, who followed me from MIT to Arizona, but the project was placed on hold after 
the first year when Drew took a one-year leave of absence. At the end of the second year, 
Drew elected to withdraw from the University of Arizona and enrolled at Cal Tech, where he 
went to work for David Stevenson. After another year. Drew left Stevenson’s employ to work 
for Yuk Yung, abandoning the project. A paper by Stevenson and Jonathan Lunine has since 
appeared, pointing out that turbulent transport of water vapor in the Solar Nebula can lead to 
efficient outward "cryopumping" of water into the Jupiter formation region, and may be 
responsible for producing high enough surface densities of solids to actually cause the forma- 
tion of Jupiter. They did not notice the important role of water depletion in determining the 
oxidation state and volatile content of preplanetary solids in the terrestrial planet region. 
Lunine and JSL have discussed this problem at some length, and one of my students, Melinda 
Hutson, has expressed strong interest in it. She is a reformed meteoriticist with a Master’s de- 
gree earned working with Bob Dodd at Stony Brook. She has extensive experience with micro- 
probe analysis, and has already begun an analytic study of the most relevant class of meteorites 
(low-petreologic-grade enstatite chondrites) in the LPL microprobe laboratory. 

Larry Lebofsky of the U of A has assisted us in interpreting the multicolor photometry of 
the asteroid belt and infrared spectra of C, D and P asteroids as powerful means to test this 
model (8, 10, 12, 18). 

One of these tasks deals with the chemical consequences of radial transport of water 
vapor in the Solar Nebula. Such transport, caused by outward turbulent mixing of water vapor 
across the water-ice condensation threshold near Jupiter’s orbit, establishes a radial gradient in 
the water vapor content of the terrestrial planet formation region. This occurs, as has been 
shown by Lunine and Stevenson, even for relatively poorly mixed nebular models of the sort 
favored by Stevenson. Prinn’s reassessment of turbulent mixing in nebular models has shown 
that radial mixing is much more important than Stevenson’s published models suggest, and 
hence that radial transport of volatiles is even more efficient than previously recognised. 

Given rapid radial mixing, ice will accumulate in greater than solar proportions and water 
vapor will be progressively depleted from the ice condensation threshold inward toward the 
Sun, with the most severe water depletion found closest to the condensation threshold. Thus 
the total C:0 ratio in the nebula increases outward from the Sun, causing a dramatic radial var- 
iation of oxidation state of minerals in the sense opposite to that usually assumed for a com- 
positionally uniform nebula, where the thermodynamic tendency toward progressive oxidation 
of iron during cooling is dominant. Highly reduced minerals, such as the nitrides, elemental 
silicon, and unusual sulfides found in the enstatite chondrites, and carbonaceous matter pro- 
duced by disequilibrating processes in the nebula, both become more stable with increasing dis- 
tance from the Sun when water transport is allowed. We are pursuing three aspects of this 
problem: detailed modeling of the geochemical consequences of water transport in the nebula, 
rough quantitative assessment of the importance of various mechanisms for disequilibration of 
the nebula, and microprobe analysis of selected meteorites (principally enstatite chondrites) to 
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test the predictions of this model regarding stable mineral assemblages. Implications of oxygen 
transport for metal-silicate fractionation during the formation of Mercury are also being inves- 
tigated. 

Melinda Hutson and JSL have begun work on the thermochemistry of the Solar Nebula 
with proper consideration of water vapor transport. Jonathan Lunine is collaborating with the 
modeling of water transport, and JSL has prime responsibility for the chemical modeling of 
mineral stability in the presence of a radial gradient in water vapor concentration. Steffi Engel 
(presently working for Jon Lunine) is further extending this model to the interpretation of com- 
position of cometary ices (19). 

In the context of the same model, we have begun to explore the thermochemistry and 
petrology of the minerals found in low-petrologic-grade enstatite chondrites. JSL takes the 
lead in applying the computer model of the nebular chemistry described above to the predic- 
tion of the mineralogy of highly reduced condensates formed in the water-depleted portion of 
the inner Solar System. Melinda Hutson has been carrying out a suite of analyses, principally of 
enstatite chondrites of low petrologic grade, to test the predictions of the theory, paying special 
attention to carriers of volatile elements (especially S, N, C and P) and to mineral assemblages 
that permit calculation of the oxygen fugacity. Because the stability of organic matter is 
enormously enhanced by reduction of the water abundance, the significance of this model for 
water transport in the nebula for the formation of the carbonaceous chondrites is also being ex- 
plored as part of the same program. 

Application of the same type of model to the zone of silicate condensation would permit 
investigation of the relationships between oxygen transport, metal-silicate fractionation, and the 
origin of Mercury. 


4. I .aboratorv Studies of Low-Temperature Condensates 

Background: Voyager 2 observations of tectonic features on the surfaces of the satellites 
of Saturn and Uranus, as well as observations of the compositions of comet comas, have stimu- 
lated new interest in the chemistry and physics of icy solar system condensates (22). It has be- 
come abundantly clear that cryovolcanic processes have been important in shaping the surfaces 
of icy satellites. Outgassing associated with cryovolcanism has also been fundamental in form- 
ing the atmospheres of Titan, Triton, and (apparently) Pluto. Thus, to understand these satel- 
lite’s atmospheres it is important that we also understand the physical chemistry of condensed 
volatiles in the outer solar system and the nature of cryovolcanic degassing processes. 

Achievements: Laboratory work by J.S. Kargel on the equation of state, rheology, and 
thermal conductivity of ammonia-water liquids has produced data that are of great value in in- 
terpreting the behavior of low-temperature "magmas" in the ice-rich satellites of the outer 
planets (6, 13, 14, 24, 32). The de-gassing history of icy satellites depends on their thermal 
evolution, and thus, the thermal conductivities of the ices and rock which constitute icy satel- 
lites. The conductivities of ice and rock are reasonably well known, but the same is not true of 
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ammonia hydrate. Our measurements indicate that ammonia dihydrate has a thermal con- 
ductivity one-fourth that of water ice for a given temperature. This unexpectedly low con- 

based^on t0 mCTe f e intenor temperatures in icy satellites above those modeled 

based on ammoma-free ice-rock compositions, thus promoting cryovolcanism and degassing. 

, have continued the equation-of-state experiments with existing equipment investieat- 

Z th n e ri :l C ? 0t H dd T "'k 3 ” me,a ' Sa,tS ’ fo ™ ldeh * <° ‘he aqueous ammoniaTiqu^ds. 
The principal product has been a thorough study of the rheological behavior of "cryomagmas" 

m ai^oma-water-x systems (32). This work has resulted, in collaboration with Bob StmmTa 
hS ? ? 11 -° f th f 6rUptlVe StyIeS and viscous behavior of cryomagmatic activity on 

If Prof nf R he J K 0Vian i an T'c °u UV Iaboratory effort was g featl y enhanced by the presence 
o Prof. D.L. Hogenboom, Head of the Department of Physics at Lafayette College who spent 

a semester o his sabbatical (Fall 1990) in our lab. He has' considerable e^erience ta smE 

p ase re ations in ammonia- and water-rich systems at pressures up to about 2 kb His 

work, conducted at no cost to this grant, is being continued in informal Collaboration with u 
since his return to Pennsylvania. ° wun us 

, S^sing history of icy satellites depends on their thermal evolution, and thus the 

thermal conductivities of the ices and rock which constitute icy satellites. The conductivities of 
£ e and rock are reasonably well known, but the same has not been true of ammonia hydrate 

one fourZha> iy . mea f uren,e " ts indicaled that ammonia dihydrate has a thermal conductivity 
<me-fourth that of water rce for a given temperature. This unexpectedly low conductivity would 

l™lT eaSe '"‘ erl0r tem P era, ures in icy satellites above those modelled based on 
ammoma-free ice-rock compositions, thus, encouraging cryovolcanism and degassing. 

Cryovolcanic de-gassing also depends on the solubilities of gases such as N-, and Cft. in 

aTsTom ?’ TT S * * ° IUt TL at l0W tem P eratures ^ and high presses Sgd 
and Strom). Limited data have been available in the literature (e.g., CH 4 in water at 273 K) 

ut not for other compositions with more compelling cryovolcanological interest (e g CH in 
ammonia-water at 176 K). We have obtained gas solubility data for ammSSS 
over a range of interesting pressures and temperatures. 

tin™ C 7° v ° lcamc and associat ed degassing processes also depend on solid-liquid phase rela- 
withD L H!lr SUreS typical of the interiors of icy satellites. As part of our collaboration 

water sv^em W, h° m hi 8 h -P ress ^ Phase equilibria in the ammonia- 

, , y ' e have extended this work to the water-magnesium sulfate system, promoted 

by the great abundance of water-soluble magnesium sulfates in carbonaceous chondrites (33 P ). 


5 - I mpact grosion and the Evolution of t h e Ma r tian Atmnsphprv* 

Background: It has seemed probable since the 1982 studies by JSL and H Watkins that 

impact erosion of the atmosphere of Mars might be the single most impor.anttactoHnt 
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evolution (9). Numerical hydrodynamic models of blast waves from cometary and asteroidal 
impacts have been reported by Vickery and Melosh. They identify a regime of high total ener- 
gy in which all impactors above a certain velocity (which depends on impactor composition) ac- 
celerate most of their own mass and much of the overlying atmosphere to above escape 
velocity. They apply this model to Mars and find that even relatively small, slow impactors can 
expel atmospheric and impactor material at speeds above the martian escape velocity. Im- 
pactors slightly too small or too slow to cause net blowoff administer severe shocks to the atmo- 
sphere and inject cometary and asteroidal volatiles into it. 

Ann Vickery and Jay Melosh have been working on the adaptation of their hydrodynamic 
code to permit simultaneous study of impactor blowoff, atmospheric erosion, and ejection of 
surface rocks, and my group has considerable experience in modeling the chemical con- 
sequences of impact events, including alteration of the chemistry of the retained atmosphere by 
both shock heating and injection of volatiles by the impactor (3, 7, 11, 27, 31; also LPI Tech. 
Rept. 86-07 (1986); GSA Special Paper 120, p. 215 (1982); EPSL22, 239 (1974); Icarus 43, 85 
(1980)). This area is ripe for further development. 


6. Monte Carlo Modeling of the Impact Histories of the Terrestrial Planets 

Achievements: We have constructed and are now testing a Monte Carlo impact model 
for the Earth in which the known physical, chemical, and dynamical properties of Earth- 
crossing asteroids and comets, and the compositions and physical properties of the meteorite 
population that falls on Earth, are all included. 

An important aspect of this work is the inclusion of new mass-frequency data on small 
(roughly 10 to 1000m) near-Earth objects produced by the Spacewatch program at LPL. Their 
discovery statistics clearly show a rapid upturn in object abundances below 100m diameter. 
Ten-meter objects are about lOOx as common as a linear extrapolation of the Belt asteroid pop- 
ulation would predict (D.L. Rabinowitz, Ap. J. submitted, 1992). The new frequency data join 
smoothly into the classical fireball data for meter-sized objects already established by meteor 
observers and photographic networks. Since these small bodies carry 10-100x as much mass 
and energy as has generally been believed, their role in short-term volatile deposition on 
Venus, atmosphere erosion from Mars, and shock -wave processing of the terrestrial atmo- 
sphere is clearly very large. 

The purpose of the model to demonstrate a reasonable time-series of impacts for the pur- 
pose of cross-correlating the impactor mass, kinetic energy, iridium signature, NO x and dust 
production, etc., as well as permitting correlation with cratering and meteorite fall statistics. 
The phenomena of ablation, deceleration, breakup in the atmosphere, explosive blowoff of 
material from high-energy impacts [A.M. Vickery and H.J. Melosh, GSA Special Paper 247, 
289 (1990)], radiative heating of the ground [I.V. Nemchinov and V.V. Svetsov, Adv. Space 
Res. _U, 95 (1991); Effects of Nuclear Weapons (1977)], skipout into heliocentric orbit, [see 
L.G. Jacchia, Sky & Tel. 4g, 4 (1974); Z. Ceplecha, Bull. Astron. Inst. Czech. M 349 (1979); J. 
Borovicka and Z. Ceplecha, Astron. Ap. 257. 323 (1992)], and capture into eccentric Earth or- 
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bit are all included in the program. The model clearly demonstrates how the velocity and com- 
position diversity of asteroidal and cometary impactors produces a poor correlation between 
the various measures of impact severity. In particular, the absence of detectable iridium layers 
associated with many large impacts is clearly shown (31). The model provides reasonable 
predictions of the flux of near-horizontal impactors [P.H. Schultz and R.E. Lianza, Nature 255 , 
234 (1992)]. We have recently begunt to incorporate into the model the injection of volatiles by 
impactors and the shock production of nitrogen oxides and other species by impact blast waves 
[2; 3; R.G. Prinn and M.B. Fegley, EPSL.83, 1 (1987); P.J. Crutzen, Nature 220, 108 (1987)]. 
Neither seems to present any serious problems. We shall then include a Monte Carlo simula- 
tion of impacts on a mathematical model of the human population density distribution of the 
Earth. David Grinspoon has already suggested several ways to improve and extend the existing 
program, and we will continue to collaborate on its operational use, extending his 1988 PhD 
Dissertation research done under this grant (11). The program has been written so that the 
target may be Mercury, Venus, Earth, Mars, or a "custom" planet. It is therefore suitable with 
little modification to modeling exogenic influences on the evolution of all of the terrestrial 
planets. We would like to include the effects of tidal waves from oceanic impactors in the 
hazard assessment, but present understanding of impact tidal wave generation and propagation 
is probably insufficient to permit this [R. Strelitz, Proc. Lunar Planet. Sci. Conf. JJ), 2799 (1979); 
W.B. McKinnon, GSA Special Paper 120, 129 (1982)]. 
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The short lifetime of water on Venus suggests that the water abundance k in a near 
steady-state balance between loss by escape and replenishment by mfaH. In addition, the 
observed deuterium-to-hydrogen ratio on Venus is consistent with a steady state and does 
no^necessarifv imply a past water excess. We present results of a model incorporating a 
^a tk cometao source and nonthermal escape of hydrogen that produces the ob- 
ser^ed waterTbunJance and D/H ratio. The stochastic variability of each «f these quanti- 
S . el ^ . . t K „ lorop We conclude that water on Venus is in a quasi-steady state 

^ or ...» on -I.. ruoo. ^ 

been obscured by a history of random impacts, e iw* Ac.*™ Fr«, i«. 


INTRODUCTION 

Venus is an exceedingly dry. place by ter- 
restrial standards, having a global water in- 
ventory five orders of magnitude lower than 
that on Earth. The aridity of Venus has of- 
ten been explained as due to the loss of a 
much greater primordial endowment of wa- 
ter, perhaps up to a full terrestrial ocean, 
which has been lost due to a runaway 
greenhouse, UV photolysis, and thermal 
and nonthermal escape of hydrogen. This 
model is supported by a deuterium-to-hy- 
drogen ratio 100 times the terrestrial ratio, 
which is interpreted as a residue left by 
mass-selective escape of at least 100 times 
the current water abundance, perhaps a 
much greater amount depending on the effi- 
ciency of deuterium escape (Donahue el al. 
1982). This scenario has been widely ac- 
cepted and has diffused into the popular 
science literature, where it is often pre- 
sented as a parable of a sister Earth gone 

astray. . 

A consideration of the current lifetime 
of water on Venus casts doubt on this 

1 Presented at the Origin and Evolution of Planetary 

and Satellite Atmospheres Conference, Tucson, AZ, 

March 1987. 


monotonic secular decline in abundance. 
Dividing the current hydrogen column 
abundance by the nonthermal escape flux 
yields a characteristic lifetime on the order 
of 10 8 years. Either we have arrived just in 
time to witness the final departure of water 
from Venus or, much more likely, there is a 
source of water, and the abundance is in or 
near steady state. Possible sources include 
volcanic outgassing (Kumar et al. 1983) and 
cometary impact (Lewis 1974, Grinspoon 
and Lewis 1986). 

Models of planetary formation differ on 
how much water and bound water should 
have been incorporated in the original ma- 
terial which accreted to form Venus. It is 
not known whether there was enough helio- 
centric mixing of planetesimals in the for- 
mation region of the terrestrial planets to 
give them similar original volatile invento- 
ries (Wetherill 1985). The observed varia- 
tion of density with heliocentric radius 
among the terrestrial planets would not be 
produced in an extremely “well-mixed ac- 
cretion process. Planet formation may have 
been a more quiescent process with low ec- 
centricities preventing a large degree of ex- 
change (Greenberg 1987). In this case, the 
temperature gradient in the nebula would 
have resulted in a compositional trend 
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among the planets which would provide 
Venus with a much smaller water endow- 
ment than Earth (Lewis 1972). 

We have pursued an extreme model in 
which endogenous water sources are ne- 
glected, in order to explore and illustrate 
the possible importance of exogenous wa- 
ter on Venus. The apparent steepness of 
the size-frequency distribution of comets 
indicates that the bulk of the mass in the 
impacting flux may be contained in a rela- 
tively small number of massive nuclei. This 
suggests the possibility of large stochastic 
variations in water abundance. (The K-T 
impactor on Earth may have been a 10 l9 -g 
comet (Lewis et al. 1982). This would pro- 
vide about 50 9c of the current water abun- 
dance on Venus! The mean interval be- 
tween such impact events appears to be 
similar to the characteristic lifetime for wa- 
ter on Venus.) Accordingly, we have con- 
structed a model incorporating stochastic 
impact injection and nonthermal escape of 
hydrogen, in order to explore the behavior 
of this interesting system. 

WATER ON VENUS 

There is a discrepancy between Venera 
and Pioneer Venus measurements of water 
abundance, but the favored interpretation is 
of a surface mixing ratio of 20 ppm rising to 
a value of 200 ppm at the cloud tops (Young 
et al. 1984). The cause of this gradient is not 
understood. 

The current hydrogen escape flux is dom- 
inated by hot O impact (McElroy et al. 
1982) and charge exchange with hot H' 
(Kumar et al. 1983), yielding a total escape 
flux of 2 x 10 7 H atoms cm' 2 sec -1 . Assum- 
ing that water is the major H-containing 
species, which appears to be the case, then 
the column abundance of 0.1 g of H per 
square centimeter has a lifetime of 10* 
years, quite short compared to the age of 
Venus. 

It is unsatisfactory simply to assume (al- 
though it is conceivable) that the present 
epoch is a unique point in the evolution of 
the atmosphere of Venus. A much more 


likely possibility is that Venus is being sup- 
plied with water at a rate which balances 
the nonthermal escape flux. Two poten- 
tially large water sources suggest them- 
selves: long-term continuous volcanic out- 
gassing, and impact of volatile rich bodies, 
among which icy comet nuclei should carry 
the bulk of the incoming hydrogen. 

Although there is good evidence for vol- 
canic activity at some point in the history of 
Venus (Barsukov et al. 1986), the amount 
of accompanying outgassing is totally un- 
known. It is not known whether this activ- 
ity continues at present, although there 
have been intriguing observations of secu- 
lar changes in the SO; abundance above the 
clouds (Esposito 1984). Estimates of the 
time dependence and total extent of out- 
gassing on Venus must be highly dependent 
on one’s preference of models for the origin 
and geologic evolution of the planet. 

Lewis (1974) found that the flux of com- 
ets and volatile-rich asteroids could be suffi- 
cient to maintain a steady-state abundance 
of hydrogen on Venus. Estimates of the 
cometary flux through the inner Solar Sys- 
tem have been made based on cratering 
records, astronomical observations, and 
dynamical calculations (Weissman 1985, 
Wetherill .and Shoemaker 1982. Ip and 
Fernandez 1987). Although there is also a 
large uncertainty in these estimates, this is 
a known hydrogen source for Venus which 
may be considered independently of geo- 
logic models. 

DEUTERIUM ON VENUS 

The large deuterium-to-hydrogen ratio in 
the atmosphere of Venus was established 
by McElroy et al. (1982), who derived a 
D/H ratio of 1 x 10“ 2 from the mass 2 ion 
detection by the Pioneer Venus ion mass 
spectrometer. Donahue et al. (1982) ana- 
lyzed the PV large probe neutral mass spec- 
trometer measurement of HDO, finding a 
D/H ratio of (1.6 ± 0.2) x 10' 2 . This mea- 
surement was regarded as a critical test be- 
tween theories in which Venus formed in its 
present dry state and those in which large 
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amounts of water were lost over the history 
of the planet. 

For computational simplicity it is as- 
sumed that the escape flux of hydrogen is 
close to the diffusion limit, and therefore 
proportional to the hydrogen abundance, 
and that the escape of deuterium is likewise 
proportional to the deuterium abundance 
multiplied by a factor /, the fractionation 
factor, which represents the relative ineffi- 
ciency of deuterium escape due to its larger 
mass. The differential equations represent- 
ing the escape of these isotopes have the 
solution (Hunten 1982, Kasting and Pollack 
1983) 


/D. / D\ 

‘h'i “ WoL(/ Hi o)|J ’ 


( 1 ) 


where the subscripts 1 and 0 can be taken to 
mean “present” and “primordial,” respec- 
tively, and /h-,o is the mixing ratio of water 
in the atmosphere which we have assumed 
represents the total hydrogen inventory. It 
can easily be seen that, using this formal- 
ism, the interpretation of a 100-fold deute- 
rium enrichment is that the water abun- 
dance has declined by at least a factor of 
100, perhaps much greater if/, the deute- 
rium escape efficiency, was significantly 
larger than zero. 

This interpretation depends upon the as- 
sumption that there was no concurrent 
source of hydrogen during the fractionating 
escape and that the original D/H ratio on 
Venus was equal to the modern terrestrial 
ratio. The latter assumption is questionable 
since the origin of this “standard” is not 
understood and thus its applicability be- 
yond the Earth is not known (Grinspoon 
and Lewis 1987). In this model the fraction- 
ating nonthermal escape was preceded by a 
period of hydrodynamic escape of much 
larger quantities of water, perhaps up to a 
full terrestrial ocean (Kasting and Pollack 
1983). There is little evidence for this hy- 
drodynamic escape phase, beyond the ar- 
gument that some models of planet forma- 
tion would have provided Venus with 


Earthlike quantities of water which then 
would have escaped in this fashion. An 
early period of massive hydrodynamic es- 
cape has also been invoked by Pepin (1987) 
to generate the observed mass fractionation 
of noble gas elements and isotopes. 

If, as suggested in the previous section, 
water on Venus may be in a near steady 
state, rather than in monotonic decline, 
then the D H ratio must be explained in a 
consistent fashion. 

When a time-averaged source term for 
hydrogen, 6. which is equal to the escape 
flux, is introduced we obtain the solution 



g -l 


( 2 ) 


where a is the D/H ratio of the hydrogen 
source, H is the steady-state hydrogen 
abundance, and t is the time elapsed. 

If we assume that D/H in the hydrogen 
source is equal to the primordial D/H, an 
assumption which is perhaps more safe for 
outgassing than for an external source, then 
this reduces to Eq. (10) in Krasnopolsky 
(1985). 

Hunten 11982) and Hunten et al. (1987) 
give the solution 

D/H = alf (3) 


for the steady-state case. This is indeed the 
limiting solution, with the exponential term 
(representing the decaying signature of the 
original D/H ratio) disappearing as t goes to 
infinity. However, it is not obvious that this 
term can be ignored. The time constant for 
decay of this exponential term is 

r = H /<*>/. (4) 


For a hydrogen inventory corresponding 
to 20 ppm water, an escape flux of 2 x 10 7 
cm -2 sec -1 , and / = 0.022 for charge ex- 
change (Krasnopolsky 1985), we find r = 
3.97 x 10 9 years. Hunten et al. (1987) sug- 
gest a weighted fractionation factor, allow- 
ing for the relative magnitude and differing 
fractionation factors of the various escape 
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.h»ni'ms of / = 0.013. Employing this 

mechanisms, o j _ , 71 x in’ years. 

value in (4) gives us r • ems 

Thus the decay time for this term seems 
t^be on the order of the age of the Solar 

^Thesolar EUV flux was more intense ;tn 
,.} TJt (Zahnle and Walker 1982) and thus 
ic seems likelv that the escape flux of H may 
have been at the diffusion limit, or closer to 
the limit than at present, or much ‘ 
nlanet’s history (Krasnopolsky 1985). ^lh 
Lreased flux for a given water abundance 
would have the effect of lowering r by an 
amount proportional to the increase m flux^ 
Uncertainties in estimating the present and 
past vares of*, H, and/make it difficult to 
determine whether or not the limiting sol - 
tion of (2), represented by (3). should 
reached in 4.5 x 10’ years of steady-state 

evolution of the D/H ratio. 

If the source of hydrogen has a i terrest 

D/H ratio and the fractionation factor is as 

modeled by Hunten et al. ( 198/) ( “ _ L 03 
0 _4 f= o.013), then (3) gives D/H - 1-3 
x 10" : The average cometary D/H ratio is 
unknown. Some models of ^ ome ^ 
lead to predictions of e nhanced D/H i 
cometarv ices (Vanysek and Vanysek 1985 
In 1984). Eberhardt et al. (1987) 

Giotto Neutr^ ^ 

x%“”or Come. Halley. Wh«he. or no^ 

this measurement represents the - bulk _ D< 'H 
in Comet Halley and whether Halley should 

^ntair^average cometary D/H are both^in- 

known. Using this range for am (3) 2 

steady-state D/H of 4.6 x 10 to 3. 
for Venus, nicely bracketing the observed 
value. Including the exponential term inU 
with a decay constant on the order of the 
age of the Solar System lowers these values 
somewhat, depending on the chosen i value 
for “original” D/H. Yet it »s sttl possiWe 
with reasonable assumptions for the un 
known parameters a, /, (D ^ )o ’ h g n steady . 
derive D/H values, employing the ste y 
state solution, which are consistent with 
the observations. It is not nece ^ y ^ f 
tulate any past or “primordial excess 


water on Venus to explain the observed 
D/H ratio. 

the stochastic impact injection 
model 

For a power law mass-frequency distri 
bu.ion with a cumulative slope of 

than 2, most of the mass rs coma, ned I m .he 
largest objects in the distribution. The siz 
distribution of comets appears to be of this 
nature (Hughes and Daniels 1980, btioe 
nature (« = 1982 ). Thus it seem> 

stochastic fluctuations in the water abu 

a ratio on Venus. In order io 
dance and D/H ratio un fllictu . 

There is a large uncertainty in our knowl- 
edge of the flux of large bodies through the 
inner Solar System and the resulting fre- 
'"uencS of impact with .he 
ets Thus any conclusion based o 
“chSng such an impact tax must 
v,iwt aeainst large variations in the map 
rnde 5 Z flux', to treat this t— « 

Tideffot tL b mos“L“ta-.ency disuv 

- 

” ofrCeUneC S is an adaptation of the 
of tapactors on the Earth in terms o the 

~t«S£5£ 

%SA — 

and Wasson (1986). We have modified this 
flux to include only the cometary impacts 
and considered the relative probabilities o 
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impact on Venus and Earth for objects in 
cometary orbits. 

Within the limits of observational uncer- 
tainty, it appears that the size distributions 
of asteroids and comets are similar (Hughes 
and Daniels 1980, Watkins, 1983). Thus, we 
represent the flux of comets as a fraction ot 
the total impacting mass, having the same 
size distribution as the overall flux. In light 
of the spacecraft observations of the 1986 
apparition of Comet Halley, which revealed 
a lower albedo and larger object than ex- 
pected, Weissman (1986) now estimates 
that 50% of the recent terrestrial large im- 
pacts are by active cometary nuclei of 
which approximately one-third are long pe- 
riod and two-thirds are short period. For 
long-period comets the impact probability 
on Venus is 77% greater than that on Earth 
fWeissman 1985). For short-period comets 
.he impact probability is 75% greater on 
Earth (Basaltic Volcanism Study Project 
1981). Thus the cometary flux on the two 
planets can be considered to be equal 
within the limits of the other uncertainties 
involved in the flux estimate. 

Based on these considerations we have 
employed a baseline comet flux which is 
50% of the total flux used by Watkins 
(1983). This distribution is shown in Fig. 1. 

A pseudo-random-number generator was 
used to simulate a stochastic flux with this 
size-frequency distribution. The largest 
comet in our baseline distribution has a 
mass of 4.2 x 10 19 g, having a characteristic 
time between impacts of 4.5 x 10 9 years. 
The integrated mass of this flux of comets 
over the age of the Solar System is 5 x 10- 
g. By comparison the mass which Ip and 
Fernandez (1987) estimated to have been 
received by Venus from mass loss from the 
Oort cloud is 6 x 10 : ° g over this same time 
period. This is independent of a much 
larger component of 10- 4 to 10“ g which 
would have been received by Venus in the 
first few hundred million years of Solar Sys- 
tem history due to planetesimals scattered 
during the formation of Uranus and Nep- 
tune (Ip and Fernandez 1987). For a come- 



Fic. 1. Estimated comet flux on \ enus. The formula 
for this distribution is logio(f) = 0.8333 • • - x '°8io 
(jrf} - 6.699, for t in years and .V/ in grams. 


tarv composition of 50% H:0 by mass, an 
influx of 5 x 10-° g over the age of Venus 
would be balanced by an average^ escape 
flux of 2.5 x 10 7 H atoms cm' : sec* , 
within 25% of the present escape flux. 

Kyte and Wasson (1986) have estimated 
the uncertainty in the terrestrial impact flux 
as a factor of 4 at the low mass end of the 
spectrum. For larger masses the uncer- 
tainty is at least an order of magnitude, 
perhaps several. Our baseline flux is con- 
servative in that we have assumed a 
time-independent flux, remaining at cur- 
rently estimated impact rates throughout 
the entire history of Venus. Yet impact 
rates were surely greater in the past. It is 
not clear what fraction of the massive early 
bombardment suffered by the terrestrial 
planets was due to comets, but dynamical 
considerations suggest that the rate of 
comet impact in the early Solar System 
should have been greater than that at 
present by many orders of magnitude (Ip 
and Fernandez 1987). An additional poten- 
tially significant source of hydrogen which 
we have not included is volatile-rich aster- 
oids, particularly that subset of Apollo as- 
teroids that may be extinct cometary nu- 
clei. These objects may have substantial 
volatile reservoirs in their interiors. 

We assume a cometary composition 
which is 5% hydrogen by mass (Whipple 
1984). We further assume that the impact- 
ing nuclei are vaporized and that the hydro- 
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Fig. 2. Hvdrosen escape flux adapted from Krasno- 
polsky (1985). At high/„,o the cold trap results in an 
escape flux which is independent of hydrogen abun- 
dance. We have extended this function to 1ow/h, 0 bv 
assuming a linear relationship. The dashed line is the 
diffusion-limited flux described in the text. 

gen enters the atmosphere and chemically 
equilibrates, mostly ending up in water 
molecules. Both Watkins (1983) and 
Walker (1986) studied impact blowoff of at- 
mospheres and concluded that it is not a 
significant process for Venus. 

~The second essential component of this 
model is the escape flux of hydrogen and 
deuterium. The function <M/h ; o*), which 
describes the dependence of the hydrogen 
escape flux on the water mixing ratio at the 
surface, is quite complex, involving many 
different escape processes which dominate 
in different regimes of homopause hydro- 
gen abundance. / H . and several regimes of 
differing functional dependence of / H on 
f H " 0 s . The latter is controlled at present by 
the sulfuric acid trap and in larger /h 2 o re ' 
gimes by the cold trap. This has been mod- 
eled in detail by Kumar et al. (1983), Kast- 
ing and Pollack (1983), and Krasnopolsky 
(1985). The results of this modeling are dis- 
played graphically in Fig. 13 of Kumar et 
al. (1983) and Fig. lb of Krasnopolsky 
(1985). These results differ somewhat due 
to different assumptions about the behavior 
of the sulfuric acid trap and the nature of 
the cold trap. We have found that these dif- 
ferences are not great enough to affect our 
results significantly. Our model employs 
the escape function shown in Krasno- 
polsky’s Fig. lb. It was necessary to extend 


this function to regimes of water abundance 
below the current abundance, where it has 
not been modeled in detail (since all pre- 
vious models have assumed a declining wa- 
ter abundance). It was sufficient for our 
purposes to assume that in this regime the 
escape flux simply declines in proportion to 
the water abundance, as the results are not 
very sensitive to the detailed shape of the 
escape function in this region. The escape 
function employed is shown in Fig. 2. 

We also performed some runs with es- 
cape occurring at the diffusion limit for all 
water abundances. For this function, also 
shown in Fig. 2, we used the diffusion-lim- 
ited flux as a function of homopause hydro- 
gen abundance 4> — 4 x 10 13 (/h) crn^ ■ 
sec -1 , given by Kasting and Pollack (198 j), 
For the sulfuric acid trap function we used 
the linear function, /h = 2 x 10 (/h; 0 5 ), 

given by Kumar et al. (1983). Krasno- 
polsky's treatment of this function is more 
complex, allowing for the dependence of/ H 
on the H;0/S0 : ratio. The use of the result- 
ing nonlinear sulfuric acid trap function ex- 
plains why the escape function shown in 
Fig. 2 exceeds the “diffusion-limited’' flux 
at the upper end of the sulfuric acid trap 
regime, as the latter was derived using the 
linear assumption. 

The horizontal portion of the curve in 
Fig. 2 represents the cold trap, with an as- 
sumed temperature of 170 C K, which is 
the case for the present-day atmosphere 
(Krasnopolsky 1985). Here the escape flux 
becomes independent of/H 2 o s until the H;0 
mixing ratio rises to =10% at the surface, 
when rapid hydrodynamic escape becomes 
possible (Kasting and Pollack 1983). The 
hydrodynamic escape regime is not treated 
in our model, as the resulting quasi-steady- 
state water abundances remain well below 
this value. 

The escaping hydrogen was integrated in 
million-year increments, with the amount 
of deuterium escaping in each increment 
equaling the product of the number of es- 
caped H atoms, the instantaneous D/H ra- 
tio, and the chosen fractionation factor. In 
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a more rigorous treatment the fractionation 
factor would also vary with water abun- 
dance as the relative importance of the non- 
thermal escape processes changed. 

The number and mass of impacts were 
calculated at 50-million-year increments. 
The cumulative mass of small comets with 
a characteristic time between impacts of 
less than 50 million years was included as a 
constant “low mass trickle" ’ in each incre- 
ment. For the slope index used, this consti- 
tuted 34% of the total flux. A simple way to 
modify this model to include the effect of 
outgassing would be to increase the value 
of this parameter, although there ]s no rea- 
son to expect that outgassed and cometary 
water would have the same D/H. For each 
50-million-year increment the water mixing 
ratio, D/H ratio, cumulative cometary 
mass, and percentage of water of cometary 
origin were recorded. 

MODEL RESULTS 

Figure 3a shows the water abundance as 
a function of time for five model runs em- 
ploying an impact flux equal to twice the 
baseline comet flux. This behavior, spikes 
of high / H ,o corresponding to large impacts, 
interspersed with periods of declining w ater 
abundance, was typical of runs with a wide 
range of integrated impact fluxes and with 
diverse assumptions regarding escape be- 
havior. Figure 3b is a histogram showing 
percentage of time spent in different bins of 
'vater abundance for these same five runs. 
Figures 4a and 4b show the behavior of the 
water abundance for a flux depleted by a 
factor of 2 relative to the baseline flux. As 
can be seen, the behavior is essentially 
identical, with the system simply spending 
a larger percentage of time in the lower wa- 
ter-abundance range. Model runs were per- 
formed with the integrated impact flux var- 
ied over two orders of magnitude, and still 
the behavior was similar with the histo- 
grams becoming more perturbed (Fig. 5). 

The histograms each exhibit a sharp cut- 
off at the low water-abundance end. This 
corresponds to the level at which the steady 


“low mass trickle“ for the chosen impact 
flux just balances the escaping hydrogen 
flux. The low' end cutoff in the histogram 
for many model runs is higher than the ob- 
served 20 ppm water abundance. We do not 
regard this as problematical, as the uncer- 
tainties in this measurement and in many 
model parameters are large. The model is 
mainly meant to illustrate the character of 
the behavior to be expected in this system. 
In most cases however the system spent a 
considerable amount of time within the 
range of uncertainty in the current water 
abundance, w'hich we consider to be 15-100 
ppm. 

Figure 6 show's the behavior of the D/H 
ratio for the same runs displayed in Fig. 5. 
As expected, the stochastic variation of 
D/H is large but the net trend is toward 
larger D/H as the steady state evolves. Fig- 
ure 7 shows the effect of varying the D H of 
the cometary hydrogen and holding all 
other variables constant. Clearly, the evo- 
lution of the D/H ratio is much more sensi- 
tive to choices of unknowm input parame- 
ters than is the / H ,o evolution, which was 
consistent w'ith observations for a very 
wide range of assumptions. Although the 
D/H ratio entered the observed range for 
some reasonable choices of the input pa- 
rameters, there w r ere many combinations of 
total flux, escape function, and cometary 
D/H for which the planetary D/H remained 
w'ell below the observed value. For exam- 
ple. an impact flux depleted by a factor of 8 
relative to the baseline flux did not produce 
sufficiently enhanced D/H unless the come- 
tary D/H was a factor of 10 greater than 
terrestrial. Thus our conclusions about D/ 
H evolution based on this model cannot be 
as robust as those about / H: o evolution. 
These results led us to consider the effect of 
a time-dependent comet flux on the D/H 
ratio. 

Although the model runs described up to 
this point have assumed no higher ancient 
cometary flux, this is counter to both the 
observational evidence contained in the 
cratering record of the terrestrial planets 
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Fig 3 (a) Five model runs for evolution of water abundance on Venus. These runs employed an 
bins of water abundance. 


and the theoretical treatments of outer 
planet formation and cometary dynamics. 
We do not wish here to effect a precisely 
detailed simulation of the time dependence 
of this flux, but merely to explore the con- 
sequences of a past excess in impact rate 
which mimics to some extent those which 
have been postulated in the literature. 
Shoemaker and Wolfe (1982) suggested that 


the late heavy bombardment of the inner 
Solar System was due to comets, but Strom 
(1987) has offered arguments to the con- 
trary based on detailed comparisons of the 
cratering records of the planets. We simply 
assume that the proportion of cometary ob- 
jects in the total impact flux was the same in 
the past as it is at present. 

According to the model of Hartmann 
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Fig. 4. (a) Results of five runs employing an impact flux depleted by a factor of 4 relative to the runs 
in Fig. 3. (b) Water abundance histogram for these runs. 


(1987). the terrestrial impact flux 4 billion 
years ago was approximately 10 3 times the 
current flux, although the earlier, postac- 
cretional flux must surely have been much 
higher. We modified our model to include 
an enhancement of the impact rate by 10 3 at 
; = -4 x 10 9 , decreasing with a 150-million* 
year half-life. This gives a total integrated 
flux of 2.8 x 10** g. Compared to the 10 24 to 
10 25 g which Venus may have received from 


scattered planetesimals during the forma- 
tion of Uranus and Neptune (Ip and 
Fernandez, 1987), this is still a conservative 
flux estimate. The effect of this time-depen- 
dent flux on the evolution of the D/H ratio 
is shown in Fig. 8. Several more modest 
time-dependent fluxes were also run. Also 
shown in Fig. 8 is a run where the baseline 
flux is depleted by a factor of 4 and the 
enhancement at t - -4 x !0 9 years is only a 
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Fig. 5. (a) Five model runs abundance histogram for these runs. 








:tor of 500, again decaying with a half-life 
150 million years. The integrated mass o 

s flux model is 3.3 x 10 :i g. 

When the effects of such time-dependent 
1X es are included, then in nearly every 
se the current D/H ratio on Venus is ex- 
uded bv the model results. In these runs, 
e large final D/H values are due to a corn- 

nation of early fractionating water loss 
id later steady-state evolution. It the 


steady state were allowed to opera e for 
long enough (several r) then the final D/H 
reached would be independent of this earlv 
historv, but since r is on the order of the 
age of the Solar System the resulting D i 
ratios are quite enhanced after 4.5 billion 
years The water abundance evolution for 
these time-dependent runs is shown in Fig. 
9 This scenario differs from the previoush 
published models of early fractionating wa- 
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ter loss in one important respect: here none 
of the water excess is endogenous. The 
large early water abundances are simply the 
steady-state abundances corresponding to 
the massive early cometary influxes. The 
peak water abundances, which are off scale 


in Fig. 9, are approximately 27,000 and 
8000 ppm for the two time-dependent fluxes 
described above. As the cometary flux de- 
clines exponentially to the present-day 
value, the water abundance declines until 
the escape flux balances this lower influx. 



Fir 7 The effect of cometary D/H on Venus D/H evolution. The lower curve is taken from Fig. 6, 
Fig. 7. The effect ot y ^ ^ difference jn scale The upper curve is the same run, 

using°a cornetary D/H value of 4.8 x 10-', the observed upper limit for Comet Halley (Eberhardt et a . 

mi). 
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Fic. 8. D'H evolution on Venus for son. 
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The final results are water abundances 
which are compatible with the observed 
abundance on Venus and D/H ratios which 
are quite enhanced. 

DISCUSSION 

A consideration of the current lifetime of 
water on Venus and our model results lead 


us to conclude that the water abundance on 
Venus is likely to be in a quasi-steady state, 
mediated by episodic large comet impacts. 
Any remnant of an early Earthlike ocean on 
Venus is obscured by a history of random 
impacts. Known processes provide an ex- 
planation of the observed water abundance 
and D/H .on Venus. There is no need to 



Fig. 9. Water abundance evolution for the same runs shown in Fig. 8. The peak abundances, which 
are off scale, are discussed in the text. 
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postulate an unknown source of early water 
^nd no evidence that there was one. Out- 
gassing may have contributed as well, but 
no endogenous source of water is necessary 
to explain the observations. The water now 
observed on Venus is quite possibly more 
than 99% of cometary origin. 

The observed deuterium to hydrogen ra- 
tio may not be evidence of a “primordial 1 
vater excess on Venus. If an early water 
xcess is responsible in part for the current 
D H, it is likely to have been the steady- 
state water abundance for an enhanced 
cometarv flux. This ratio is coupled in a 
complex relationship to the average come- 
tary D/H, the time-integrated cometary flux 
and any time dependence in the flux, and 
the actual mass and D/H of any large com- 
ets which may have impacted on Venus in 
ne last several hundred million years. 

In a similar fashion, the current water 
abundance on Venus may depend critically 
on the recent comet impact history of the 
planet as well as on the average flux and 
outgassing history. This stochastic variabil- 
ity. as illustrated in Figs. 3 through 5, must 
add an unavoidable level of uncertainty to 
any deterministic interpretations of these 
observations. 

It is possible that these stochastic varia- 
tions in water abundance will be damped by 
surface-atmosphere buffering reactions of 
the kind described by Nozette and Lewis 
(1982). Characterization of the time scale 
for this buffering is difficult as it depends on 
such poorly known parameters as diffusion 
rates through the surface material, particle 
sizes, and the rates of many heterogeneous 
reactions which have not been well studied. 
If there is not significant damping, these os- 
cillations in water abundance may have im- 
portant implications for the climate history 
of the planet; water is an important IR ab- 
sorber which helps to maintain the strong 
greenhouse. 

When a consideration of the likely time 
dependence of the comet flux is included in 
the model, then the problem of explaining 
the current D/H ratio becomes one of ex- 


plaining why it is so low compared to model 
results. This could put constraints on the 
magnitude of any past increase in the come- 
tary flux. However, the attempt to define 
such constraints will be hampered by the 
difficulty in determining how the fraction- 
ation factor varies with escape flux. The 
fractionation factor is likely to be higher 
than at present for enhanced escape fluxes 
(Krasnopolsky 1985, Kasting et al. 1984). 
This would lower the final D/H for the time- 
dependent case, possibly reconciling the 
enhanced values shown in Fig. 8 with the 
observed values. Chyba (1987) has used 
scaling arguments to conclude from the lu- 
nar impact record that if a small fraction of 
the late heavy bombardment impact flux 
was composed of cometary objects then the 
Earth's oceans could have been fully sup- 
plied by cometary water. Under such a sce- 
nario Venus should have received similar 
amounts of water, yet the low D/H may be 
problematical. A continuous outgassing of 
juvenile water, as proposed by Kumar et al. 
(1983) and Kasting and Pollack (1983), 
could help to lower the D/H. As can be 
seen by examining Eqs. (2) and (4), this 
would have the effect of lowering the time 
constant for the steady-state D/H value to 
be reached. In this case the steady-state 
value would be that given in Eq. (3) with a 
being a weighted average of the outgassed 
and cometary D/H values. Another possi- 
ble explanation is that Venus has recently 
(within 10 s years) suffered a large comet im- 
pact which lowered the D/H ratio to the 
observed value. If the terminal cretaceous 
event on Earth was caused by a comet 
shower, rather than by a single large im- 
pact, then the water abundance and D/H 
ratio on Venus should bear the signature of 
this same catastrophe. 
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ttc presemresuits ot an extensive laboratory and telescopic investigation of HvO distribu* 
hon among the low-albedo, outer belt asteroids (2.5-S.2 AU). The water distribution was 
determined by surveying asteroids in that region for the 3- M m reflectance absorption of 
molecular H,0 and structural OH ions; the 3- jam band is the spectral signature of meteorite 
and asteroid hydrated silicates. The 19 asteroids observed in the survey, augmented bv 
ear '^reflectance data, yielded 3-/zm band depth measurements that reflect the present 
outer belt H,0 distribution. Of the 32 C-class asteroids in this sample, 66% have hvdrated 
silicate surfaces, evidence of a mild aqueous alteration episode earlv in solar system history 
the C class is not a primitive asteroid group in the mineralogical sense, but its anhvdrous 
members appear little altered, as do the P and D classes beyond 3.5 AU. In addition'to this 
pronounced hydration state difference between Cs and the more distant Ps and Ds, hvdrated 
silicate abundance declines gradually among the Cs from 2.5-3.S AU. Coupled with the 
apparently anhydrous P and D surfaces, this points to an original outer belt asteroid composi- 
tion of anhydrous silicates, water ice, and complex organic material. Earlv solar wind induc- 
tion heating of proto-asteroids declined in intensity with heliocentric distance, and produced 
the observed radial decrease in hydrated silicate abundance. The mild thermal process- 
ing of outer belt asteroids is a continuation of the intense heating and differentiation that oc- 
curred sunward of 2.5 AU, and both events are consistent with the induction heating mech- 
anism. Water ice may be preserved in the interiors of P and D objects, and in some larger 

US* © 1990 Academic Press, Inc. 


I. INTRODUCTION 

Asteroids whose orbits have semimajor 
axes, a, > 2.5 AU are predominantly mem- 

| Visiting Astronomer at the Infrared Telescope Fa- 
cility which is operated by the University of Hawaii 
under contract to the National Aeronautics and Space 
Administration. 


bers of the C, P, and D taxonomic classes 
We refer to these as outer belt asteroid.' 
since these dark objects dominate the popu 
Iation from near the center of the main bei 
at 2.5 AU all the way to Jupiter's orbit a 
5.2 AU. Spectral observations reveal tha* 
some have surfaces dominated by low-tem- 
perature, water-bearing minerals. We alsi 
infer a substantial organic component from 
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Fig. 1. Spatial distributionofthe 38 low-albedo aster- 
oids included in this 3-p.m reflectance survey. The sam- 
ple spans the mid- to outer portion of the main belt, a 
region dominated by dark asteroids. 


th - uniformly low albedos. Both traits are 
typical of relatively unaltered meteorites 
like the carbonaceous chondrites. If these 
dark asteroids do retain their original miner- 
alogy, then they should yield direct informa- 
tion on both the solids produced in this por- 
tion of the solar nebula, and the early 
chemical and physical processes operating 
in the asteroid belt. 

This paper discusses the results of a multi- 
year reflectance survey of low-albedo aster- 
oids, using the 3-^m hydrated silicate ab- 
sorption to define the H,0 distribution in the 
outer belt, (see Fig. I). Hydrated silicates, 
which are thermodynamically stable at rela- 
tively low temperatures f<500 K), are trac- 
ers of early thermal activity on primitive 
asteroids. The detection of hydrated sili- 
cates (via 3-/*m spectrometry) on individual 
as: .roids, and their overall outer belt distri- 
bution should reveal important characteris- 
tics of heat sources that acted on those bod- 
■es just after formation. 

Previous Work on Low-Albedo 
Asteroid Composition 

This work is part of a larger effort to dis- 
cover asteroid compositional trends and es- 
a t :ish firm links between meteorite classes 
and likely parent bodies. Such composi- 
uonal links are important: if remote sensing 
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can reveal an asteroid's present surface min- 
eralogy, via analogy to meteoritic “ground 
truth,” one can establish a tentative “phys- 
ico-chemical history" for that body (Gaffey 
1987). Such an approach will be, of neces- 
sity, the only method of establishing individ- 
ual compositions for large numbers of indi- 
vidual asteroids. 

The asteroid belt displays a systematic 
variation of spectral class with heliocentric 
distance (Gradie and Tedesco 1982). The S 
class dominates the inner belt, the Cs the 
middle and outer part of the main belt, and 
the Ds the Trojan region. When linked to 
surface mineralogy by spectrophotometry, 
the class variation is clear evidence for over- 
lapping “rings” of different asteroid compo- 
sition at increasing distances from the Sun. 
If, on the basis of meteorite analogies, one 
equates the S, C, and D classes, respec- 
tively, with volatile-poor silicates, carbo- 
naceous chondrite-Iike material, and vola- 
tile-nch “ultracarbonaceous” material 
(Gradie and Veverka 1980), the radial struc- 
ture represents the mineralogical record of 
belt formation temperatures, or the effects 
of some radially varying thermal event early 
m solar system history, or both. 

Because of their spectral similarity to the 
chemically primitive carbonaceous chon- 
drites, the C asteroids have been regarded as 
substantially unaltered, volatile-rich bodies, 
darkened by opague organic material that 
dominates their reflectance properties 
(Chapman et al. 1975). In the outer belt, the 
P and D classes were thought to be even 
ncher in organics and hydrated clay miner- 
als, unaltered ultracarbonaceous objects 
formed at low temperatures in the solar neb- 
ula (Gradie and Veverka 1980). Mid-IR 3- 
£im observations of low-albedo asteroids 
confirmed the presence of surface hydrated 
silicates on several Cs (Lebofsky 1978). On 
the basis of this work, Feierberg et al. 
(1985a) suggested that the varying abun- 
dance of hydrated silicates among the Cs 
would be typical of a population “of colli- 
sional fragments: the Cs were pieces of orgi- 
nally “wet” (hydrated silicate) parent bod- 
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ies whose cores had been heated and 
dehydrated, then shattered. We set out to 
test this premise in 1986, but our telescopic 
observations and 3-pun reflectance measure- 
ments of carbonaceous chondrites led us to 
a quite different conclusion, that the varia- 
tion in hydration state among the dark aster- 
oids reflected different thermal histories 
controlled by composition and solar dis- 
tance (Jones et al. 1987, 1988a, b, Jones 
1988). 

Rationale for 3-fjLrn Asteroid Observations 

Hydrated silicates are the dominant min- 
eralogical component of carbonaceous 
chondrite matrix material, or CCMM (Wil- 
kening 1978). Because CCMM is black, and 
comprises the bulk of Cl and CM chondrite 
matrices, it is thought to be prevalent on the 
surfaces of primitive outer belt asteroids. 
Ordinary clays have strong H : 0 absorption 
bands at 1.4 and 1.9 pun in the near-IR. Un- 
fortunately, the reflectance properties of 
CCMM are dominated at visual/near-IR 
wavelengths by an opaque, complex organic 
material known as kerogen. The opaque ker- 
ogen lowers the visual albedo and reduces 
spectral contrast. Shortward of 2.5 pan, or- 
ganic absorptions effectively overwhelm the 
relatively weak hydrated silicate overtone 
and combination bands. Despite the pur- 
ported abundance of hydrated silicates in 
the outer belt, visual and near-IR observa- 
tions have been incapable of identifying 
them. Only recently have high quality, high 
resolution reflectance spectra (0.5-0. 9 pun) 
revealed weak Fe 2 "-Fe 3 ~ charge transfer 
features (indicating phyllosilicates) on sev- 
eral dark asteroids (Vilas and Gaffey 1989). 
The H : 0 and structural hydroxyl absorp- 
tions in hydrated silicates at 1.4 and 1.9 pirn 
remain invisible in dark asteroids, due to the 
presence of organics. 

Fortunately, beyond 2.5 pt.m, the spectral 
signature of water is sufficiently strong that 
it dominates the spectrum and allows us to 
distinguish mineralogical differences typical 
of the major carbonaceous chondrite 
classes. The 3-pim hydrated silicate band is 


a combination of (1 ) a sharp 2.7-ju.m fea: 
due to structural OH ions, and (2) a m 
broader 2.9-pun absorption due to interlav 
H : 0 molecules. The wing of the 2.9-pt: 
band overlaps the centers of the 3. 1-ptm w, 
ter ice feature and the 3.4-ptm C-H stretc 
vibration absorption. Since Lebofsky's in 
tial discovery of a 3-pim band on 1 Cere 
(Lebofsky 1978). this feature has been tf 
target of increasingly sophisticated obs * 
ing programs, both to determine the exur 
of hydrated silicate presence in the belt, an 
to characterize the surface mineralogy 
individual asteroids (Larson and Veede 
1979, Lebofsky 1980. Lebofsky et al. 198; 
Feierberg et al. 1981. 1985a, Larson et a 
1983). 

Feierberg et al. (1985a) concluded that 
of the 10 largest C asteroids have a hvdr. : 
silicate surface component. With hydra. e 
silicates apparently a common surface mate 
rial on the Cs, we decided to establish 
clearer picture of the H : 0 distribution 
through a broader survey of low-albedo a^ 
teroids. In addition, use of a higher spectr; 
resolution would pay dividends in greats 
characterization of the 3-pcm band and po<s 
ble detection of ice and organic absorptu ^ 

We obtained higher resolution and spee 
using the new Cooled Grating Array spei. 
trometer (CGAS) at NASA's Infrared Tek 
scope Facility (IRTF). (See Tokunaga an 
Smith [1985] for details). Instead of the sir 
gle-detector CVF system, the CGAS uses 
linear 32-detector InSb array to sample th 
infrared spectrum simultaneously over 
span of nearly 0.6 pun. CGAS can prod c 
a respectable spectrum faster and at thre 
times the resolution of a CVF system O' 
AX *■ 170, vs 25-50). The improvement i 
resolution and SNR with CGAS is apparer 
in spectra of 1 Ceres taken with both instru 
ments (Figs. 2 and 3). CGAS is ideal for r 
pun surveys of moderately bright asteroid' 
scanning 0.6 pun at a glance, the array c a 
measure reflectance on both sides of the c 
luric “water gap." where radiation froi 
2. 5-2. 8 pirn is almost totally absorbed b 
atmospheric H : 0. With such a simultaneot 
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Fig. 2. Spectra of I Ceres using the circular variable 
filter (CVF) photometer and the Cooled Grating Array 
spectrometer (CG AS) at the NASA IRTF. Mauna Kea. 
The marked drop in reflectance from 2.2 to 2.8 ptm is 
due to hydroxyl and water of hydration absorptions. 
Signal loss from 2.5 to 2.8 y . m is caused by atmospheric 
H : 0. Note increase in resolution and spectral range 
w : 'h CGAS. (Spectra are normalized to 1.0 at 2.2 ^m 
(C Ft and 2.5 fim {CG AS), then offset for clarity. 
CVF spectra have had a calculated thermal component 
removed. Error bars = I cr). 


measurement, CGAS yields a far more accu- 
rate 3-M-nrt band depth measurement than a 
sequential CVF scan. 

II. CVF and CGAS OBSERVATIONS— 1986-88 

The observing program goal was to collect 
3-Mm reflectance data for as large a group of 
low-albedo asteroids as possible over a span 
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Fig. 3. CVF spectra of 1 Ceres (G class). Reflectance 
difference between 2.2 ptm (K filter) and 2.9 is the 
basis for band depth calculations. The 3.1-/x m absorp- 
tion band is also seen in CGAS spectra, and is probably 
a n H ; 0 frost feature (Lebofsky 1980). 
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TABLE I 

CVF Observation Log: NASA Infrared 
Telescope Facility, Mauna Kea. Hawaii 


Instrument: RC-2 photometer. CVF 

Observers: M. A. Feierberg 

L. A. Lebofsky 

Date Object Spectral Standard 

class" star 

02/18/86 511 Davida 

C 

Hyd 64 

1 Ceres 

G 

35 Leo 

02/19/86 511 Davida 

C 

Hyd 64 

1 Ceres 

G 

35 Leo 

Instrument: RC-1 

Observers: L. A. 

T. D. 
E. F. 

photometer. CVF 
Lebofsky 
Jones 
Tedesco 

12/08/86 19 Fortuna 

G 

Hyd 64 

10 Hvgtea 

C 

Hyd 64 

704 Interamnia 

F 

35 Leo 

12/09/86 lUKassandra 

T 

Hyd 64 

10 Hygiea 

C 

Hyd 64 

704 Interamnia 

F 

35 Leo 


a The C class includes subclasses B. F. G. and T. 


of approximately 2 years. We conducted all 
observations from the NASA IRTF, on the 
dates indicated in Tables I and II. 

Table III lists the low-albedo asteroids 
that have been observed to date at 3 m m. 
(Piscitelli et ai [1988] report CGAS 3 -fxm 
observations, but their 3. 1-3.6 fim coverage 
does not include the entire water of hydra- 
tion band.) With CGAS we were able to 
cover the brightest of the Cs (and its related 
subclasses) down to a visual magnitude 
~ 12. Unfortunately, members of the P and 
D classes are usually even fainter, and sky 
background saturation of the CGAS detec- 
tors prevented observations of these inter- 
esting groups. To increase the sample size 
and better characterize the hydrated silicate 
distribution, we added sensitive CVF obser- 
vations by Lebofsky (1980) and Feierberg et 
al (1985a, b). 

The data reduction scheme is common to 
reflectance spectrophotometry (see Lebof- 
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TABLE II 


CGAS Observation Log: NASA Infrared Tele- 
scope Facility, Mauna Kea, Hawaii 


Observers: T. D 

Date Object 

. Jones, L. A. Lebofsky 

Spectral Standard 

class star 

04/21/87 

532 Herculina 

S 

70 Vir 


1 Ceres 

G 

16 Cyg B 


511 Davida 

C 

70 Vir 


2 Pallas 

B 

\ Serp 

04/24/87 

5 Astraea 

S 

35 Leo 


511 Davida 

c 

70 Vir 


18 Melpomene 

s 

SAO 120107 


2 Pallas 

B 

SAO 65083 


1 Ceres 

G 

16 Cvg B 

04/25/87 

704 Interamnia 

F 

35 Leo 


16 Psyche 

M 

35 Leo 


511 Davida 

C 

70 Vir 


2 Pallas 

B 

SAO 65083 


1 Ceres 

G 

16 Cyg B 

Observers: T. D. Jones, M. S. Martey 

10/01/87 

92 Undina 

M 

SC 8167 


1 1 1 Ate 

C 

SC 8167 

10/02/87 

324 Bamberga 

CP 

Hyd 142 

10/03/87 

375 Ursula 

C 

51 Peg 


423 Diotima 

C 

SAO 147237 


51 Ncmausa 

CU 

KCeti 


52 Europa 

CF 

Hyd 142 


13 Egeria 

G 

Hyd 142 

Observers: T. D. 

Jones, L. A. Lebofsky 

02/01/88 

55 Pandora 

M 

Hyd 106 


13 Egeria 

G 

Hyd 106 


10 Hygiea 

C 

35 Leo 

02/02/88 

313 Chaldaea 

C 

Hyd 64 


130 Elektra 

G 

SAO 136389 


505 Cava 

FC 

35 Leo 

02/03/88 

55 Pandora 

M 

Hyd 64 


173 Ino 

C 

SAO 136389 


409 Aspasia 

cx 

70 Vir 

02/04/88 

88 Thisbe 

CF 

Hyd 64 


324 Bamberga 

CP 

Hyd 64 


sky 1978, 1980, Chapman and Gaffey 1979, 
Feierberg et al. 1985a). We chose standard 
stars whose spectra were close to solar 
(Hardorp 1980, 1982): the ratio of asteroid 
reflectance to ( solar analog spectrum re- 
moves solar absorptions and displays only 
mineralogical features. The resulting spec- 
tra plot relative (i.e., scaled) reflectance vs 
wavelength; the reflectance is “relative” to 
that of the standard star, and normalized 
to 1.0 at 2.2 or 2.5 jam. Spectral quality 
depended mainly on intrinsic asteroid 
brightness and viewing geometry; objects 


near M v = 12 proved to be difficult targe’s. 
For example, 13 Egeria produced a g>. v 
spectrum in October 87 (Fig. 4). but by Feo- 
ruary 88 it had moved away from opposi- 
tion, and 3-/xm band depth measurements 
became impractical. Note the similar loss of 
quality for 324 Bamberga in Fig. 4. 

The CVF error bars represent the stan- 
dard deviation (cr) of the brightness at a 
given wavelength, propagated through the 
averaging of several complete CVF spec a 
on the same night. For the CGAS spectra, 
the 1 -<t error bars represent the standard 
deviation of the mean detector voltage (an- 
alagous to brightness) at a given wave- 
length, similar to the CVF case. 

A significant data reduction problem en- 
countered at 3-/xm is the removal of asteroid 
thermal flux. Beyond 3 /xm, emitted ther-al 
energy from the asteroid begins to rival e- 
flected sunlight at the telescope. A “stan- 
dard thermal model" developed by Lebof- 
sky et al. (1986) yields an estimate of the 
emitted/reflected flux ratio at a given wave- 
length. The emitted flux for main-belt ob- 
jects is negligible at 3 /xm, contributes 
5-10% of the signal at 3.4 /xm, and domi- 
nates beyond 3.5 /xm. Removal of the calcu- 
lated therrpal flux does a reasonable jot of 
restoring the spectrum to a normalized re- 
flectance value near 1 .0 at 3 .5 /xm (similar to 
carbonaceous chondrite spectral behavior), 
and so we subtracted a calculated thermal 
component from all CVF spectra (Fig. 5). 
However, in the case of 511 Davida and 19 
Fortuna, the model appears to overcompen- 
sate for thermal emission and unrealistica'ly 
lowers the 3.4-/xm reflectance (see Feier- 
berg et al. 1985a). Most of the CGAS spectra 
were free from significant thermal effects at 
3 /xm, so we chose to leave the CGAS spec- 
tra uncorrected. (The problem is not severe 
for most objects in the outer belt, where a 
lower solar irradiance reduces overall ther- 
mal emission and shifts the peak radiation to 
longer wavelengths.) In the few cases where 
thermal contributions were important, we 
did correct the computed 3-/xm band depth. 
To reiterate, the CVF spectra published 
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TABLE III 

Asteroids Observed at 3 /im 


1 Ceres 

2 Pallas 
10 Hygiea 
13 Egeria 
■9 Fortuna 

► Themis 
jl Euphrosyne 
36 Atalante 

51 Nemausa 

52 Europa 
65 Cybele 
70 Panopaea 
72 Feronia 
74 Galatea 

87 Sylvia 

88 Thisbe 

1 Kassandra 
j-0 Elektra 
139 Juewa 
148 Gallia 
173 Ino 
233 Asterope 
247 Eukrate 
308 Polyxo 
313 Chaldaea 
324 Bamberga 
344 Desiderata 
" "5 Ursula 
j 86 Siegena 

409 Aspasia 

410 Chloria 
423 Diotima 
505 Cava 
511 Davida 
554 Peraga 
704 Irueramnia 
721 Tabora 
773 Irmintraud 
"^6 Berbericia 


2.768 
2.773 
3,138 
2.576 
2.442 
3.129 
3.148 
2.749 
2.366 
3.095 
3.428 
2.614 
2.266 

2.781 
3.483 

2.769 
2.678 
3.109 

2.782 
2.771 
2.745 
2.660 
2.742 
2.749 
2.376 
2.685 
2.591 
3.134 
2.894 
2.577 
2.724 
3.069 
2.685 
3.181 
2.374 
3.060 
3.554 
2.857 
2.934 


Tholen 

class 


Tedesco a 

class 


Controls 


4 Vesta 

2.362 

V 

r 

5 Astraea 

2.577 

S 

S 

16 Psyche 

2.922 

M 

M 

18 Melpomene 

2.296 

S 

S 

55 Pandora 

2.761 

M 

E 

92 Undina 

3.205 

M 

M 

532 Herculina 

2.774 

S 

S 


Band depth 

W 

21 - 1.3 
18 1,2 

12 3.7 

40 15 

25 6 

4-2 
0 1 
14 1 

37 15 

0 6 
0 - 1 
35 12 

8 5 

2 5 

1 10 
0 r 10 
0 5 

22 13 

0 2 
0 9 

0 ± 7 
0 2 
2 2 
14 2 

22 15 
6 ± I 

18 4 

18 14 

23 2 

11 6 

24 ± 3 

0 60 

5 11 

12 4 

0 

6 ± 4 

0 8 

0 6 

23 2 


37 22 

45 12 


Note. Source notes: (I) this work; (2) Feierberg et a!. 1985b; (3) Feierberg ti al. 1985a; (4) Lebofsky 1980 
0 See Tedesco et al. 1989 
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F'g. 4 CGAS spectra of 13 Egeria and 324 Bam- 
rga. showing the deterioration in signal quality due 
to changes ,n observing geometry. Errors in band depth 
measurement become large, and in some cases only a 
qualitative determination of hydrated silicate presence 
possible. Observation are dates listed in Table II 


molecular H : 0 absorptions. We see it ,, 
these spectra as a distinct reflectance dror 
from 2.5 to 2.9 /un (e.g.. Fig. 2). Telluric 
water vapor obscures the signal from 2.5 to 
-.8 fim. The CVF filters precluded coverage 
between 2.2 and 2.5 (Fig. 5), but carbo- 
naceous meteorite spectra are generally flat 
in this region, so we assumed that the J 
0.25 Atm), H (1.65). and K (2.2 M m) O “ 
measurements (only K is shown) define a 
continuum whose peak reflectance corres- 
ponds to the 2.5-p.m CGAS value. While 
most of the CVF spectra lack the S/N 
and/or resolution (A/AA = 50) to reveal de- 
tailed absorptions, the broad 3- M m band 
lends itself to easy detection. We define the 
3-Atm band depth as 


here have a thermal component removed, 
the CGAS spectra do not, and the 3-um 
band depths listed in Table III are corrected 
tor thermal contributions where appro- 
priate. FF 

. H,0 band is a combination of 

the _. 7-A<.m structural hydroxyl and 2.9-Aim 



asteroid* Va ? abie fi!ter sp€ctra of low albedo 

asteroids. The thermal contribution to the flux at longer 

wavelengths has been modeled and removed. Gap from 

Tahir if ’“’l ‘ S dUC t0 te " Uric Water abs °rption. See 
Table I for observation dates. 


CGAS: d = /?:; - * 

25 Aim K2.9fj.mi 

where R is the reflectance at the specified 
wavelength. 

In calculating d, we excluded data be- 
tween 2.6 and 2.9 ,um, since they are either 
noisy or unusable there due to telluric H.O. 
The same water vapor also obscures the 
deepest portion of the intense 2.7-p.m h 1 - 
droxyl banded thus underestimates the re- 
flectance drop from 2.7 to 3.0 M m (i.e. the 
true band depth is > d). Our calculated 3- 
A*m band depths and associated errors (the 
latter based on the uncertainty in 2.9-p.m 
reflectance data) appear in Table III. 

SelecTed Asteroid Spectra 

The bright asteroids I Ceres and 2 Pallas 
yielded the most detailed 3- M m Ipectra 
(Figs. 6 and 7). 

I Ceres (Fig. 6) displays what may be 
a distinct, optically thin, 3.1-Atm water ice 
absorption superimposed on the long-wave- 
length wing of the 3-p.m band (Lebofsky et 
al. 1981). This feature persists from night to 
night, but we saw no real structure corre- 
sponding to the 3.4-pm organic (C-H 
stretch) absorption. The declining contin- 
uum at 2.5 /4m, not seen in carbonaceous 
chondrite spectra, is peculiar to Ceres. The 
3-/xm band depth d is 0.21 (±0.01). 
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Fig. 6. Cooled Graling Array spectrometer spectra 
of I Ceres taken over several nights in Apni 1987. Note 
the prominent 2.5-3.0 /xm absorption, probable ice 
band at 3. 1 /xm. and the methane artifact at 3 33 
We observed no 3.4-/xm C-H stretch absorptions. On 
-I April we were unable to measure the 2.5-/xm contin- 
uum. (All CGAS spectra are normalized to 1.0 at 2.5- 
,xm ' offset for clarity. Thermal emission is in- 
u<ied. with error bars = I (t). 



Ftc. 8. Seven CGAS spectra of low-albedo asteroids. 
These asteroids display little or no reflectance 
drop; their surfaces are anhydrous. See Table II for 
observation dates. 


Ceres spectra taken over several nights 
show remarkable consistency. Rotational 
variations in the Ceres and Pallas spectra 
are not significant, because on a given night 
we averaged spectra over at least one-third 
of a complete asteroid rotation. (Ceres’ pe- 
riod. 9.08 hr; Pallas: 7.88 hr, from Bums and 
Tedesco [1979]). 

Several features in the Ceres CGAS spec- 
tra are common to other asteroids: 

(/) 3.3-ixm methane artifact . Successful 
averaging of many spectra taken over a wide 
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Fig. 7. Cooled Grating Array spectrometer spectra 
of 2 Pallas taken over several nights in April 1987. 
Pallas shows a prominent 2.5-3.0 /xm absorption, but 
no ice or organic features. 


range of airmass requires an accurate ex- 
tinction correction for both standard star 
and asteroid. Limits on the precision of this 
airmass compensation prevented complete 
removal of the deep 3.3-/xm terrestrial meth- 
ane absorptio n, p roducing a corresponding 
spike in the asteroid reflectance spectrum. 

(2) Thermal emission. Thermal emission 
is apparent in Ceres CGAS spectra at and 
beyond 3.5 //.m (especially compared with 
corrected CVF spectra). This thermal com- 
ponent is noticeable to varying decree in 
other CGAS spectra. 

2 Pallas spectra show similar repeatability 
(Fig. 7). The spectral shape is similar to that 
of CM chondrites. The 2.5-/xm continuum is 
more nearly level, no 3.1-/xm ice feature is 
present, and no organic feature is evident at 
3.4 /xm. Pallas’ H.O band depth d is 0 18 
(± 0 . 01 ). 

The high quality Ceres and Pallas CGAS 
spectra are directly comparable to carbo- 
naceous chondrite reflectance data (Jones 
1988; Jones et al. 1991, in preparation). 
Ceres’ 3-/xm spectrum is unlike that of any 
carbonaceous chondrite, and only distantly 
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resembles the Cl chondrites. Pallas’ spec- 
trum is quite similar to those of CM chon- 
drites, but its H,0 band is not as deep (d = 
0.18 vs 0.35 for CMs). A possible explana- 
tion is that Pallas surface has a significant 
anhydrous component mixed with CM-like 
hydrated silicates (Larson et al. 1983 Jones 
1988, Jones et al. 1990). 

Other Asteroids 

We observed a number of non-C asteroids 
to verify that 3-/im bands are not typical of 
other, anhydrous classes of asteroids. For 
example, 5 Astraea and 16 Psyche (repre- 
senting classes S and M, respectively) show 
no 3-jim reflectance drop (Fig. 8). However, 
two of the “experimental control” aster- 
oids, 55 Pandora and 92 Undina, did have 
definite 3-/xm absorptions (Fig. 9). Both are 
M-class, supposedly anhydrous metallic as- 
teroids. With featureless near-IR spectra, 
neither object was well-categorized under 
the "TRIAD” Tucson Revised Index of As- 
teroid Data system (Zellner 1979). The am- 
biguity arose because visual/near-IR spectra 
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Fig. 10. CGAS spectra of four low-albedo asteroi'.*- 
from 1987 to 1988. See Table II for observation da: 
The thermal emission and noisy spectra make be d 
depth estimates uncertain, but determination of hy- 
drated silicate presence is still possible. 



Fig. 9. CGAS spectra of seven low-albedo asteroids 
displaying 3-^m absorptions of varying strength. 55 
Pandora and 92 Undina are M-ciass objects that never- 
theless have prominent hydrated silicate features. See 
I able II for observation dates. 


from the C, M, E, and U classes cannot be 
discriminated without albedo data. Thoien 
(1984) reclassified both objects as M on the 
basis of their relatively high albedos. Hov. - 
ever, the IRAS albedos for 55 and 92 (0.32 
and 0.20, respectively) were dependent on 
assumptions in the thermal model about the 
surface thermal inertia. If these objects have 
rocky rather than the assumed dusty sur- 
faces, the resulting higher thermal inertia 
would lower the emitted thermal flux and the 
calculated albedo in turn. A lower albedo 
would be consistent with the observed 3- M m 
absorptions, more typical of C asteroids. 
Interestingly, a Soviet polarimetry study of 
M-class asteroids found anomalous polar- 
ization properties for both 55 and 92, and 
concluded that their surfaces were some- 
how atypical and/or had a significant silicate 
component (Lupishko and Belskaya 19891. 
From our observations, those silicates pres- 
ent are certainly hydrated. 

We note here that in our examination of 
individual asteroid spectra (Figs. 6-12) we 
found no trace of any 3.4-/xm organic ab- 
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sorption. We are thus unable to confirm the 
direct detection of organic kerogen on any 
asteroid (Cruikshank and Brown 1987). 

III. HYDRATED SILICATE DISTRIBUTION IN 
THE ASTEROID BELT 

We consider a sample of 38 3-jtxm spectra 
from low-albedo asteroids, derived from 
-VF, CGAS, and earlier spectrophotomet- 
nc data (Lebofsky 1980, Feierberg et al. 
1985a,b). The measured 3-/xm band depths 
when compared to asteroid physical and or- 
bital characteristics, offer new insight into 

the original H,0 distribution in the outer 
belt. 

First, we examined the relationship be- 
tween asteroid class and the presence of a 
-Atm absorption. The sample includes all 
>ow-aIbedo classes: C, P, D, (and the C sub- 
classes G, T, F, and B). Overall, 66% of 
the C and C subclass objects show a 3-Aim 
absorption diagnostic of hydrated silicates 
(Fig. 13). The percentage differs for each 
class (C — 68%, G— 80%, T — 50% 

F— 66%, and B— 100%), but clearly a ma- 
jority of Cs have hydrated silicate surface 
components. By contrast, of the seven Ps 
and Ds in the sample, only one shows a 3- 
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V rapidly nsmg thermal component at 

and beyond 3 M m precludes searches for absorption 
features other than water of hydration. Removal of 
the modeled thermal component revealed Thisbe as an 
anhydrous asteroid. 


Mm band, though 324 Bamberga’s feature is 
very weak. Recent observations by Lebof- 
sky et al. (1990) reinforce this contrast: of 
an additional P and two Ds observed in 1988 
none exhibited a 3-/tm band. To repeat, only 
1 of 10 Ps and Ds observed to date is hy- 
drated. We conclude that low-albedo aster- 
oids fall into two predominant groups, char- 
acterized by different hydration states: the 
carbonaceous chondrite-like Cs, and the ul- 
tracarbonaceous Ps and Ds. 



ASTEROlO C_ ASS 


F 'C- 1 1- CGAS spectrum of 423 Diotima (C class) 
rom Oct. 1987. The thermal contribution to the flux at 
mese wavelengths is quite strong, but removal of that 
component at 2.9 M m still permits a qualitative search 

hydrated" Mm band ’ D '° t0ma does not a PP ear t0 


Fig. 13. Observed distribution of the 3- M m hydrated 
silicate band over the various low-albedo asteroid 
classes. About two-thirds of the C and C subclasses 
are hydrated while the P and D classes are nearly 
anhydrous. If the Lebofsky et al. (1990) observations 
are included, the contrast is even more striking. 
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(A~) 

Fig. 14. Observed low-albedo asteroid reflectance at 
3-^m vs solar distance. There is a highly significant 
statistical correlation between increasing reflectance 
and increasing distance, suggesting that hydrated sili- 
cates decline in abundance with distance. Data are from 
this work, Lebofsky (1980), and Feierberg et al. 
(1985a. b). 


A quantitative measure of this class differ- 
ence can be seen in a comparison of each 
asteroid’s 2.9-/u.m reflectance value level (or 
band depth d) to orbital semi-major axis (so- 
lar distance). Because most Ps and Ds reside 
beyond 3 AU in the outer belt, a class-wide 
absence of hydrated silicates there would be 
recognizable in the form of higher 2.9-^m 
reflectance values. Figure 14 confirms this 
trend: in our sample, there are no hydrated 
asteroids beyond 3.4 AU (the latest 3 Ps 
and Ds from Lebofsky et al. [1990] are all 
anhydrous and orbit beyond 3.9 AU). 

A* weighted linear least-squares fit to the 
sample (Fig. 14) produces a positive slope 
and a correlation coefficient of R = 0.43. 
For N = 38, the probability that an uncorre- 
lated set of data points would yield R = 0.43 
is <1%; the probability that the data are 
correlated is >99%. These data indicate a 
statistically significant decline in hydrated 
silicates (and surface H 2 0 content) outward 
from the C portion of the belt (Jones et al. 
1988b). They also imply that the supposedly 
volatile-rich Ps and Ds are the most water- 
poor objects (spectrally) in the outer belt! 

The relationship between solar distance 
and 2.9-p.m reflectance does not appear to 
be a simple linear one (Fig. 14). The rela- 
tively large errors in the reflectance data 



Fig. 15. Observed low-albedo asteroid reflectance a: 
3-ptm vs geometric albedo. Reflectance declines with 
increasing albedo. The correlation coefficient. R 
0.49, indicates that the positive statistical correlation 
between albedo and hydrated silicate abundance i- 
highly significant. 

combine to produce such a poor \ l val - 
that we use the linear model here only to 
show correlation. A more complicated pro- 
cess must have been at work to create the 
observed H : 0 distribution. 

Plots of 3-/xm reflectance vs diameter, al- 
bedo, eccentricity e, and inclination i are 
shown in Figs. 15-18. Band depth appear 
unrelated to diameter, e, and i, argui 
against <£’-collisional origin for the prese 
population of anhydrous and hydrated aste: 
oids. If, as Feierberg et al. (1985a) sug 
gested", collisions had shattered mildl; 
heated asteroids into fragments of hydrate^ 



DIAMETER (KM) 

Fig. 16. Observed low-albedo asteroid reflectance 

3-^m vs diameter. Reflectance and diameter also 
pear strongly correlated, but the linear fit is hea^ 
influenced by Pallas and Ceres at diameters >500 k 
No clear relationship between diameter and hydra: 
silicate abundance can be demonstrated from the 
data. 
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Fig. 17. Observed low-albedo asteroid reflectance 
at 3-ptm vs orbital inclination. The parameters appear 
statistically uncorrelated. Higher inclinations may be 
a result of asteroid collisions: the data suggest that 
collisions are unimportant in controlling hydrated sili- 
cate abundance. 


mantle and anhydrous core remnants (dehy- 
drated from an original phyllosilicate com- 
position), then the wetter asteroids (the 
proto-asteroid outer layers) should be in 
generally more inclined or eccentric orbits. 
Such a trend is not evident. Further, one 
anhydrous asteroid, 52 Europa. is 312 km in 
diameter; it -is difficult to imagine the total 
dehydration of such a large core and the 
subsequent loss by impact of all hydrated 
surface components, while still retaining the 
typical C-class spectrum and low albedo. 
The hydration state of the outer belt must 
be largely a function of original composition 



Fig. 18. Observed low-albedo asteroid reflectance at 
3-/xm vs orbital eccentricity. As in Fig. 17, there is 
no correlation. Higher eccentricities may be due to 
collisions as well as gravitational perturbations: these 
data suggest that collisional history does not control 
hydrated silicate abundance. 


and thermal evolution rather than collisional 
history. 

An apparent inverse relationship between 
2.9-/xm reflectance and albedo supports this 
conclusion (Fig. 15). Here the probability 
that these two parameters are correlated 
(R = 0.49, A’ = 38) is >99%. Deeper H ; 0 
bands tend to occur on objects with higher 
albedo, perhaps due to aqueous emplace- 
ment of hydrated salts, carbonates, or sul- 
fates, all light-colored meteoritic phases. 

The 2-pim reflectance data help to inte- 
grate the low-albedo asteroids into the con- 
text of the entire asteroid belt: beyond the 
differentiated inner belt asteroids (Gaffey 
1988). most of the Cs display traces of mod- 
erate internal heating and subsequent aque- 
ous alteration. The Ps and Ds, however, 
despite a supposedly volatile-rich heritage, 
display no surface H : 0 (Lebofsky et al. 
1990). Our observations support a model for 
the formation of low-albedo asteroids that 
explains this seeming paradox. 

IV. ORIGINAL OUTER BELT ASTEROID 
COMPOSITION 

The decline in hydrated silicate abun- 
dance with increasing heliocentric distance 
runs counter to the expected volatile en- 
hancement in the ‘'primitive” P and D 
classes. This paradox must be considered in 
the context of the original composition these 
objects probably acquired during accretion 
from the solar nebula. 

A starting point (and only that) for esti- 
mating the probable raw materials available 
for asteroid accretion can be found in the 
chemical equilibrium condensation models 
of Lewis (1972) and Grossman and Larimer 
(1974). In these models proto-asteroids ac- 
crete from solid grains condensed under the 
prevailing nebular conditions in the vicinity 
of the asteroid belt. In an ideal solar nebula 
model, condensed solids at all distances 
equilibrate with the surrounding gas, with 
conditions near the Sun favoring the closest 
approach to chemical equilibrium between 
the two phases. Abundant meteorite evi- 
dence suggests that thermodynamic equilib- 
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riuffi was not attained in the meteorite (and 
presumably asteroid) formation zone. 
(Wood and Chang 1985). 

We are interested in particular in the ef- 
fects of disequilibrium on the production of 
asteroidal hydrated silicates. These low- 
temperature minerals are expected products 
of equilibrium nebular condensation, and 
are stable against water loss under asteroi- 
dal surface conditions. However, the idea of 
hydrated silicates as “primordial” nebular 
material has fallen from favor. The present 
consensus among meteoriticists is that hy- 
drated silicates are of secondary origin (Ker- 
ridge and Bunch 1979, McSween 1979, 
Bunch and Chang 1980, Fredriksson and 
Kerridge 1988). In addition, there is strong 
theoretical evidence arguing against a nebu- 
lar origin for hydrated silicates (Lewis and 
Prinn 1984, Prinn and Fegley 1987, 1988, 
Fegley 1988). 

The fundamental implication for our work 
is that asteroidal hydrated silicates, like 
those in meteorites, must also be alteration 
products. Two corollaries are readily appar- 
ent: (1) There was sufficient water available 
on parent asteroids to produce hydrated sili- 
cates abundant enough to survive impact 
processing and be detectable today; and (2) 
an early heat source must have warmed 
these objects enough to drive an aqueous 
alteration episode. 

Ice as the Asteroidal Water Component 

Without hydrated silicates as an accretion 
feedstock, the C asteroids (and the Cl and 
CM parent bodies) must have originally ac- 
quired their present H ; 0 in the form of ice. 
Bunch and Chang (1980) discussed the pos- 
sibility of ice as the original asteroidal water 
reservoir, but did not propose an explicit 
origin (see also Nozette and Wilkening 
1982). Ice condensation at 3 AU is not an 
absolute requirement for supplying the C 
region with water; there was undoubtedly 
some inward transport of material from Ju- 
piter (5.2 AU) during its formation. (Ice 
must have condensed at Jupiter’s distance 
to have survived to the present on the Gali- 


lean satellites [Lebofsky 1980]). It may have 
migrated inward via turbulent transport of 
dust and ice grains (Prinn 1989), scattering 
of icy Jupiter planetesimals, or via extensi' n 
of the outer belt accretion (or “feeding ) 
zone toward Jupiter, where ice had already 
condensed. If accretion were rapid and the 
material sufficiently icy. significant quanti- 
ties of ice could be incorporated into early 
asteroids and protected from sublimation. 

Surface ices on asteroids do not have a 
long life expectancy (Lebofsky 1980). How- 
ever, any ice incorporated into asteroid in - 
riors during accretion should have survived 
for the few million years required (at most) 
for aqueous alteration (DuFresne and An- 
ders 1962, Kerridge and Bunch 1979). Water 
ice is present today on at least one aster- 
oid— 1 Ceres (Lebofsky et al. 1981, this 
work). Given the high equilibrium tempera- 
tures of low-albedo asteroids at 3-4 A T J 
(-200 K), only large ones like Ceres, v. i 
a generous H 2 0 inventory, may be able :o 
support surface ices. But near-surface ice 
may persist if initial sublimation leaves be- 
hind a lag deposit similar to those proposed 
for comets (Brin and Mendis 1979. Houpis 
and Mendis 1981, Fanale and Salvail 1984). 
Such a lag deposit would probably make 
even near-surface ices spectrally invisible. 

Low-Albedo Asteroid Composition 

Our 3-ptm reflectance data draw a picture 
of hydrated silicate abundance that declines 
outward through the C. P, and D classes 
(Jones etal. 1988a,b. Jones 1988, this work). 
Hydrated silicates are conspicuous by their 
absdnce on the few observed Ps and Ds 
(Lebofsky et al. 1990). The observations 
contradict the proposed trend toward prin..- 
tive compositions (H : 0 and organics) in the 
outer belt (Gradie and Veverka 1980). We 
resolve this apparent dilemma by proposing 
that the hydrated silicate abundance does 
not reflect a proto-asteroid’s original water 
content, but only the extent to which aque- 
ous alteration occurred there. The original 
outer belt asteroid composition was neces- 
sarily a mixture of anhydrous silicates, w v 
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ter ice, and relic interstellar organic mate- 
rial. The absence of hydrated silicates (on a 

P, or D object) implies either (1) a vola- 
ie-poor parent body, or (2) early interior 
temperatures and pressures too loss to sup- 
port an alteration episode. In the latter case 
the asteroidal silicates svould retain their 
original anhydrous character, resembling 
the olivine-rich matrix of the CV chondrites 
(McSsveen 1987). 

V. NEW SUPPORT FOR ASTEROID 
ELECTRICAL INDUCTION HEATING 

The observed decline in hydrated silicates 
with heliocentric distance over some 3 AU 
is inconsistent with :6 Al heating. Instead, 
the range of heating events we observe in 
the belt seems to fit much more closely with 
a wide-ranging episode of solar wind induc- 
tion heating (Sonett et al. 1968, Lebofsky et 
al. 1989, Herbert 1989). 

The C asteroids as a class do not represent 
an unaltered mineralogical suite: some two- 
thirds of our sample have experienced sec- 
ondary aqueous alteration to some degree. 
While hydrated and (apparently) anhydrous 
asteroids are well-mixed inside =3 AU, the 
decline in hydrated silicate abundance with 
distance seems to be real. Oxygen isotopic 
evidence from meteoritic hydrated silicates 
suggests parent bodies achieved tempera- 
tures of 0-20°C (Clayton and Mayeda 1984). 
If we assume hydrated Cs once experienced 
interior temperatures at least that high, then 
the hydrated silicate distribution pictured 
in Fig. 14 corresponds to a crude 300 K 
isotherm. For example, the fraction of aster- 
oids heated to that temperature is close to 
unity at 2.5 AU, and near zero in the P region 
near 3.5 AU. This profile strongly echoes 
me apparent decline in differentiated aster- 
oids from 2 to 4 AU (Gaffey 1988). The radial 
dependence common to Gaffey’s igneous 
“isotherm” and our aqueous alteration pro- 
file in the C region leave little doubt that 
the heating intensity also declined with solar 
distance. 

That dependence is a key argument 
against the action of :6 AI as a heat source in 


the outer belt. The decay of this short-lived 
radionuclide has been proposed as the cause 
of both igneous differentiation and mild 
aqueous alteration on the meteorite parent 
bodies (Lee 1979). There is no doubt as to 
its heating potential. For example, if 26 A1 
was present at the concentrations implied 
by 26 Mg excesses in the Allende CV chon- 
drite, it would have melted bodies of radius 
> ~ 10 km (Lebofksy et al. 1989). That effi- 
ciency is a drawback among primitive, low- 
albedo asteroids, however. Homogeneous 
distribution of 26 Al from 2 to 5 AU would 
have led to wholesale melting of asteroids 
with no radial dependence. Both the succes- 
sion of outer belt asteroid classes (C-P-D) 
and the decline in hydrated silicate distribu- 
tion prove that this did not occur. 

It is a well-known failing of 26 Al that it is 
most efficient at melting large bodies, yet 
the two biggest asteroids. Ceres and Pallas, 
show no evidence of anything beyond low- 
temperature aqueous alteration. The ob- 
served hydrated silicate distribution (Fig. 
14) requires even more capricious action by 
26 Al, because at small scales in the C region 
both heated and apparently unaltered aster- 
oids are thoroughly mixed. 26 Al‘s heating 
characteristics are inconsistent with both a 
general decline in aqueous alteration over 
a 3 AU span and starkly different thermal 
outcomes on a local scale. 

Electrical induction heating (wherein a 
vigorous early solar wind deposits kinetic 
energy in a planetary interior via resistance 
to induced currents) seems to explain both 
the overall belt compositional structure and 
our 3-jU.m hydrated silicate data. Applied to 
small bodies, this mechanism can account 
for the igneous differentiation of inner belt 
asteroids, the apparent decline in hydrated 
silicate abundance with distance, and the 
variety of C-class hydration states (Lebof- 
sky et al. 1989, Jones et al. 1988b, Herbert 
1989). 

Asteroid induction heating is particularly 
adept at explaining the production of hy- 
drated and anhydrous surfaces over short 
distances in the belt. For example, anhy- 
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TABLE IV 


Heat Source Comparison 


Observed 

: ‘AI 

Induction 

characteristics 

heating 

heating 

Beit-wide variation 

Yes 

Yes 

(melted vs unaltered) 
Intensity decline with 

7 

Yes 

distance 

Small-scale heterogeneity 

7 

Yes 

(0.1 AU) 

Ineffective at large 

No 

Yes 

diameters 




drous C asteroid 52 Europa (r = 156 km, 
albedo < 0.06) has a mean solar distance of 
3.1 AU, greater than that of many hydrated 
asteroids, so it probably formed with at least 
as much ice as its sunward companions. Col- 
lisions were probably incapable of dehydrat- 
ing 52 Europa while leaving it with a typical 
C albedo and spectrum. The present anhy- 
drous surface must be a consequence of 
minimal heating. Could a relatively large as- 
teroid like 52 Europa avoid heating and 
aqueous alteration, unlike many neighbor- 
ing Cs? 

Such different thermal outcomes are plau- 
sible, and stem from the known range of 
bulk compositions among the carbonaceous 
chondrites. While similar in bulk chemistry, 
their measured electrical conductivities 
span over 6 orders of magnitude (Brecher 
et al. 1975). Similar variation among the C 
asteroid precursors would have left them 
wildly different in their susceptibility to 
electrical induction heating. Some bodies 
would have been warmed easily, while slight 
compositional variations would have left a 
neighboring object cold. Such selectivity 
would explain qualitatively the variety of 
hydration states within the C region of the 
belt. The more general decline in hydrated 
silicate abundance was probably caused 
more simply by the reduced intensity of in- 
duction heating farther from the Sun. 

We compare 26 Al and induction heating as 
suitable asteroid heat sources in Table IV. 


One can concede that 26 A1 may have b 
more abundant in one region of the aster 
belt, causing widespread heating over, s. 

1 AU. But it is hard to see how 26 Al co 
produce such different thermal outcon 
over <0.1 AU (Fig. 14), or cause a grad 
decline in heating intensity with distan 
Moreover, 26 Al should be most effective 
larger asteroid diameters, a prediction 
supported by our 3-ptm observations ( 
drated silicate abundance is not stron 
correlated with diameter; Fig. 16). In c 
trast, electrical induction heating satisfies 
of the observed characteristics of the m 
to outer belt heat source (Bell 1986, Jo: 
1988). 

VI. OUTER BELT ASTEROID ORIGIN 
AND EVOLUTION 

Our spectra of outer belt asteroids . 
consistent with their early accretion fr 
planetesimals made of anhydrous silicat 
water ice, and organic material. The prop 
tions of those ingredients varied with he 
centric distance, with ice and organics m 
plentiful farther from the Sun. Only the 
dial dependence of the early induction he 
ing episode(s) drove the C, P, and D clas 
down different evolutionary paths: the m 
eratfely heated Cs lost nearly all their wa 
ice, the meltwater supporting the com 
sion of surface minerals into hydrated s 
.cates. To be complete, our scenario m 
also explain the anhydrous surfaces of ' 
supposedly volatile-rich P and D astero: 
which do not exhibit a 3-/xm H ; 0 absorpt 
(Lebofsky et al. 1990). 

Why Do the Ps and Ds Appear Dry ? 

Scenarios that propose mineralogical 
teration of the Ps and Ds to produce an? 
drous surfaces do not square well with av 
able evidence. Heating intense enough 
devolatilize a hydrated silicate object rai 
a host of objections: externally driven de 
dration of just the surface would reqi 
both an improbable original composit 
(phyllosilicates) and, because surface Ia\ 
cool most efficiently, an intense heat sour 


WATER ON LOW-ALBEDO ASTEROIDS 


187 


Neither condition is consistent with the mild 
modification we see among the Cs. Dehy- 
dration by a strong interior heat source is 
even more unlikely, given the great heat sink 
provided by H 2 0 and other volatiles (Lewis 
and Prinn 1984); such a thorough outgassing 
would also result in extensive spectral and 
‘bedo changes. The possibility that organic 
kerogen is masking the 3-^tm H 2 0 absorp- 
tion has been proposed, but laboratory re- 
sults show that suppressing such a strong 
feature is difficult: adding organic material 
lowers the albedo to some limiting value 
0.04), but the water band remains visible 
(Jones 1988, Jones et ai 1989, 1991 , in prep- 
aration). For example, the 3-/xm absorption 
quite strong in spectra of Cl and CM chon- 
drites, despite the intimate mixing of hy- 
drated silicates and organics in the fine- 
grained matrix. The Ps and Ds must have 
trulv anhydrous surfaces. 

Lebofsky et ai (1990) specifically discuss 
the absence of 3-/xm features among the Ps 
and Ds, but a moment spent placing these 
objects in the context of the other low- 
albedo asteroids is worthwhile. The most 
reasonable explanation for the missing H 2 0 
band is incorporation of water in the form 
of ice (Bunch and Chang 1980). The proto- 
Ps and Ds may have been loose agglomera- 
tions of ice and anhydrous silicates, whose 
outer layers quickly lost their ice to sublima- 
tion or gentle solar insolation, while their 
interiors never achieved temperatures high 
enough to drive H 2 0 vapor or fluid to the 
surface (Fig. 19). Left unscathed by electri- 
cal induction heating in the outer belt, the 
Ps and Ds lost their near-surface ice to subli- 
mation over billions of years of exposure 
to vacuum and impacts. A lag deposit of 
anhydrous silicates and carbonaceous mate- 
rial blanketing the surface, unaltered by liq- 
uid water, would explain the lack of a 3-^ 
band, the low albedos, and reddish near-IR 
>pectra reminiscent of comet nuclei (Hart- 
mann et al. 1987). Spectrally, this lag de- 
posit is organic-rich, but H 2 0-poor. No ob- 
servational evidence would preclude the Ps 
and Ds from retaining some interior ice; 



Fig. 19. Proposed scenario for producing low-albedo 
asteroids. Accretion of ice and anhydrous silicates from 
the solar nebula, combined with relic interstellar or- 
ganic kerogen. formed very dark, volatile-rich aster- 
oids, of which the Ps and Ds may be representative. 
Selective induction heating (dependent on solar dis- 
tance, electrical conductivity) produces aqueous alter- 
ation at the surface and perhaps even melting in the 
interior, giving rise to the carbonaceous chondrite par- 
ent bodies and the present C class. Some Cs are unaf- 
fected by heating and may preserve interior ices, along 
with the Ps and Ds. The P and D classes should preserve 
nearly unaltered nebular solids and volatiles. 


these asteroids would resemble comets in 
some aspects of structure and composition 
(Hartmann et al. 1987). Judging from the ice/ 
rock ratios of the Galilean satellites, the D 
asteroids (Trojans in particular) may have 
had a substantial (= 50%) ice component, 
blurring the distinction between outer belt 
asteroids and comets. 

The Outer Belt Imaged at 3 nm 

Our 3-fim reflectance data suggest the fol- 
lowing interpretation of the outer asteroid 
belt (Fig. 20). The proto-asteroids beyond 
2.5 AU were composed originally of anhy- 
drous silicates, H 2 0 ice, and relic interstel- 
lar organic material. Electrical induction 
heating differentiated the inner part of the 
belt (inside 2.5 AU), while the C region was 
subjected to mild aqueous alteration. Be- 
yond 3.5 AU, induction heating was largely 
ineffective, leaving the Ps and Ds unaltered. 
Those asteroids are the best candidates for 
preservation of the “primitive” belt compo- 
sition, including volatiles and the original 
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anhydrous silicate mineralogy. At the top 
of Fig. 20 we show how the percentage of 
“melted” asteroids (differentiated S types) 
declines with heliocentric distance, declin- 
ing to near zero beyond 3.5 AU (Gaffey 
1988). We show qualitatively how the per- 
centage of aqueously altered Cs also dimin- 
ishes markedly with distance past 3 AU. 
Together with Gaffey's 1988 S asteroid 
work, our 3-ju.m hydrated silicate observa- 
tions of low-albedo asteroids strongly sup- 
port Bell’s (1986) explanation for the present 
radial structure of the belt. 

Primitive Asteroids and 

Carbonaceous Chondrites 

The P and D asteroids should represent 
the belt’s least altered mineral assemblages. 
Meteorites appear not to have sampled the 
Ps and Ds, but McSween (1987) proposes 
that the olivine-rich, relatively unaltered CV 



Fig. 20. The 3-pim reflectance of the outer belt aster- 
oids in our survey permits us to trace the hydrated 
silicate profile, where asteroids were heated to roughly 
300 K during the belt-wide heating event. The decline 
in hydrated silicates with heliocentric distance joins the 
radial decrease in the abundance of “melted" asteroids 
in strongly supporting electrical induction heating as 
the mechanism for early asteroid heating. Differenti- 
ated asteroid trend is from Gaffey (1988); the radial 
succession of asteroid classes is redrawn from Gradie 
and Tedesco (1982). 


chondrite matrices are suitable precursors 
for the hydrated CI/CM chondrites; ti _ 
same CV material is a good candidate fur 
the P and D anhydrous silicate component 
as well. Across a span of 2 AU in the belt, we 
would expect a range of primordial asteroid 
compositions, with anhydrous Cs. Ps, and 
Ds all preserving their original mineralogy. 
We note that Bell et al. (1988) nominated 
221 Eos (C/S class) as the CV/CO parent 
body. Its hydration state is unknown, b t 
Eos’ semi-major axis (— 3 AU) would be 
consistent with a low-hydrated silicate 
abundance and an anhydrous, olivine-rich 
surface. 

In contrast, the hydrated C, B. F, G, and 
T asteroids represent proto-asteroidal mate- 
rial modified by mild induction heating. 
These classes probably include the parent 
bodies of the Cl and CM chondrites. 

Our 3-pi.m data fit very well into Bell’s 
interpretive framework for the mineralogi- 
cal evolution of the entire asteroid belt (Bell 
1986, Bell et al. 1988). He proposed three 
“superclasses” that relate taxonomy to the 
observed radial compositional structure. 
The superclasses (igneous, metamorphic, 
and primitive) represent the effects of the 
heliocentrically varying electroinductrie 
heating episode. Igneous types dominate the 
belt sunward of 2.7 AU, the metamorphic 
group occupies a narrow zone centered 
around 3.2 AU, and the primitive class dom- 
inates outside 3.4 AU. 

We would modify Bell’s framework only 
by noting that P and D mineralogy is defi- 
nitely not dominated by “clays” (a common 
assumption prior to the 3-ju.m survey) and 
that the Cs are not entirely a primitive class. 
They are chemically primitive, but that label 
perpetuates the image of carbonaceous 
chondrites as totally unaltered, when we 
know that their hydrated silicate mineralogy 
is secondary. 

VII. FUTURE WORK AND CONCLUSIONS 

The unfortunate failure in 1989 of both 
Soviet Phobos spacecraft eliminated near- 
term prospects for a detailed look at an ob- 
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;ct apparently similar to those observed in 
our study. Phobos has a visual/near-IR spec- 
trum much like 1 Ceres (Pang et al. 1978), 
and may be a captured main belt fragment; 
on the basis of the abundance of hydrated Cs 
in our 3-fj.m survey, discovery of hydrated 
silicates on Phobos would not be surprising. 
Linear crater chains associated with the sat- 
ellite’s parallel grooved fracture planes sug- 
est that it once outgassed water after a ma- 
jor impact (Veverka and Thomas 1979, 
Hartmann 1980). On the other hand, recent 
3-/xm observations of Deimos by Bell et al. 
(1989) suggest that it is much less hydrated 
than the “average” hydrated C spectra we 
present here. On the basis of 3-/i.m and 
shorter wavelength reflectance data, Bell 
and co-workers conclude that Deimos may 
'e an anhydrous P asteroid fragment from 
,.ie Cybele region just outside the main belt. 

Because of their primitive composition 
and possible links to comets, a rendezvous 
with a low-albedo asteroid should be a high 
NASA priority, particularly in light of the 
Phobos failures. While Galileo and later the 
CRAF and Cassini missions will return valu- 
able asteroid data, a great deal could be 
learned cheaply and quickly about low-al- 
bedo objects by launching a small, relatively 
inexpensive spacecraft to a rendezvous with 
an Earth-crossing C (Cole et al. 1990). 
Ground-based 3-yam spectrophotometry 
could identify an appropriate hydrated or 
anhydrous target; a conservative estimate 
indicates that dozens are more accessible 
than the lunar surface in terms of the total 
velocity change required (Lewis and Lewis 
1987). The expected compositional and min- 
eralogical data returned would test the ther- 
mal evolution scenario presented here, nar- 
row the search for carbonaceous chondrite 
parent bodies, and enhance any ground- 
based resource survey of the near-Earth 
population. 

j SUMMARY 

? (1) We obtained 3-yum spectra of 19 low- 

T albedo and 5 S- or M-cIass asteroids. The 

\ IRTF cooled grating array spectrometer 


proved efficient at recording the 3-yu.m re- 
flectance of moderately bright, low-albedo 
asteroids. 

(2) We confirmed a 3.1-/am H : 0 ice fea- 
ture in the spectrum of 1 Ceres. However, 
no 3.4-ju.m organic features were visible in 
even the best asteroid spectra. Two M-class 
objects displayed 3-jixm H ; 0 bands that indi- 
cated a hydrated surface component on 
these “metal” asteroids. 

(3) 1 Ceres has no spectral match among 
the carbonaceous chondrites. Its 3-/xm band 
depth suggests a surface H ; 0 content of 
about 5 wt%. The 2 Pallas spectrum closely 
resembles that of CM chondrites, indicating 
that Pallas’ water content is roughly 2 wt%. 

(4) Two-thirds of the C, G. T. F, and B 
asteroids in the 3-yxm sample display a hy- 
drated silicate absorption, indicating that 
the Cs (and C subclasses) underwent a wide- 
spread, though mild, heating and aqueous 
alteration episode. Only l of 10 P and D 
spectra show even a weak H ; 0 absorption, 
indicating a major difference in the hydra- 
tion states of the C and P/D classes. The 
latter appear to have escaped extensive 
heating and aqueous alteration. 

(5) Hydrated silicates gradually decline in 
abundance through the C region of the outer 
belt, disappearing entirely among the P and 
D classes. The presence of hydrated sili- 
cates correlates positively with albedo, but 
is unrelated to eccentricity, inclination, or 
diameter. The water distribution among the 
C, P, and D classes is consistent with an 
original composition of anhydrous silicates, 
water ice, and complex organic kerogen. 
Complex carbon compounds and H : 0 were 
probably more abundant in the outer belt. 

(6) The declining abundance of hydrated 
silicates with heliocentric distance, the vari- 
ety of C asteroid hydration states at near- 
identical solar distances, and mild warming 
of even large asteroids is consistent with 
electrical induction heating just after aster- 
oid formation. Anhydrous Cs are not dehy- 
drated core fragments, but primitive bodies 
that escaped induction heating. 

(7) The volatile inventory of the outer belt 
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was dominated by water ice. The inferred 
composition of the Ps and Ds suggests that 
they may still retain interior ices. 
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Solar Nebula Origin for Volatile Gases in Halley’s Comet 1 
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Modeling of physical and chemical processes in solar nebula environments is applied to 
the present data base on the composition of gases Iq Halley r s comet to infer the conditions 
under which that comet formed. Key molecular ratios — CH* to CO% CO^ to CO, and N"Hj 
to — are compared to predictions of solar nebula models. Nebular thermochemistry is 
quantified for a range of solar elemental compositions corresponding to varying degrees 
of depletion of water in the inner part of the nebula. The edge of the chemically active 
region is fixed by rates of chemical reactions in the presence of grains. The mixing of 
molecules from chemically active regions to the outer nebula is parameterized using recent 
studies of this process* The incorporation of volatile molecules into water ice by clathration 
and absorption is then applied to yield predicted ratios of volatile species for comparison 
with the data. We find that, assuming chemistry in the inner nebula is catalyzed by 
reactions on grains* the abundances of the volatile carbon species CH 4 , CO, and COj in 
Halley could have been supplied by the solar nebula. The NH^ abundance in Halley is, 
however, too high to have been derived from the nebula; an alternative source must be 
invoked. We suggest that the surrounding giant molecular cloud is a possible source for 
ammonia. This requires that nebular mixing not be 100% efficient, also consider the 
heating and processing of grains entering the nebula by gas drag, c nw Academic Pr*ss. lac. 


INTRODUCTION 

The origin of the solar system is now rec- 
ognized as part of the larger issue of the 
phenomenon of star formation in the Gal- 
axy. Consequently, the chemistry of the 
disk which was to become the planets (that 
is, the solar nebula) must have been inti- 
mately tied up in the chemistry of the sur- 
rounding giant molecular cloud which was 
the likely supply for the nebular gas and 
dust. While meteorites are regarded as pri- 
marily products of solar nebula processes, 
comets are variously considered to be pris- 
tine remanants of interstellar processes, mo- 
lecular cloud material partially altered in the 
solar nebula, or products entirely of solar 
nebula condensation and chemistry. Given 

1 Contribution 89-32, Theoretical Astrophysics Pro- 
gram, University of Arizona. 


the diversity of materials thought or kno" " 
to be contained in comets, ranging fr.m 
highly refractory silicates through organics 
of moderate volatility to ices and gases of 
high volatility, the middle view must be the 
appropriate one: comets contain records or 
processes stretching from those in the mo- 
lecular cloud (and perhaps even prior ep- 
ochs), through infail and entry into the solar 
nebula, to the solar nebula itself. 

To understand the history of comet r\ *- 
mation, distinct epochs must be isolated :> 
selecting materials within a fixed range ot 
volatility which may say something about a 
particular set of events in one (or a small 
subset of) primordial environments. The 
most volatile species in comets, namely the 
small molecules with low polarizability 
(methane, carbon monoxide, noble gaser. 
likely reflect the final stages of chemist’} 
and physics occurring in the solar nebu;a. 
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TABLE I 

Abundances in Comet Halley 
Used in this Study - 


CO 

0.05-0. 1 

N, 

<0.05 

ch 4 

0.01-0.04 

NHj 

0.003-0.01 

CO; 

0.02-0.04 

H : CO 

0 . 01 - 0.1 

HCN 

0.00 1 


a Molecular abundances are 
mole fractions relative to water. 
Inequality indicates an upper limit 
(nondetectionj. 


although some of these materials may have 
had their origin in the surrounding molecular 
cloud. The comprehensive ground- and 
s~ice*based observations of Halley's comet 
i, 1985 and 1986 provided the first opportu- 
nity to assemble an inventory of gases con- 
tained in comets, and interpretation of this 
inventory in terms of cometarv formation 
environments is now appropriate. Table I 
from Lunine (1989) lists the abundances of 
key molecular species in Halley’s comet rel- 
ative to water. The table is generally consis- 
tent with that in Weaver (1989). The inter- 
p/etation of the ground-based and 
spacecraft data is an involved and contro- 
versial process; the aforementioned papers 
discuss the inventory in terms of observa- 
tional uncertainties in some detail. We 
therefore do not repeat the discussion here; 
original references for the data can also be 
found in those two papers. 

In the present effort we consider a range 
■"f processes in the solar nebula which may 
have determined the budget of volatile gases 
in the outer solar system. The model is tuned 
to, and compared with, the budget of gases 
shown in Table I. While we find that many 
of the carbon-based molecules have abun- 
dances which can be fitted by a plausible 
solar nebula model (with certain modifica- 
tions and caveats), other species, such as 
formaldehyde and ammonia, decidedly can- 
not. They indicate that the chemical history 


of at least a subset of volatiles stretches 
beyond the solar nebula to other sources. 
We consider these sources toward the end 
of the paper. 

Prinn and Fegley (1989) were the first to 
consider the gas budget of Halley in terms 
of solar nebula models. They concluded that 
substantial mixing of material from the vi- 
cinity of giant planets was required to ex- 
plain the data. Why then the present analy- 
sis? The new features of this study include 
(1) consideration of the partitioning of gases 
between the nebula and ice phases, which 
alters the inventory of gases trapped in the 
ice, and hence the comet, compared to what 
was in the nebula; (2) a model for the ther- 
mochemistry of molecular species in the 
chemically active zone of the nebula, which 
considers the elemental carbon-to-oxygen 
ratio to be a varible, based on the model of 
diffusive redistribution of water by Steven- 
son and Lunine (1988): (3) consideration of 
different efficiencies of mixing of gases in 
the solar nebula, which affects the degree of 
chemical coupling between the outer and 
inner regions of the nebula; and (4) calcula- 
tion of the thermal history of icy grains fall- 
ing into the solar nebula and subjected to 
drag heating. 

In what follow s we assume that the com- 
position of Comet Halley is reflective at 
least in part of the physical and chemical 
processes associated with the origin of the 
solar system and planets. A different view- 
point has been put forth, that Halley is ge- 
netically completely unrelated to the solar 
system but was captured well after the ep- 
och of planetary formation. Support for this 
model has come in the form of a nonterres- 
trial ,3 C/ ,2 C ratio (Wyckoff er a/., 1989). 
While we acknowledge the possibility that 
Halley is an “oddball” (the prevailing view 
is that most comets could not have been 
captured in this fashion), we also point out 
that understanding of the carbon isotope 
data is far from complete. Until and unless 
this issue is resolved, we make the assump- 
tion that formation of Halley was coeval 
and colocated with that of the solar system. 
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Further observations ot other comets and 
the upcoming CRAF mission may help re- 
solve this critical issue. 

Section 2 describes the modeling of nebu- 
lar physical and chemical processes; section 
3 zives the model results. Section 4 de- 
scribes the implications of the results tor 
the ability of solar nebula models to explain 
cometary compositions: this section also de- 
scribes models of the heating ot grains tail- 
ing into the nebula. Section 5 concludes the 
paper with prospects for future measure- 
ments and their impact on theory. 

:. NEBULAR PHYSICAL AND 
CHEMICAL PROCESSES 

A. Thermochemical Calculations 

In the inner portion of the nebular disk, 
temperatures are sufficiently high that a 
suite of reactions interchanging CH* with 
CO and CO ; . and NH, with N : . proceed at 
rates such that equilibrium is reached on 
time scales much shorter than the lifetime 
of the nebular disk. Because temperatures 
fall with increasing radial distance from the 
proto-Sun. there must exist a region, for 
each suite of reactions, at which the time to 
achieve equilibrium among molecular spe- 
cies becomes comparable to the lifetime of 
the nebula (or. more stringently, the time 
scale for dynamic mixing). This region fixes 
the ratios of the molecular species for any 
parcels of gas mixed outward from that re- 
gion (i.e., transported to lower temperature 
regions). This concept was first quantified 
for the nebula by Lewis and Prinn ( 1980) and 
Prinn and Fegley (1981). Provided nebular 
mixing is sufficiently vigorous that there is 
significant interchange of gas between the 
cold outer nebula and the hot. chemically 
active inner nebula, this “quench zone” will 
fix the inventory of gas-phase molecular 
species in the entire solar nebula. Prinn and 
Fegley (1989) argue that thermochemistry 
of this sort predominates over chemistry 
driven by other energetic processes, such as 
lightning and nebular flares. We therefore 
assume here that thermochemistry is the pri- 


mary determinant of the molecular inven- 
tory in the inner portion of the solar nebula. 

To model the gas-phase chemistry occur- 
ring in the inner disk, we constructed a com- 
prehensive chemical model which computes 
equilibrium among 50 molecular species of 
the elements H. 0. C. N, and S. Thermody - 
namic data were taken from the J.-.NAF 
bles (1971)- Equilibrium was evaluu. J 
along a hydrogen-helium adiabat such that 
a temperature T of 2000 , K occurs at 10 
bar total pressure and 600°K at 10 bar. 
Based on the nebular model of Cameron 
(1978). 680°K occurs at 1 AU from the proto- 
Sun. Recent models (reviewed in Wood and 
Morfill 1988) indicate that such a tempera- 
ture profile does not reproduce well the .:- 
dial temperature profile in a necula w n 
varyins opacity due to grains, different 
mechanisms of angular momentum trans- 
port. etc. However, the model results are 
not sensitive to variations of pressure at a 
given temperature over a coupie ot oruers 
of magnitude. Our application of the results 
to Hailey is sensitive to the temperature as 
a function of radial distance from the S n. 
through the models of mixing of gas to :.:e 
outer part of the nebula. The reader w *11 
see in Section 3 that it is straightforward to 
assess the effect of such uncertainties. 

For a given set of elemental abundance;, 
and a fixed temperature and total pressure 
in a given region, the partial pressures for 
the various molecular species are calculated 
iteratively. The results are shown in Fig. ' 
The temperatures correspond to the adiat. : 
described above. Consistent with the result: 
of Lewis and Prinn (1980) and Prinn and 
Fealey ( 198 1 ) , CO and N : dominate at higher 
temperature, and CH^ and N; at lower .em- 
perature. The water abundance is tightly 
coupled to the abundance of CO. especially 
for the O/C = 1 case. The CO ; abundance 
is a strong function ot the elemental abu 
dance as well. The nitrogen-bearing specie- 
while sensitive to temperature, are not sen- 
sitive to the oxygen abundance expresses 
by the oxygen-to-carbon ratio. 
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0/C = 1.66 



Fig. 1 . Calculated equilibrium abundances for vola- 
tile species of interest in this paper along a H-He adia- 
bac are shown. The O/C ratio was assumed to be solar. 


B. Reaction Kinetics and the 
Quench Temperature 
A continuing uncertainty which plagues 
the modeling of solar system chemistry con- 
cerns reaction rates. In particular, limita- 
tions in understanding the gas-phase reac- 
tions which interconvert CO and CH 4 have 
been well documented (SPASE 1989). We 
accept the rates used by previous authors 
[reviewed very recently by Prinn and Fegley 
(1989)]. The numbers we give below refer to 
a nebular lifetime (or most sluggish mixing 
time) of 10 13 sec, based upon astrophysical 
evidence on the lifetime of disks around 
young stars (Strom et al . 1989). In the gas 
phase, the CO-CH 4 reaction system 
quenches at 1470°K. the N^-NHj intercon- 
version at 1600°K. and CO-CO : at 830°K. 
Reference to Fig. 1 then gives the abun- 
dance ratios at those temperatures. 

Anticipating the discussion in Section 3, 
these temperatures are so high that very 
little CH 4 or NH } is expected anywhere in 
the solar nebula, except in the high-pressure 
nebulas surrounding the giant planets, 
which are highly reducing environments. 
This has been noted by Prinn and Fegley 
(1989), in connection with their analysis of 
Halley's comet. However, the presence of 
a significant population of refractory grains 


in the inner nebula opens the possibility of 
catalyzed reactions, as first noted by Lewis 
and Prinn (1980) and Prinn and Fegley 
( 1981 ). These reactions effectively lower the 
quench temperature, provided they occur. 

There are currently no experimental data 
at the proper temperatures and hydrogen 
pressures to determine directly the cata- 
lyzed rates of reaction. Vannice (1975) 
comes closest, with high (but not solar) ra- 
tios of hydrogen to CO. The dominant prod- 
uct of the catalysis is C[ hydrocarbons, with 
smaller amounts of higher hydrocarbons. 

On the basis of those experiments, Prinn 
and Fesiey (1989) calculate a quench tem- 
perature of 520 ; K. Similarly, they calculate 
a quench temperature tor iron-catalyzed 
conversion of N : to NH ; , of 480 3 K (termi- 
nated in this case by formation of iron oxides 
in the nebula). Figure 1 show's that while 
still dominates over other nitrogen-bearing 
species, the lower quench temperature pro- 
duces a CH 4 -dominated carbon inventory. 

Prinn and Fegley (1989) have pointed out 
that the conditions under which effective 
grain catalysis may occur are not well under- 
stood. They propose that if grain surfaces 
are not relatively pure (that is. become 
coated with carbonaceous or sulfur-bearing 
material) the catalysis becomes highly inef- 
ficient. They also propose that production 
of heavy-organic material is a more probable 
outcome than that of methane. While the 
experimental data do not exist to answer 
this question definitively one way or the 
other, we note that the available data do 
show positive results in the CO-CH 4 con- 
version, but that the issue of what consti- 
tutes an effective catalysis surface remains 
open. Our model results for methane are 
predicated on efficient catalysis, and the 
reader is reminded of this in later sections. 

C. Diffusive Redistribution of 

Water Vapor 

The decision to let the O-to-C ratio vary 
rests on the realization that nebular trans- 
port processes may act to produce strong 
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regional variations in elemental abundances 
in the solar nebula. Evidence for this is pres- 
ent in chondrites (Rubin et aL 1988), al- 
though much of this effect may be due to 
dust in the nebula and possibly environ- 
ments other than the nebula itself. Steven- 
son and Lunine (1988) modeled the effect 
of transport and condensation processes on 
water. Under the assumption that viscous 
processes transport mass and angular mo- 
mentum, the water distribution through the 
nebula as a function of time was tracked. 
The condensation boundary in the region 3 
to 5 AU acts as a “cold finger,” to actively 
and continually remove water from the gas 
phase. Consequently, the profile of nebular 
water vapor may be strongly altered on a 
time scale comparable to the evolution of 
the disk itself. In the later stages of the neb- 
ula, when material is not being resupplied 
from the surrounding molecular cloud, the 
inner part of the disk could be strongly de- 
pleted of water, and hence have an average 
oxygen-to-carbon ratio of unity. The cosmo- 
chemical implications of this for inner solar 
nebula bodies is potentially important but 
complicated, as the water depletion may 
vary both spatially and temporally [see the 
discussion following Fegley and Prinn 
(1989)]. Figure 2 shows equilibrium abun- 
dances of molecules for the case in which 


O/c = i.o 



600 500 1000 1200 1*00 1500 1800 2000 

Temperature (K) 

Fig. 2. As in Fig. I . but the O/C ratio is unity. Oxygen 
is depleted due to the removal of water vapor to the 
water condensation boundary. CH 4 is more abundant 
than in the solar case and CO and H ; 0 are depleted. 


O/C = I. Methane is now much more abun- 
dant at a given temperature than for the 'ur 
case, and carbon dioxide is depleted. ' ^ 
nitrogen species are not affected. 

D . Nebular Mixing 

The calculations considered above focus 
on active chemistry in the hot inner par: of 
the disk, inward of 1 AU [for hot models of 
the nebula, such as those of Boss (1983). 
active thermal chemistry could occur a- - ar 
out as 3 AU]. If such chemistry affe. d 
the abundances in Halley, which probacly 
formed no closer in than 15-20 AU [see the 
review on cometarv origins by Weissman 
(1985)], radial mixing of material from the 
outer and inner portions of the disk would 
have had to occur. Traditionally, studies of 
nebular thermal chemistry as in Lewis and 
Prinn (1980) have assumed that the in er 
and outer parts of the disk were fully mi .d 
over the lifetime of the nebula, and hence 
the molecular composition was determined 
entirely by that at the radial quench zone. 

Good models relating the transport of 
nebular angular momentum to the mixing of 
trace species do not yet exist; the problem 
probably requires treatment by three-di- 
mensional hydrodynamic codes. Stevens :n 
(1989) considered a semianalyticai mode; A 
a viscous accretion disk. By linearizing the 
stress term associated with angular momen- 
tum transport, it is possible to solve the con- 
tinuity equations for the viscous disk. Ste- 
venson assumes D = kv. where D is the 
diffusion coefficient for mass transport and 
v is the eddy viscosity (diffusion coefficient 
for angular momentum transport). He ar- 
gues that numerical models place the pr. - 
portionality constant k between 1 and 3. 
which severely limits the degree of radial 
mixing in the nebula. Approximately, the 
mixing ratio of a species produced at radial 
distance R p has a mixing ratio at R diluted by 
the factor (R ? IRy ,Zk . For a chemical quench 
zone at 1 AU and a comet formation region 
at roughly 20 AU. the mixing ratio is dilute I 
by 100 for k = 1 and by 5 for * = 3. Wtu: 
is occurring physically is that because the 
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e '.dv and mass diffusion coefficients are as- 
. med proponionai, the total amount of 
mass mixing over the nebular lifetime (de- 
fined by the redistribution of a significant 
fraction of the disk angular momentum) is 
fixed. A higher mass diffusion coefficient is 
accompanied by a higher angular momen- 
tum diffusion coefficient and hence shorter 
nebular lifetime . 

Prinn (1989) argues that because the eddy 
• ress terms are quadratically nonlinear, a 
unearized solution does not capture all of 
the processes involved in mixing. He points 
out that natural systems exhibit nonlinear 
acceleration processes which may produce 
an effective k of 10 or higher. For a k of 10, 
a trace species mixing ratio at 20 AU is over 
60 % that at 1 AU where it is assumed 
formed. In this concept the inner and outer 
-onions of the nebula are nearly fully 
nixed. 

As noted above, the resolution of this 
problem requires a three-dimensional, hy- 
drodynamic investigation of various nebular 
models, in which angular momentum and 
mass are transported not only by viscous 
(local) processes but also by waves. Here 
we consider three cases in which material is 
mixed from 1 AU out to the region 15 to 20 
AU, corresponding to 100%, 10%, and 1% 
mixing. From the scaling given above the 
reader can select others. 

E. Trapping of Volatile Molecular Species 

in Water Ice 

The spacecraft and ground-based investi- 
gations of Halley did not give information 
as to how the various parent molecules were 
locked into the nucleus originally. Direct 
condensation is unlikely for methane and 
more volatile molecules based on solar neb- 
ula temperature-pressure profiles. In the in- 
terstellar cloud, attachment of molecules 
such as N 2 and CO to grains, followed by 
hydrogenation to form H ; 0, CH 4 , and NH 3 , 
has been suggested (Blake et al. 1987). We 
show below that such grains probably evap- 
orated upon entry into the solar nebula. The 
most likely mechanism in the solar nebula 


for getting very volatile species into comets 
is trapping of gases in water ice, either by 
adsorption or by clathration. Trapping on 
media other than water ice is possible, but 
in the outer solar nebula the abundance of 
water is such, and its ability to trap species 
great enough, that we assert it to be most 
important in this process. 

Highly uncertain is the dominant mode of 
trapping. Physical adsorption will occur on 
icy surfaces, but for anything but very small 
ice grains (i.e., less than 0.1 /im), the 
amount of gases trapped on surfaces is 
small, not enough to explain the ratio of 
volatile species to water ice of over 10% 
seen in Table I. Volumetricallv significant 
adsorption (caused either by codeposition 
of water vapor and volatiles or by exposure 
of fresh ice to the gas) can trap sufficient 
volatiles, provided the water ice is amor- 
phous (Mayer and Pletzer 1986). Clathrate 
formation is, roughly speaking, a three-di- 
mensional adsorption process which in- 
volves reorientation of the ice crystal struc- 
ture to produce large voids, capable of 
trapping volatile species such as methane 
and nitrogen. The ratio of trapped species to 
ice is about 13%, consistent with the Halley 
data, but requires physical exposure of all 
of the ice to the gas, since diffusion is very 
slow (Lunine and Stevenson 1985). 

Codeposition is an unlikely way to trap 
volatiles in the nebula, in spite of it being a 
favorite for laboratory simulation (Bar-Nun 
and Kleinfeld 1989). Reference to condensa- 
tion plots shows that water ice condenses, 
for most solar nebula models, at tempera- 
tures around 140-160°K. well above the 
temperature (roughly 50°K) at which clath- 
rate formation or significant adsorption of 
volatiles would occur. Trapping of gases in 
amorphous ice runs into the difficulty that 
condensation temperatures for water ice in 
the nebula are well above 130°K, which is 
generally recognized as that required to 
make amorphous ice from the vapor (Hobbs 
1974). Note that amorphous ice cannot be 
made by cooling crystalline ice. Nonethe- 
less, some amorphous ice may have sur- 
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vived entry into the nebula and we do not 
dismiss its presence in primordial comets. 

Clathrate formation, while thermodynam- 
icallv plausible, requires good exposure of 
the ice to the gas (as does volumetrically 
significant adsorption). We can think of only 
one plausible process by which fresh ice 
may be exposed to trap gas: relatively slow 
collisions among planetesimals to expose 
fresh ice. Lunine and Stevenson (1985) con- 
sidered this process, and while the numbers 
are marginal, they do not by any means rule 
out this possibility. We adopt as a working 
hypothesis that the volatiles in Comet Hal- 
ley were trapped in water ice in the form of 
clathrate hydrates. 

Lunine (1989) considered the chemical 
fractionation processes which occurred dur- 
ing trapping to alter the abundances of mole- 
cules from' those in the nebular gas. We 
briefly sketch the argument here, consider- 
ing the species methane, carbon dioxide, 
carbon monoxide, nitrogen, and ammonia. 

In the carbon species, methane incorporates 
into clathrate much more readily than does 
carbon monoxide; at low temperatures car- 
bon dioxide thermodynamically prefers to 
form dry ice but some will be trapped in 
the clathrate. We expect a lower CO-to-CH* 
ratio in the clathrate than in the nebula. How 
much lower depends upon the availability 
of cage sites for trapping, and hence on the 
amount of ice exposed to the nebular gas. 
In the case of a small number of available 
cage sites, and hence extreme competition, 
the CO/CH, value is depleted in the clath- 
rate by a factor of IQ 6 or more. The CH 4 
and CO abundances relative to water ice in 
Halley given in Table I constrain for us the 
amount of ice exposed to the gas: it must be 
nearly the maximum possible. Assuming a 
gaseous nebula in which most of the carbon 
is in the form of CO. with CH 4 secondary, 
essentially all of the nebular methane is 
trapped, while only a fraction of the carbon 
monoxide can be accommodated. Under 
such conditions, the CO/CH 4 value is de- 
pleted in the clathrate by a more modest 
factor, 10 to 100. 


Bar-Nun and Kleinfeld (1989) find a si- :- 
lar CO depletion factor for their experiments 
involving codeposition of amorphous ice 
with CO-CH, gas mixtures. While their ex- 
perimental situation is rather complex, the 
results imply that (aj the amorphous ice .s 
able to trap amounts of CH 4 and CO compa- 
rable to clathrate, and (b) the molecular in- 
teractions between the trapped moiecu.e i 
and water ice are similar in clathrate c c 
amorphous ice, which is not unexpected. 

Since carbon dioxide is able to condense 
out as dry ice for the abundances predicted 
by Fig. 1 and plausible quench tempera- 
tures, we assume here that essentially all c . 
the carbon dioxide in the region of Come: 
Halley's formation was available for incor- 
poration into the nucleus. Theretore. -e 
CO-to-CO ; ratio in Halley is predicted tv e 
10 to 100 times less that in the nerula. ag...r. 
because a significant fraction of the CO dees 
not make it into the clathrate. 

Turning to the nitrogen species, we r.r.c 
a slightly more complex situation. Because 
ammonia attaches to water in the form o. 
hydrosen-bonded ammonia hydrates. i-- 
abundance in the ice is decoupled from *.e 
molecular-.nitrogen which only weakly 
corporates into clathrate hydrate. The only 
constraint on the N ,-to-NHj ratio in the 
lev formation zone is the measured ratio o: 
ammonia to water. If it were 15% in Halley , 
then arguably the NHj-to-N; ratio in .r.e 
nebula would be infinite, i.e., all of the ele- 
mental nitrogen was locked up as ammon.d- 
(This assumes, of course, solar eleme. -■ 
N-to-0 abundance and condensation o: .... 
of the water and ammonia.) The abundance 
of ammonia of roughly 0.1% in Halley con- 
strains the nebular N : -to-NH 3 ratio to c= 
roughly 10 to 100. 

F. Postformation Processing 

We have not attempted to model th- 
chemical evolution of species in the nuc.* - 
after formation. While temperatures of m. 
comets stored in the Oort cloud remain low 
[less than 30K (Stem and Shull 1988)]. slow 
diffusion of species within the ice over time 
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v .r during earlier encounters with the Sun. 
cannot be ruled out. The nucleus of Halley 
below the thermal wave is at approximately 
25’K based on observations of the spin tem- 
perature in water (Mumma et «/. 1989). It is 
possible that the molecular inventory ob- 
served in Comet Halley has been somewhat 
modified by postformation processes, but 
-ot to a great extent given the very low 
. .temal temperature. 


3. MODEL RESULTS 


We use the nebular chemistry model 
described above to predict the ratios of 
key molecular species in the region of 
formation of Comet Halley. We assume 
that Halley formed at roughly 20 AU from 
the Sun; recent support for the agglomera- 
jn of at least the short-period comets in 
lids region was presented by Duncan 
(1989); whether this is relevant to Halley 
is problematic. We also assume that the 
chemically active zone in the nebula ex- 
tends out to roughly 1 AU from the Sun. 
We then compute CO/CH*. CO/CO;, and 
N-/NHj at 20 AU as a function of the 
elemental O-to-C ratio in the chemically 
r.ttive zone, assuming that heterogeneous 
catalysis on grains was effective in reduc- 
ing the quench temperature to below 
800°K. Three mixing efficiencies are con- 
sidered. corresponding to full mixing of 
outer and inner nebular gas. 10 Fc mixing, 
and Wo mixing. The outer solar nebula is 
assumed to start with a gaseous inventory 
of CO and any CH 4 , NH 3 . and CO; 
are supplied by inner nebula chemistry. 
V. e revisit this assumption after presenting 
our results. 


Figures 3. 4. and 5 show the results. Be- 
cause we desire to compare nebular gas- 
phase abundances, we do not compare our 
model directly with the Halley data. Rather, 
we apply the fractionation corrections dis- 
cussed in Section 2E to the Halley data. 
Thus, the Halley CO/CH 4 , CO/CO : . and 
N. NHj ratios are all increased by a factor 
of 10 to 100 over the values in Table l. The 
range in the figures is actually larger to ac- 
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Fig. 3. For three different quench temperatures. CO/ 
CH 4 abundance is plotted versus O/C ratio. The thick 
lines are the Halley data corrected for fractionation 
effects, like clathration and adsorption, and uncertain- 
ties in the measurements. The thin lines represent the 
values from our calculations assuming different degrees 
of mixing into the outer solar nebula. It can be seen 
that the amount of CH 4 found in Comet Halley is consis- 
tent with an inner solar nebula origin for quench tem- 
peratures of 600 to 800° K and for a range of nebular 
mixing efficiencies and elemental 0»C values. 
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count for uncertainties in the measured val- 
ues themselves. Results for three different 
quench temperatures are shown. 

Figure 3 indicates that the Halley data are 
consistent with methane manufactured in 
the solar nebula at quench temperatures of 
600 to 800°K for a range of nebular mixing 
efficiencies and elemental O/C values. 
These quench temperatures are all appro- 
priate for catalyzed reactions on grain sur- 
faces; if gas-phase chemistry only were op- 
erative insufficient methane (by orders of 
magnitude) would be available to fit the Hal- 
ley data. Figure 4 paints a similar picture 
for CO : ; although gas-phase reactions can 
occur with quench temperatures as low as 
830°K without catalysis, it appears that 
lower temperatures (and hence reactions on 
grains) are required to reproduce the Halley 
abundance. Note that to fit both the Halley 
CO/CH 4 and CO/CO : ratios for a given in- 
ner nebula elemental abundance and nebula 
mixing efficiencies, somewhat different 
quench temperatures are required in Figs. 3 
and 4. This is not surprising, since the reac- 
tion schemes (and hence rates) are different, 
so that for a fixed nebular lifetime or mixing 
time, different quench temperatures for 
CO-CH 4 and CO~CO : would be expected. 

Figure 5 shows the result for N : /NH 3 for 
the 100% nebular mixing case, which yields 
the most NH 3 in the outer nebula. Clearly, 
insufficient ammonia is produced to account 
for the Halley observations. Quench tem- 
peratures below any reasonable range 
would be required to fit the Halley data. 
It is interesting to note that Marconi and 
Mendis (1988) argue that the Halley results 
are actually consistent with zero ammonia, 
when variations in the solar UV flux (which 
is required for the data reduction) are con- 
sidered. However, the corroborating data of 
Wyckoff et al. (1989) appear to' place the 
ammonia abundance on a firm footing. 

New elemental abundances determined 
by Anders and Grevesse (1989) change the 
O-to-C ratio significantly, from 1.66 to 2.36. 
Running our chemical model for O-to-C ra- 
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Fig. 4. Fracionated Hailey values and data from c r 
own model were calculated in the same way as de- 
scribed in Fig. 3 . A solar origin for CO : in Comet Halley 
requires quench temperatures slightly lower than those 
for CH 4 . 


tios higher than those given in the figures 
results in an increase in the H : 0 and CO^ 
abundances while the CH 4 and CO abur- 
dances decrease. Because CO decreases 
slightly less rapidly than CH 4 , the predicted 
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Fig. J. Distribution of nitrogen between N ; and NHj 
is independent of the O/C ratio. 


CO/CH, value would increase for a given 
:ench temperature. However the CO/CO ; 
number decreases significantly. The reader 
can see these tendencies in Figs. 3 and 4. For 
the lower initial O/C of 1.66 two different 
quench temperatures for C0-CH 4 and 
CO-CO : interconversion were required to 
fit the Halley data. For O/C greater than 
1.66. the quench temperatures for these two 
reaction schemes begin to approach each 
• ■ther. Otherwise, our conclusions are not 
qualitatively affected by the new abundance 
determination. 

In summary, the Halley CO/CH 4 and CO/ 
CO^ ratios are reproduced by a nebular 
chemical model in which heterogeneous ca- 
talysis occurs efficiently and gases are 
trapped in water ice which eventually forms 
Halley’s nucleus. The model NH 3 abun- 
dance is orders of magnitude below that 
jeen in Halley, however. 

4. IMPLICATIONS 

A. Solar Nebula Origin for Volatiles 

in Halley 

The results presented above indicate that 
the abundances of the carbon-bearing spe- 
cies CO. CO : , and CH 4 can be accounted 
for by more or less traditional models of 
solar nebula chemistry involving equilib- 
rium chemistry in the inner disk, quenching, 


and transport of gas parcels outward to the 
region at which Halley formed. The work- 
ability of this model hinges, however, on 
several assumptions or assertions: 

1. Halley formed in the Uranus/Neptune 
region, and not at several hundred AU (this 
can be relaxed if nebular radial mixing is 
assumed to be very efficient: but see below). 

2. Catalysis of chemical reactions on 
grains in the inner nebula was an efficient 
and widespread process. 

3. Gases were incorporated into Halley 
by being trapped in water ice which favored 
incorporation of CH 4 over CO. CO ; was in- 
corporated with high efficiency largely as 
dry ice. 

Note that constraints on nebular mixing 
or the radial distance of the quench zone are 
not very tight: Figs. 3 and 4 show that a 
range of mixing efficiencies could be played 
off against quench zone location. Likewise 
the results do not place firm constraints on 
the inner nebula elemental abundances. 

The inability to explain Halley’s ammonia 
abundance with this model forces us to con- 
sider alternatives. We favor the hypothesis 
that the ammonia in Halley was derived 
largely from the surrounding giant molecu- 
lar cloud. Chemical models and observa- 
tions of giant molecular clouds (Blake et al. 
1987) suggest that NHj may be abundant. 
The amounts of methane and carbon dioxide 
in such clouds are not known. Inability to 
detect methane, even with tight observa- 
tional constraints, implies that either (a) 
CH 4 is locked up in grains or (b) CO vastly 
dominates over CH 4 . 

To argue that Halley incorporated inter- 
stellar ammonia requires that nebular mix- 
ing be less than 100% efficient, otherwise 
essentially all of the ammonia which found 
its way into the nebula would have been 
cycled to high-temperature regions and re- 
cycled outward as molecular nitrogen. The 
recent evidence that formaldehyde is pres- 
ent in Halley (Mumma and Reuter 1989) 
tends to argue in favor of this hypothesis. 
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since that molecule is not expected to form 
in abundance in solar nebula chemical pro- 
cesses. Moreover, there is evidence from 
the Mumma and Reuter data of composi- 
tional heterogeneity across the nucleus of 
Halley. 

We therefore postulate the following situ- 
ation: The nebula is radially mixed only to 
a limited extent. CH a and CO : are manufac- 
tured in the inner part of the disk, trans- 
ported outward past the quench zone to the 
region of comet formation, and trapped in 
ice grains which went to form the nucleus 
of Halley. Along with this gas was incorpo- 
rated some relatively pristine giant molecu- 
lar cloud material, containing NH 3 , which 
was incorporated into Halley either in sepa- 
rate grains or through the gas phase. We 
cannot rule out the possibility that some 
CH 4 was also present in the unprocessed 
molecular cloud gas. Our results indicate to 
us that Halley is a hybrid, containing giant 
molecular cloud gas partially reprocessed in 
the inner solar nebula. 

Prinn and Fegley (1989) have proposed 
that the reduced gases in Comet Halley are 
the result of chemistry occurring in the gas- 
eous subnebulas surrounding Jupiter and 
Saturn, followed by sweep-up of this mate- 
rial (in gaseous or solid form) by solar orbit- 
ing objects. As noted by the authors, the 
dynamics of such processes needs to be 
quantified. We offer the present model for 
Halley's volatiles as an alternative hy- 
pothesis. 

B. Drag Heating or Crains Entering the 

Solar Nebula 

It is of interest to assess whether icy 
grains falling into the solar nebula from the 
surrounding molecular cloud could have 
found their way unaltered to the nucleus 
of Halley. Several processes may modify 
grains during infall: (1) grain-grain colli- 
sions, (2) accretion shock heating at the neb- 
ular boundary, and (3) drag heating of 
grains. We explore the third process here. 

The calculations are based on a paper by 
Wood (1984), in which drag heating was pro- 


posed as a possible means of forming chc.;- 
drules. Here we use the formalism for heat- 
ing of icy grains, and add a term for 
sublimation of water ice. Heat gained by 
aerodynamic drag will be lost by radiation 
and sublimation; the energy balance for the 
grain can be described as 


dt 


4rrr; f o!p,u 3 
me 8 
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where the sublimation rate of ice is given by 


dm _ p 
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T p , r p , v, m, and c are particle temperature, 
radius, velocity, mass, and specific heat: ,u. 
and p, are the gas molecular mass and den- 
sity (the latter set here to a value of 10~ :! 
g/cm 3 ). T 0 is the background radiation tem- 
perature (10°K); L and P s are the latent heat 
of sublimation and the vapor pressure. Fi- 
nally a, s, or, and k are the accommodation 
coefficient (set to 0.5), the thermal emissi' - 
ity, the Stephan-Boltzmann constant, ar.d 
the Boltzmann constant. The equations are 
solved iteratively with the temperature T. 
as the independent variable. After each time 
step the sublimed mass was subtracted from 
the original mass and the new particle radius 
and*ve!ocity were calculated. The radiative 
properties of a panicle depend strongly or. 
the particle size through the emissivitv. Fc 
a 1-mm particle (probably larger than most 
interstellar grains) an emissivitv of 0.75 w as 
assumed (Wood 1984); for r p = 50 nm. e = 
0.01 (Greenberg 1978). Thus small particles 
are poor radiators and hence sublime more 
rapidly than larger particles. 

Particles falling onto the accretion disk 
have infall velocities dependent on their fi- 
nal radial distance from the disk center. A 
final orbital distance of 15-20 AU corres- 
ponds to an infall velocity of 6 km/sec. 
whereas 1 km/sec corresponds to several 
thousand AU. These velocities are conser- 
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Fig. 6. For two different panicle velocities, the heat- 
ing characteristics of these panicles during infail onto 
the accretion disk are shown. 


vative, in the sense that they assume zero 
initial velocity contributed by random mo- 
tions (or energetic processes) within the mo- 
lecular cloud itself. Figure 6 shows the tem- 
perature rise versus time for the two 
different particle sizes given above. For 
each size, curves corresponding to the two 
infall velocities l and 6 km/sec are shown. 
The amount of mass remaining after heating 
and sublimation is shown for each curve, 
expressed as a percentage of the original 
mass of the particle. 

The results indicate that while significant 
release of volatiles accompanies drag heat- 
ing of infalling particles, some mass survives 
in both the large- and small-particle cases. 


However, peak temperatures are high 
enough that molecules loosely bound to the 
interstellar grains would be released into the 
gas phase. Sublimation moderates the tem- 
perature rise to a great degree, such that 
some amorphous ice could have survived 
entry into the solar nebula. Other mecha- 
nisms for heating infalling grains, such as 
accretion shock, have yet to be explored. 

5. CONCLUSIONS AND PROSPECTS 

A model of equilibrium chemistry in the 
solar nebula was combined with considera- 
tions of nebular mixing, elemental abun- 
dance variation, and fractionation effects 
associated with trapping of gases in water 
ice to predict the inventory of major carbon- 
and nitrogen-bearing molecules at 20 AU, 
assumed here to be roughly the region 
within which Hailey formed. For a range of 
nebular mixing efficiencies and elemental 0- 
to-C ratios, the Hailey daca on the CO, CO : , 
and CH 4 abundances could be reproduced 
fairly well. However, the predicted NH 3 
abundance fails far short of that inferred for 
Halley. We adopt the alternative hypothesis 
that the ammonia is preserved from the giant 
molecular cloud within which the solar neb- 
ula is embedded. Infall of grains and conse- 
quent drag heating are sufficient to release 
volatile gases into the nebula, as well as 
sublime much of the water; however, rem- 
nants of interstellar grains, including amor- 
phous ice, could have in part survived. 

The comparison with Halley data is to a 
large extent stymied by uncertainties in the 
true parent molecule abundances in that 
comet and the lack of knowledge of the man- 
ner by which the volatile gases were trapped 
in the nucleus. It is also unclear how typical 
the composition of Halley may be. Perhaps 
not unexpectedly, significant developments 
in ground-based observations of comets ac- 
companied the International Halley Watch 
program and associated spacecraft flybys. 
It is hoped that in the next few years these 
techniques will be brought to bear on a num- 
ber of comets to assemble a good sample 
of inventories of water, ammonia, methane. 
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carbon monoxide, formaldehyde, and other 
volatile species accessible to remote 
studies. 

The next breakthrough in nucleus studies 
may have to await the Comet Rendezvous 
Asteroid Flyby Mission (CRAF). The 
CRAJ 7 penetrator will carry a suite of instru- 
ments designed to determine the ice struc- 
ture and the manner in which volatiles are 
trapped in a subsurface nucleus sample. An 
extensive inventory of parent molecules will 
be accumulated. Other experiments on the 
Mariner Mark II spacecraft will involve de- 
tailed chemical studies of the dust panicles 
emitted from the nucleus. These types of 
observations are required to tightly con- 
strain the type of modeling attempted here. 
It is to be hoped that the calculations pre- 
sented above give at least a flavor of the 
kinds of information which could be derived 
on cometary formation from the CRAF mis- 
sion, and others such as the Rosetta comet 
nucleus sample return, to follow. 
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Vertical Distribution of Water in the Atmosphere of 
Venus: A Simple Thermochemical Explanation 

John S. Lewis and Davi d H. Grinspoon* 

Several lines of evidence concerning the vertical abundance profile of water in the 
atmosphere of Venus lead to strikingly unusual distributions (the water vapor 
abundance decreases sharply in the immediate vicinity of the surface) or to serious 
conflicts in the profiles (different infrared bands suggest water abundances that arc 
discrepant by a factor of 2.5 to 10). These data sets can be reconciled if (i) water 
molecules associate with carbon dioxide and sulfur trioxide to make gaseous carbonic 
acid ;11K I sulfuric acid in the lower atmosphere, ami (ii) the discrepant 0.94-niicmmcter 
water measurements arc due to gaseous sulfuric acid, requiring it to be a somewhat 
stronger absorber than water vapor in this wavelength region. A mean total water 
abundance of 50 ± 20 parts per million and a near-surface free water vapor abundance 
of 10 ± 4 parts per million arc derived. 


T he surface temperature of Ve- 
nus is roughly 750 K, far coo high 
for the presence of liquid water: any 
water present would at once boil away and 
exist as water vapor in die atmosphere. But 
it is also well established diat die atmo- 
sphere of Venus is extremely dry: whereas 
heating Earth to die temperature of Venus 
would produce a very wet atmosphere with 
a surface partial pressure of water vapor of 
about 300 atm, typical estimates of die 
water partial pressure on Venus arc near 
3 x 10' 3 atm, fully 10“ s as much. But 
almost as interesting as die extraordinary 
aridity of Venus is die serious disagreement 
among die results of die several experiments 
diat have attempted to measure the water 
abundance below the clouds. I his problem 
is reviewed in detail by von Zahn et al. (/). 

Very high water concentrations of several 
thousand parts per million were reported as 
a resulf of the Pioneer Venus gas chromato- 
graph (l’VGC) experiment (2), whereas all 
other direct measurements of die water con- 
tent lie in die range from 20 to 200 ppm. 
Tile PVGC experimenters pointed out, 
however, that dicir high value could be 
caused bv die slow evaporation of sulfuric 
add (H 2 S0 4 ) droplets ingested during pas- 
sage dirougii die clouds. Anodicr, less di- 
rect, measurement was carried out by die 
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infrared net flax radiometers on die four 
Pioneer Venus (PV IR) entry probes (J). 
The absolute flaxes measured below die 
clouds were much too large to be reasonable 
and increased widi depth. Large corrections 
were applied to die data on the basis of a 
proposed error source, but even the correct- 
ed data show many cases in which die flux 
increases with depth. The results for the four 
probes were found to be extremely discor- 
dant, suggesting water-mixing ratios from 
10 to over 1000 ppm. For diese reasons, we 
shall attempt no quantitative explanation of 
die PVGC or PV IR radiometer data. 

The most sensitive measurements bearing 
on die water abundance in die lower atmo- 
sphere of Venus, infrared spectra returned 
by die Soviet Venera 1 1 entry probe (4), 
cover a spectral region (0.7 to 1.25 g.m) 
within which one weak (0.82 g.m) and two 
strong (0.94 and 1.13 \Ltn) water bands arc 
found. The intensity of transmitted sunlight 
was measured from an altitude of over 61 
km down to die planetary surface. Below an 
altitude of about 44 km, just below die base 
of die H 2 S0 4 cloud layer, we would normal- 
ly expect die mole fraction of water vapor to 
be independent of altitude. Preliminary in- 
terpretation of these bands (4), however, 
showed clearly diat such homogeneous 
models did only a poor job of explaining the 
data. Models in which the water mole frac- 
tion increased with altitude, from about 20 
ppm at the surface to about 200 ppm at the 
cloud base, gave an improved but nor con- 
vincing fir. The water vapor profiles that 
fitted the 0.94-g.m data lnrst had the unfor- 
tunate property that they seemed quire un- 


physical. In a search for “hydrogen-hiding" 
reactions with other known or suspected 
atmospheric constituents (C, N, O, S, F, Cl, 

As, Sb, and P compounds), no explanation 
for such a strange water distribution could 
be found. 

A much more thorough interpretation of 
the Venera 1 1 data (5) refined die dilemma. 

It was found diat (i) if it was assumed that 
all die opaciry in die Venera 11 spectral 
region was due to carbon dioxide (C0 2 ) and 
H 2 0, die data for die 0.94- jxm band re- 
quired a vertical variation in the water abun- 
dance from abour 20 ppm at the surface to a 
peak value of about 200 ppm near the cloud 
base at 44 km, widi die abundance dropping 
off again above the clouds to below 30 ppm 
near 60 km (Fig. 1 ); (ii) the 0.84- and 1.13- 
p.m bands, although not analyzablc widi the 
same precision, were compatible with aver- 
age water abundances of about 20 ppm, an 
amount insufficient to provide even halt die 
equivalent vviddi of die 0.94-jim band. This 
result clearly implied a contribution from an 
absorber other than water in the 0.94-jcm 
band, but no such species could be identi- 
fied: however, the correlation of the gaseous 
absorption maximum with the cloud base 
suggested the possibility of a gaseous ab- 
sorber formed by evaporation of H 2 S0 4 
cloud droplets (5). If, radicr unreasonably, 
we ignore the estimate of 20 ppm given by 
the 0.84- and 1.13-p.m bands and attribute 
ail of die 0.94-ji.m opacity to water vapor, 
then the average water abundance (integrat- 
ed over the altitude range from 0 to 70 km) 
would be 43 ppm. A lower limit on the 
water abundance would be set by attributing 
die excess (discordant) opacity in the 0.94- 
jLm band near an altitude of 44 km to some 
other species besides water. Then the alti- 
tude- averaged water abundance required to 
explain the rest of the 0.94-g.m opacity is 
near 24 ppm. 

We have developed a model of die chemi- 
cal behavior of water on Venus (6), based on 
laboratory data for the thermodynamics of 
the gns-pha.se reaction 

H 2 0 + SO) = H 2 S0 4 (1) 

on the vapor pressure of H 2 S0 4 (7, 5), and 
on die stability of die gaseous carbonic acid 
(H 2 C0 3 ) molecule, 

H 2 0 + C0 2 * H 2 CO) (2) 

which has recently been characterized (9). 
We use a standard enthalpy change of 
* -22.7 keni mol" 1 and an entropy 
change of AS 0 = —34.42 cal mot 1 K. 1 for 
reaction 1 (6), and AH 4 ^* s 6 kcal mol 
and AS 0 = 2 cal moP 1 K' 1 for reaction 2 
(W). 

Venus atmospheric models were calculat- 
ed for water-mixing ratios ranging from 20 
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Fig. 1. Vertical distribution 
of water on Venus inferred 
from the Venera il and 12 
infrared scanning photom- 
eter for a ten -layer mode) 
atmosphere (5). The Venera 
data (4, 5) show a wide vari- 
ation from layer to layer, 
with a striking minimum at 
the surface and a strong 
“spike’’ near an altitude of 
45 km. Calculated profiles 
of the gases H ? CO,, FLO, 
and hLSO* arc shown (as- 
suming a rotal of 90 ppm 
warer) for the model dc- 
scri I >ct I in the rear. The 
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smoothed cloud density is given for reference on the left side of the diagram (261). 


model predicts that there should be an alti- 
tude interval (19 to 35 km) within which at 
lease 70% of die total water is in the form of 
free H 2 0 vapor. It is satisfying to find that 
die data from these six experiments arc 
concordant widiin (and only within) this 
aitirude interval. We therefore use the four 
reported detections of water in this region, 
plus the lower limits discussed above, to 
estimate a total water content in the range 
from 30 to 70 ppm. 

At die surface die nominal percentage of 
total water bound in gaseous H2CO3 is 
80%. Allowing an uncertainty of a factor of 
3 in the equilibrium constant for reaction 2 
[in accord with the estimates of gaseous 


ro 200 ppm. In our lust models we made 
die simplifying assumption char in die vicin- 
ity of die clouds the abundances of both 
warer and sulfur trioxide (SO\) were derived 
from the evaporation ami dissociation of 
H2SO4 and hence were equal (an excess of 
water is more probable). The results of such 
a calculation for a warer abundance of 90 
ppm are shown in Fig. 1. 

Several generalizations emerge from diis 
and other calculated models. First, the de- 
pletion of water vapor near the ground can 
be ascribed ro die formadon of gaseous 
HjCOj. Second, SO3 and water vapor have 
negligible affinity near the surface, and thus 
H2SO4 vapor reaches I ppm only near an 
altitude of 21 km. Third, gaseous H2SO4 is 
extremely stable in the vicinity of the cloud 
base, reaching 63 ppm (versus 27 ppm for 
water vapor) near 44 km. Hie anomalous 
(excess) 0.94- p.m absorption peak corre- 
sponds almost perfectly with die calculated 
profile of gaseous H 2 S0 4 . This finding 
strongly suggests dur H : S0 4 vapor, or 
some odier gas derived directly from it, is 
responsible for die very high spike in die 
0.94-jxm absorption. Fourth, the column- 
average free water vapor abundance for tills 
model is 36 ppm, compared to the as- 
sumed total water abundance of 90 ppm: 
60% of the rotal water is bound to C0 2 or 
SO3 (principally die former). The percent- 


thc falloll of the warer abundance above 
die clouds is governed by the vapor pres- 
sure of water over concentrated H 2 S0 4 , as 
has been long accepted (/), and no new 
modeling of this familiar phenomenon was 
attempted. 

The altitude at which H 2 S0 4 saturates to 
make constant-boiling (98% by weight) 
I US() 4 droplets is a sensitive function of die 
gaseous H 2 S0 4 abundance. Steffes (8) has 
shown from vapor pressure arguments that 
15 to 30 ppm H : S0 4 would not only 
provide saturation above altitudes of 46 to 
48 km, but that this amount of H 2 S0 4 
vapor would also account for the microwave 
absorption seen in spacecraft occultation 
experiments. 

In light of the present model, it is interest- 
ing to reexamine the other sources of data 
on the water abundance below the Venus 
clouds. Figure 2 summarizes die results 
from the Venera 11 and 12 mass spectro- 
meters (12), the Venera 12 gas chromato- 
graph (7J), die Pioneer Venus Large Probe 
Mass Spectrometer (LMS) (74), and, for 
comparison, the Venera 11 and 12 infrared 
photometer data for the same altitudes (4, 
5). Our model suggests diat the striking 
discordance between these experimental re- 
sults is due to the complex and unexpected 
airitude-dependenr spociarion of water. Our 


H2CO3 abundance given in (9)J, the frac- 
tion of total warer bound as HjCOj at the 
surface of Venus must lie in the range 
80220%- The correcr water partial pressure 
to use in equilibrium treatments of die 
chemical interaction of atmospheric gases 
widi crustal minerals is then lO^ ppm. 
Discussions with Donahue (15) have re- 
vealed char the Pioneer Venus LMS experi- 
ment did not provide dirccr evidence for the 
H2CO3* ion near die surface of Venus, but 
die speciacion of gaseous H 2 CO] in the 
LMS inlet system and ionization region is 
quite unknown. The abundance of free hy- 
drogen seen by the LMS varies rapidly and 
(apparently) inexplicably with altitude in the 
same altitude range where we predict abun- 
dant H2CO3 vapor. These data remain un- 
published. 

We conclude chat die 0.94- p.m band mea- 
sures opacity due not only to C0 2 and water 
vapor but also to some odier cloud-related 
species such as gaseous H 2 S0 4 . The general 
trends evident in the Venera 11 infrared 
data, including both die complex altitude 
dependence of opacity and die impossibility 
of reconciling the 0.84- and 1.13-jxm data 
with free water abundances much larger 
than 20 ppm, arc all readily comprehensible 
when formation of gaseous H 2 S0 4 and 
H2CO3 is considered. Laboratory study of 
the near- infrared absorption spectrum of 


age of free warer is model- independent be- 
cause the equilibrium expressions for reac- 
tion 2 gives a FljOrF^C.Oj ratio that is 
independent of the total water abundance. 
Therefore, die lower limit on free water set 
above by omitting the 45-km "spike'" then 
suggests a minimum total warer abundance 
of 24/0.6 = 40 ppm, and the free water 
abundance of 20 ppm deduced from the 
0.82- and 1.13-jjlui bands implies* a total 
water abundance of 20/0.6 = 33 ppm. Pre- 
cisely at die surface, folly 80% of the total 
water is bound to CO;. All these results arc 
extremely insensitive to plausible errors in 
the thermodynamic data used (7/). Finally, 


Fig. 2. Comparison of water abun- 
dance determinations for the atmo- 
sphere of Venus beneath the main 
cloud layer. The Venera 11 (VII) 
and Venera 12 ( V 1 2^ entry probe 
experiments that provided these 
data arc the infrared photometer 
(4>) {4, 5), the mass spectrometer 
(MS) (72), and die gas chromato- 
graph (GQ ( 15) experiments. The 
I* ii nicer Venus LMS data are given 
with the original authors’ error bars 
(14). Because of the unanticipated 
complexity of the speciacion of wa- 
ter, these experiments yielded high- 
ly tli.witrd.inr results. However, the 
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region IVom |9 ro 35 km, where free water va|M>r dominates and complexes arc minor, yielded 
concordant estimates of the water abundance. 
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HiSOj vnj'Hir and related compounds could 
permir a quantitative interpretation of the 
Venus spectrum, and is srronglv encour- 
aged. 'Hie ncar-surtace association of CO> 
and water vapor is a direct result of our 
acceptance of the data of Tcrlouw et al. (9). 
It is critically important that the equilibrium 
behavior of this system be confirmed by 
further experiments. One can also directly 
test this model by using Venus water pro- 
files to Ik measured by the Near Infrared 
Mapping Spectrometer experiment on the 
Galileo spacecraft (16). Further, die condi- 
tions of temperature and pressure and the 
high CO 2 and water vapor abundances con- 
ducive to extensive gaseous H>C0 3 forma- 
tion are closely similar to those hypothesized 
for early Karth. An improved undemanding 
of the rhermochcmical and optical proper- 
ties of diis species may permit substantial 
improvement of these models. 
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Negative Pressure-Temperature Slopes for Reactions 
Forming MgSi0 3 Perovskite from Calorimetry 

E. Ito, M. Akaogi, L. Topor, A. Navrotsky 


A new and sensitive differential drop solution calorimetric technique was developed 
for very small samples. A single experiment using one 5.18-milligram sample of 
perovskite, synthesized at 25 gigapascals and 1873 Kelvin, gave 110.1 ± 4.1 kilojoules 
per mole for the enthalpy of the ilmenitc- perovskite transition in MgSiOj. The 
thermodynamics of die reaction of MgSiOj (ilmenite) to MgStOj (perovskite) and of 
Mg 2 Si0 4 (spine!) to MgSiOj (perovskite) and MgO (periclasc) were assessed. Despite 
uncertainties in heat capacity and molar volume at high pressure and temperature, 
both reactions clearly have negative pressure -temperature slopes, -0.005 ± 0.002 and 
-0.004 ± 0.002 gigapascals per Kelvin, respectively. The latter may be insufficiently 
negative to preclude whole -mandc convecdon. 


M gSiCVricii perovskite is 
diouglu to dominate die earth s 
lower mantle (/). Knowing its 
stability is crucial for understanding the 
seismic discontinuity near 650 km and for 
predicting lower mantle properties. The 
pressure-temperature (P-T) slope (dP/dT) of 
reactions forming perovskite may help de- 
termine whether die mantle is composition- 
ally layered or whether whole mantle con- 
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vecrion can occur; a strongly negative slope 
would inhibit convection across die bound- 
ary (2, J). The nature of convection has 
important implications for die evolution of 
die earth, its heat budget, and die coupling 
of convection in the core to processes in the 
mantle and crust. 

Ito and Yamada (4) first demonstrated 
thac for die transformations of MgSiOj 
ilmenitc ro perovskite and of Mg : Si0 4 spi- 
nel to perovskite and periclasc, have nega- 
tive slopes (<lP/dT < 0) and that both reac- 
tions occur at similar pressures and tempera- 
tures. fro and Takahashi (5) refined these 
observations. They concluded that the pro- 
nounced sharpness of the 650-km disconti- 


nuity can be reasonably interpreted as the 
dissociation of silicate spinel into perovskite 
plus magnesiowustice. Although their data 
suggest that JP/dT is negative, pressure and 
temperature determinations in multianvil 
devices still have considerable uncertainties, 
and phase boundary reversals arc difficult. 
Thus, the slope of the perovskitc-forming 
reaction has remained somewhat uncertain. 
Because of die importance of this slope in 
considerations of whole mantle convection, 
it is desirable to determine dP/dT by a 
different and independent method, in order 
to test whether the phase synthesis runs 
were indeed close to equilibrium, to confirm 
that the slope is negative, and to obtain a 
reliable value for dP/dT. 

High- tempera cure oxide melt solution 
calorimetry (6) offers a means of obtaining 
thermochcmical data for high-pressure sili- 
cates (7-9). However, in order to apply this 
technique to the very small amounts of 
MgSi0 3 perovskite obtainable by synthesis 
at 25 GPa in a mukianvii apparatus, its 
sensitivity had to be improved. In this re- 
port, wc describe a new differential drop- 
solution technique and report a calorimetric 
derermmarion, using 5 mg of perovskite, of 
the enthalpy difference between MgSiO* 
pyroxene and perovskite. Combination with 
earlier thcmuxlynamic data (7 -10) allows 
calculation of die stability fields of MgSi0 3 
perovskite and comparison with high-pres- 
sure studies. 


14 SEPTEMBER 1990 


REPORTS 1275 


If Ai«;sH9, 93-112(1991) 


Rheological Properties of Ammonia-Water Liquids and 
Crystal-Liquid Slurries: Planetological Applications 

J. S. Kakci-l. S. K. Croft, J. 1. Lunini:. and J. S. Lewis 

Lunar ami Planetary Laboratory, University of Arizona. Tucson. Arizona X572I 
Received April 25. 1990: revised August 24. 1990 


The viscosities of aqueous mixtures plausibly representing a 
range of cryovolcanic substances seen on the icy satellites have 
been measured in the laboratory or obtained from literature 
sources. These viscosities range from 10~ 2 poise (pure water} to 
10 2 * poises (ammonia-water peritectic) to about 10 5 poises (ammo- 
nia-water-methanol peritectic). The viscosities of the liquid mix- 
tures examined in this work are much greater than would be 
expected based on the assumption that the end member molecules 
are noninteractive. This observation supports others based on mo- 
lar volumes and vapor pressure relations indicating that strong 
molecular interactive forces exist and have an important bearing 
on the physical properties of the mixtures. With supercooling and/ 
or partial crystallization, these substances may attain viscosities 
several orders of magnitude greater than those given above. The 
rheological effects of partial crystallization parallel the same effects 
in silicate lavas, so it is reasonable to interpret cryovolcanic mor- 
phologies on the icy satellites in the same ways that we interpret 
remotely observed silicate volcanic morphologies on the Earth and 
terrestrial planets, after accounting for differences in surface gravi- 
ties and lava densities, and allowing for uncertainties in surface 
slopes and extrusion rates. Given the wide range of viscosities for 
simple aqueous mixtures, and the rheological effects of realistic 
thermal states, the characteristics of observed cryovolcanic hows 
and resurfaced plains on the icy satellites can be understood within 
the framework of conventional magmatic processes working on 
exotic icy substances. < 1991 Academic Prcvs. Inc. 

1. INTRODUCTION 

The reconnaissance of the Voyager spacecraft through 
the outer solar system revealed an unexpected diversity 
and intensity of internally driven modifications to the sur- 
faces of icy satellites. Pre-Voyager theoretical arguments 
suggested that some of the icy satellites might have been 
volcanically (“cryovolcanically”) active at some point, 
where the “lava” was water or some cryogenic liquid 
such as ammonia-water (Lewis 1972, Consolmugno and 

Lewis 1978). Voyager’s observations have confirmed 

these early suspicions, but have turther indicated that 

cryovolcanic processes arc complex and varied and have 


combined with tectonic and other processes to produce a 
bewildering variety of bizarre landscapes, each unique in 
its own expression. 

Our present understanding of icy satellites does not lend 
itself to a good predictive understanding of the geological 
styles of icy satellites. The genesis of the many distinctive 
landforms and landscapes must, however, relate in some 
way to important intersatellite variable parameters includ- 
ing satellite radius, surface gravity, rock : ice ratio, heat 
source, and ice composition. Ice and cryolava composi- 
tions and composition-dependent rheologies are likely to 
be of particular importance for satellites with histories 
of cryovolcanism. However, our present knowledge of 
solid-liquid phase equilibria and other physicochemical 
properties of volatile assemblages relevant to cryovolcan- 
ism lags at least 60 years behind our equivalent knowledge 
of silicate properties relevant to terrestrial igneous pro- 
cesses. 

This report is the second from our group attempting to 
redress some important deficiencies. In an earlier paper 
we presented the equations of state for ammonia-water 
liquids (Croft et al. 1988). Here, we present rheological 
data for ammonia-water liquids and partially crystallized 
slurries and apply these data to matters pertaining to cryo- 
volcanism. We extend our report to include ammonia-wa- 
ter-methanol liquids, taking methanol to be representa- 
tive of many possible ammonia-water-soluble chemical 
components likely to be present in minor or trace amounts 
in icy satellites. (Other likely components include 
(NH 4 ) : S, NH 4 CN, <CH 2 )*NV NaCL NH 4 C1, (NH 4 ) : S0 4 , 
amino acids, and many others; Kargel 1988. 1990.) We 
also summarize data from the literature for other sub- 
stances so as to provide an overview of the extreme rheo- 
logical range of plausible cryovolcanic agents. 

Previous work on the viscosities of superliquidus am- 
monia-water liquids include two reports from the Russian 
literature (Pleskov and Igamberdyev 1939, Pinevich 1948) 
and, more recently, an abstract by McKinnon and Mead- 
ows (1984). The Russian workers were concerned with 
viscosities at temperatures above 200 K, although McKin- 
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non and Meadows extrapolated their data down in temper- 
ature and through composition to the ammonia-water 
eutectic and estimated viscosities similar to those re- 
ported here from direct measurements. 

2. EXPERIMENTAL 

The rheologies of ammonia-water and ammonia-wa- 
ter-methanol liquids and liquid-crystal slurries were in- 
vestigated using the falling sphere and rotating disk 
methods. 

2.1 Rotational Vis come try 

These measurements were obtained with an HB Model 
Brookfield viscometer. This instrument measures the re- 
sistance encountered by cylindrical and disk-like spindles 
rotating in a viscous medium contained in a coaxial cylin- 
drical vessel. All of the measurements reported here were 
made with Brookfield’s largest spindle (“HBr*). This 
spindle is classified as a disk spindle, although it actually 
is a short cylinder, 5.6 cm in diameter, 2.3 cm high, and 
mounted on a coaxial cylindrical shaft 0.3 cm in diameter. 
The geometry is rather complex, so the conversion from 
rotation rates and applied torque to absolute viscosities is 
by the manufacturer’s calibration rather than from first 
principles. 

Viscosity was measured as functions of temperature, 
composition, strain rate, and degree of partial crystalliza- 
tion (Table I). Measurements at three or more different 
spindle rotation rates for each composition and tempera- 
ture allowed departures from Newtonian viscous behavior 
to be detected and quantified (Newtonian behavior was 
observed in all samples except for some partially crystal- 
lized slurries, as described below). For a consistent and 
simple presentation, all viscosities reported in Table I 
were either obtained with the spindle rotating at 20 revolu- 
tions per minute, or were reduced to 20 rpm using a power- 
law extrapolation or interpolation of data obtained at other 
rotation rates; to complete the rheological description of 
non-Newtonian slurries the power-law exponent, /?, is 
also given in Table I (these details are irrelevant when the 
substance rheology is Newtonian, in which case n — 1 and 
the viscosity alone is sufficient to describe the rheology). 

The range of reliable viscosities obtained with this in- 
strument is from about 1 to I0 3 poises (10 poise = I Pa- 
sec). Instrumental precision may, in principle, approach 
1% at full scale torque. In practice, measurements were 
repeatable to about ±3% to ± 15%, 8% being most typical. 
These larger figures include the fact that torque was most 
typically not at full scale, and include uncertainties in the 
temperature and composition of the liquid (see below) as 
well as other experimental sources of error. 

Each run used about 1400 ml of solution contained in a 
closed and insulated cylindrical polyethylene vessel with 


TABLE I 

Ammonia— Water Viscosities Obtained by Rotational Viscometry 


Composition 
(wt.% NH 3 ) 

Temperature 

(K) 

% 

Crystallization 

Viscosity 

(poise) 

Power law 
exponent, n 

12.4 

254.05 

5 

2.0 

1.6 

12.4 

252.95 

9 

2.4 

1.4 

12.4 

249.95 

16 

8.0 

1.8 

12.4 

248.95 

18 

7.6 

1.9 

12.4 

245.55 

26 

10.1 

1.5 

14.0 

251.15 

? 

1.0 

1.9 

14.0 

249.75 

6 

1.4 

2.1 

14.0 

244.95 

18 

3.0 

1.8 

14.0 

244.15 

19 

4.0 

2.1 

18.1 

240.05 

5 

9.2 

1.7 

18.1 

238.95 

7 

12.0 

2.1 

18.1 

233.35 

12 

22.0 

2.2 

18.1 

231.45 

17 

44.0 

1.9 

18.1 

227.55 

21 

63 

1.8 

18.1 

221.15 

25 

77 

2.2 

21.6 

226.95 

6 

2.6 

2.1 

21.6 

224.55 

8 

3.6 

2.1 

21.6 

220,45 

13 

4.6 

2.2 

21.6 

213.45 

18 

! 5.4 

2.4 

21.6 

199.55 

25 

51 

1.9 

25.2 

214.75 

3 

2.6 

1.6 

25.2 

211.65 

5 

6.2 

1.8 

25.2 

206.15 

9 

12.4 

1.6 

25.2 

199.35 

13 

17,0 

1.3 

26.7 

209.65 

0 

1.2 

1.0 

26.7 

204.35 

5 

6.6 

1.4 

26.7 

195.65 

10 

15.0 

1.2 

28.4 

207.85 

0 

1.2 



28.4 

202.85 

0 

2.0 

1.0 

28.4 

199.85 

I 

2.8 

1.0 

28.4 

196.65 

3 

5.3 

1.1 

28.4 

191.45 

6 

13.6 

1.2 

28.4 

184.65 

10 

40.0 

1.2 

28.4 

183.75 

10 

40.0 

1.2 

29.4 

209.15 

0 

1.0 



29.4 

205.75 

0 

1.4 

— 

29,4 

205.45 

0 

1.3 

— 

29.4 

202.85 

0 

1.8 

— 

29.4 

202.55 

0 

1.8 



29.4 

198.95 

0 

3.0 

— 

29.4 

198.05 

0 

3.2 

1.0 

29.4 

197.55 

0 

3.3 

1.0 

29.4 

197.35 

0 

4.0 



29.4 

196.35 

0 

4.2 

1.0 

29.4 

195.75 

0 

4.8 

1.0 

29.4 

193.05 

2 

7.2 

1.0 

29.4 

192.15 

2 

8.2 

1.0 

29.4 

190.35 

4 

10.6 

1.0 

29.4 

188.65 

5 

16.0 

1.0 

29.4 

188.05 

5 

16.5 

1.0 

29.4 

T84.45 

6 

40.0 

1.0 

29.4 

183.55 

7 

54 

1.0 

29.4 

182.05 

8 

68 

1.0 

29.4 

177.65 

9 

102 

1.0 

29.4 

176.25 

10 

116 

1.0 

29.4 

175.75 

10 

166 

1.0 

29.6 

195.75 

0 

4.2 

1. 1 

29.6 

185.75 

5 

25.3 

1.0 
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TABLE I — Continued 


Composition 
(wt.% NH ? ) 

Temperature 

(Kj 

% 

Crystallization 

Viscosity 

(poise) 

Power law 
exponent, n 

29.6 

IH 1.05 

8 

61 

1.1 

29.6 

176.15 

— 

488 

1.0 

33.0 

202.25 

0 

1.7 

— 

33.0 

190.15 

0 

6.6 

1.0 

33.0 

185.15 

0 

13.9 

1.0 

33.0 

180.75 

0 

22.4 

l.l 

33.0 

175.85 

— 

97 

1.0 

33.0 

176.25 

— 

78 

1.0 

33.0 

176.15 

— 

93 

1.0 

35.5 

204.55 

0 

1.3 

— 

35.5 

200.25 

0 

1.9 

— 

35.5 

195.85 

0 

3.2 

— 

35.5 

190.25 

0 

6.2 

1.0 

35.5 

183.95 

0 

14.0 

1.0 

35.5 

183.15 

0 

15.0 

1.0 

35.5 

182.25 

0 

17.4 

1.0 

35.5 

180.45 

0 

24.8 

1.0 

35.5 

179.15 

0 

29.4 

l.l 

35.5 

177.95 

0 

32.6 

1.0 

35.5 

177.65 

0 

39.4 

— 

35.5 

176.45 

0 

35.6 

1.0 

35.5 

175.55 

0 

43.2 

1.0 

35.5 

173.85 

0 

50 

l.l 

35.5 

176.45 

0 

40.0 

— 

35.5 

173.45 

0 

70 

1.0 

35.5 

171.55 

0 

74 

1.0 

42.0 

203.35 

0 

1.2 

1.0 

42.0 

198.45 

0 

1.9 

l.l 

42.0 

194.25 

0 

2.9 

1.0 

42.0 

189.45 

0 

5.3 

l.l 

42.0 

188.75 

0 

6.4 

1.0 

49.1 

199.35 

0 

1.2 

1.0 

49.1 

196.75 

0 

1.6 

1.0 

49.1 

194.85 

0 

1.8 

1.0 

57.0 

190.65 

0 

U 

— 

57.0 

187.95 

0 

1.6 

1.0 


Note. All viscosities were obtained with the HB l spindle operating at 
or reduced to 20 rpm. A dashed entry under the crystallization column 
indicates that partial crystallization occurred at an invariant point pre- 
venting use of the lever rule to calculate % crystallization. A dashed 
entry in the power law exponent column indicates that the data needed 
to calculate n were not obtained. 

a 12.2-cm internal diameter. This vessel was slowly cooled 
in liquid nitrogen, and was opened for short periods of 
time to stir and thermally homogenize the liquid, to make 
viscosity and temperature measurements, and to observe 
crystallization. The total duration of each run was typi- 
cally 15 hr, with a mean cooling time on the order of 2- 
4 K per hour. 

Temperature was measured with chromel-alumel ther- 
mocouples and platinum resistance thermometers; these 
instruments were calibrated at the freezing points of wa- 
ter, methanol, and acetone, the sublimation point of car- 
bon dioxide, and the boiling points of nitrogen and water 


(the latter three temperatures corrected to lake into ac- 
count the laboratory air pressure). The accuracy of the 
temperature measurements was about ±0.2 K. However, 
it was difficult to thermally homogenize the more viscous, 
low-temperature mixtures: this increased the effective 
temperature error to as much as ±0.5 K for measurements 
below about 1 85 K, and ±Q.S K for measurements below 
170 K. 

2,2 , Falling Sphere Method 

The classical analytical approximation of Stokes gives 
the viscosity, /x, of a viscous medium if the settling rate 
of a solid sphere is measured 

p * 2 A p«V9t/, (I) 

where A p is the density of the solid sphere minus that of 
the liquid, a is the radius of the sphere. # is the gravita- 
tional acceleration, and U is the terminal velocity of the 
sphere. Equation (I) is a very good approximation if the 
flow around the sphere is a very slow creeping motion 
where the Reynolds Number. Re, is much less than I . Re 
is given by 

Re * ZtipfUlfi. (2) 

where p, is the density of the fluid. Various theoretical 
treatments have offered analytical corrections to the 
Stokes equation for non-creeping flow at Re of order I and 
higher, although none of these solutions closely reflects 
reality. Therefore, it is the usual practice to apply empiri- 
cal corrections. The present data have been reduced ac- 
cording to the empirical correction factors of Le Clair et 
al. (1970). Because the correction factors are Re-depen- 
dent and Re is p-dependent. and (i is not initially known, 
it was necessary to apply an iterative routine starting 
with the Stokes approximation to converge on the true, 
corrected values of Re, the correction factor, and the true, 
corrected viscosity. 

The Stokes Eq. (1) and empirical Re-correction factors 
apply only in an infinite fluid medium. Rigid boundaries, 
in this case those of the container, retard the motion of 
the sinking sphere, so wall corrections must also be ap- 
plied. The wall effect was accounted for based on the 
work of Fidleris and Whitmore (1961). 

To extend viscosity measurements over many orders 
of magnitude spheres of various sizes and densities were 
utilized. These included Cr-steel, glass, nylon, and ultra- 
high molecular weight polyethylene spheres ranging from 
0.6 to 1 .3 cm in diameter. Plastics were utilized in the runs 
at high temperatures in order to minimize the density 
contrast of the sphere relative to the liquid. Steel and glass 
were selected to maximize the density contrast for runs 
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at low temperatures where the liquid became extremely 
viscous. The temperature-dependent densities of these 
materials were determined to within ±0.0004 g cm' 3 . 
Small air bubbles sometimes were entrained in the liquid 
or adhered to the spheres; in some instances this problem 
may have contributed significant error to measured vis- 
cosities. 

It was important to thermally equilibrate the sphere 
with the liquid before measuring its descent rate, both to 
achieve the correct density contrast with the liquid, and 
to avoid having the sphere disturb the temperature and 
viscosity of the sample adjacent to the sphere. 

Large hemispherically headed steel bolts, normally 
used on railroad ties and in wooden bridges, were used in 
place of spheres to obtain two semi-quantitative viscosity 
estimates of order !0 6 poises. The effective density of the 
bolt heads was calculated as the mass of the bolt divided 
by the volume of a sphere with equivalent cross-sectional 
area. 

The uncertainty in the liquid temperature caused by 
thermal heterogeneity effectively limited the precision of 
viscosities measured in excess of about 100 poises. The 
uncertainties in the density contrasts between the plastic 
spheres and the viscous media, and the combination of 
large and somewhat uncertain empirical corrections for 
the wall effect and finite Re limited the precision of viscos- 
ities below about 1 poise. Hence, the falling sphere 
method yielded its best results between 1 and 100 poises, 
the same viscosity-temperature-composition region 
where most rotational viscometric data were obtained. In 
general, the quality of the falling sphere data is not as 
good as the quality of the rotational viscometry except 
near I poise. However, the falling sphere method yielded 
useful results over a broader viscosity range, from 1CT : 
to 10 6 poises. 

23. Liquid Preparation 

The compositional range for the ammonia-water liquids 
investigated in this study is from about 12.6 to 56.0 wt% 
NH 3 . Ammonia-water solutions with compositions rang- 
ing up to 29.9 wt.% NH 3 were prepared by diluting analyti- 
cal reagent grade ammonium hydroxide (Maffinckrodt, 3 
ppm total impurities) with distilled water (maximum 1 
ppm impurities). Ammonium hydroxide was precooled to 
about 260 K to reduce ammonia vapor loss during transfer 
from the bottle. The composition of the ammonia-water 
reagent was determined by two independent methods: 
measurement of the liquid's specific gravity (Croft et al. 
1988) and measurement of its liquidus temperature. 

Solutions with NH 3 contents greater than 29.9% were 
prepared by mixing ammonium hydroxide solution with 
anhydrous liquid ammonia (Mutheson, maximum 100 ppm 
total impurities, mostly water); since mixing was done at 


ambient pressure it was necessary to mix these solutions 
at low temperatures to prevent boiling. Ammonia was 
cooled by the Joule-Thompson effect to the normal boil- 
ing point (240 K) upon extrusion from the high-pressure 
cylinder. Ammonia and ammonia-water reagent were 
then cooled to the freezing point of pure ammonia (196 
K) prior to mixing. The quantities of the liquids used in 
preparing the solutions were measured volumetrically, 
and their masses were determined using the tempera- 
ture-composition-density relationship of Croft et al. 
(1988). Determinations of the liquidus temperature before 
and after experimental runs indicated that ammonia vapor 
loss and atmospheric water vapor condensation over the 
course of each run amounted to a reduction of the am- 
monia concentration from the initial value by about 
0. 1-0.3%, absolute. The uncertainty in the initial compo- 
sitions is ±0.2% NH 3 . Therefore, one may take ±0.4% 
as the uncertainty in the liquid composition for a given 
viscosity measurement. Ammonia-water-methanol liq- 
uids were prepared in the same general fashion, using 
analytical reagent methanol (Mallinckrodt, maximum 
0.07% impurities). 

3. RESULTS 

3.1 . Ammonia-Water Liquids above the Liquidus 

Figure 1 and Table I show temperature-dependent vis- 
cosities of ammonia-water solutions obtained by rota- 
tional viscosity. Figure 2 and Table II show the viscosities 
obtained by the falling sphere method. Figure 2 also shows 
a curve, discussed below, fit to all superliquidus data 
obtained in this study (including rotational viscometry) 
and a fit to the data of Pinevich (1948), both curves for 
30% NH 3 solution. 

Temperature-composition-viscosity fit. Least- 
squares numerical fitting procedures (Bevington 1969) 
were applied to the viscosity measurements to obtain in- 
terpolation equations. The basic form of the fitting equa- 
tion adopted is: 

p. = exp{A + BIT }, (3) 

where fi is the viscosity. T is temperature, and A and B 
are material constants. The form of Eq. (3) is the one 
commonly adopted to express viscosity data for both liq- 
uids and gases (Touloukian et al. 1975), and has both 
theoretical and empirical support. Effectively, exp{A} is 
a reference viscosity at very high temperatures and B 
is an "activation temperature" related to an activation 
energy. Both A and B vary with ammonia content of 
ammonia-water liquids, thus in our fit, A and B are repre- 
sented by polynomial expressions in the mass fraction, X , 
of ammonia in the mixtures. Both earlier work (McKinnon 
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FIG. 1. Effective viscosities of ammonia-water liquids and slurries 
obtained by rotational viscometry using Brookfield's largest spindle 
operated at 20 rpm. 


and Meadows 1984, Kargel 1987) and our initial attempts 
to fit both the high-temperature data (Pinevich 1948, 
Pleskov and Igamberdyev 1939) and our data (Tables I 
and II) indicated that B was a function of temperature 
also. Attempts to fit the values of B as a smooth function 
of temperature were unsuccessful. Measurements brack- 
eting the apparent transition near 240 K indicated a sharp, 
rather than a gradual transition (Fig. 2). Therefore, two 
separate fits were made: one to data above 240 K (includ- 
ing the mixed NH^-H : 0 data of Pinevich 1948, and 
Pleskov and Igamberdyev 1939; pure H : 0 data of Haar et 
al. 1984; and pure NH 3 data of Alei and Litchman 1972), 
and a second to data below 240 K (Tables I and II and 
supercooled H : 0 data from Weast and Selby 1967). The 
least-squares numerical fits are 

A = - 10.8143 + 0.71 1062 X - 22.4943 AT : 

+ 41.8343 X } - 18.5149 3f 4 (4a) 


B = 1819.86 + 250.822 X + 6505.25 X 2 

- 14923.4 X' + 7141.46 A' 4 (4b) 

for T > 240 K, and 

A = -13.8628 - 68.7617 X + 230.083 X 2 

- 249.897 X- (5a) 

B = 2701.73 + 14973.3 A' - 46174.5 A'- 

+ 45967.6 X- (5b) 

for T < 240 K. RMS deviations from Eq. (4) for the data 
used in the fit are about 5%. RMS deviations from Eq. (5) 
for the viscometer data (Table I) used in the fit is about 
10%. Deviations of the Stokes Flow data (Table 11) from 
Eq. (5) is about 15%. 

Isoviscs computed from Eqs. (3), (4), and (5) are shown 
in Fig. 3. The computed transition curve between the two 
solutions is also shown. As with density and bulk modulus 
(Croft et al. 1988), viscosities of ammonia-water liquids 
show an anomalous maximum near the eutectic composi- 
tion (see below). Approximate mean activation tempera- 
tures are about 2050 K for T > 240 K. and about 4200 K 
for T < 240 K. The corresponding activation energies are, 
respectively, about 17 kJ/mole and 35 kJ/mole, which 
agree reasonably well with the results of McKinnon and 
Meadows (1984). Further discussion of this transition is 
given in Section 4. 1 . 

3.2 . SuperCooIed Ammonia-Water Liquids 

The ability to chill liquids below their normal freezing 
points without crystallization is a common metastable 
phenomenon (supercooling). The experimentalist must 
use some measure of patience to prevent significant super- 
cooling of liquids close in composition to the ammo- 
nia-water peritectic. In fact. Rupert (1909). in the earliest 
determination of the complete solid-liquid phase diagram 
of the ammonia-water system, extrapolated supercooled 
“liquidus” data from higher temperatures leading to the 
erroneous conclusion that a eutectic occurred in this sys- 
tem near 158 K, 17 K below the actual value. Several 
subsequent determinations also erred on the low side. 

The correct eutectic composition and temperature (35.4 
wt.% NH 3 , 175.4 K) and the existence of ammonia dihy- 
drate and of the ammonia-water peritectic (32.6% NH 3 , 
176.2 K) were first established by Rollet and Vuillard 
( 1956). The early errors, due to strong metastability in the 
ammonia— water system, have caused a fair amount of 
confusion in the planetary science community as to the 
actual minimum melting temperature, with published tem- 
peratures of 170, 173, 175, and 176 K. For primitive, 
previously undifferentiated ammonia-water ices, the ac- 
tual temperature of melting at low pressures is the peritec- 
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FIG. 2. Viscosities of ammonia-water liquids obtained by the falling sphere method. Curves represent all data obtained by both methods in 
this study and data of Pinevich (1948). 


tic, 176.2 K, rather than the eutectic as is commonly 
stated. The supercooling phenomenon continues to frus- 
trate many experimentalists when crystallization is de- 
sired. However, under certain circumstances supercool- 
ing may occur in natural cryovolcanic settings and may 
be an important factorjoveming the rheologies of these 
cryolavas. 

Figures I and 2 show the viscosities of supercooled 
ammonia-water liquids near the peritectic in composi- 
tion. It is possible to maintain ammonia-water liquids 
several K below the normal freezing point for several 
hours as long as they are not too vigorously stirred and 
as long as ice nucleii are not introduced by atmospheric 
condensation or from moist thermocouple probes or vis- 
cometer spindles. It was, therefore, very easy to get good 
viscosity measurements by either viscometric technique 
for slightly supercooled liquids. 

The more highly supercooled mixtures become highly 


viscous and thermally inhomogeneous, requiring ener- 
getic manual mixing, long thermal equilibration times, 
or small sample sizes in order to obtain more thermally 
homogeneous mixtures. Invariably, vigorous manual mix- 
ing triggers crystal nucleation within a few minutes, often 
before the liquid can be thermally homog enized. Some - 
times crystallization may be triggered merely by inserting 
a thermocouple probe or dropping a steel sphere into the 
liquid. The great difficulty encountered in trying to stir 
and otherwise manipulate these mixtures severely limited 
the number and quality of viscosity determinations for 
highly supercooled mixtures and made them impossible 
to study by rotational viscometry. The measured viscosi- 
ties for the strongly supercooled mixtures (Fig. 2 and 
Table II) were all obtained by the descending sphere 
method using liquids lacking any observable signs of crys- 
tallization. The greatest uncertainty in these measure- 
ments is the temperature, which typically varied by 1 K 
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TABLE II 

Viscosities of Ammonia-Water Mixtures Obtained by the 
Descending Sphere Method 


Composition 
(wl.% NH 3 ) 

Temperature 

(K) 

Stokes 

viscosity 

(poises) 

Corrected 

viscosity 

(poises) 

21.4 

238.6 

0.121 

0.069 

21.4 

239.6 

<0.117 

0.065 

21.4 

240.0 

0.120 

0.067 

21.4 

240.7 

0.122 

0.069 

21.4 

241.5 

0.144 

0.060 

21.4 

241.5 

0.110 

0.057 

21.4 

243.5 

0.105 

0.052 

21.4 

244.5 

0.098 

0.046 

21.4 

245.6 

0.095 

0.044 

21.4 

248.7 

0.092 

0.042 

21.4 

251.6 

0.085 

0.035 

21.4 

251.7 

0.085 

0.035 

21.4 

255.2 

0.081 

0.038 

23.2 

234.1 

0.162 

0.097 

23.2 

235.3 

0.161 

0.095 

23.2 

236.8 

0.149 

0.084 

23.2 

237.3 

0.136 

0.069 

23.2 

238.0 

0.140 

0.073 

23.2 

238.4 

0.132 

0.065 

23.2 

240.5 

0.124 

0.058 

23.2 

244.6 

0.118 

0.050 

23.2 

245.6 

0.102 

0.038 

23.2 

247.0 

0.110 

0.045 

23.2 

247.6 

0.104 

0.040 

23.6 

223.6 

0.210 

0.146 

23.6 

230.0 

0.180 

0.117 

23.6 

234.1 

0.160 

0.097 

23.6 

239.7 

0.130 

0.065 

23.6 

249.6 

0.112 

0.046 

23.6 

253.3 

0.089 

0.027 

24.7 

250.6 

0.124 

0.052 

24.7 

244.6 

0.134 

0.062 

24.7 

238.0 

0.134 

0.063 

24.7 

233.0 

0.195 

0.121 

24.7 

225.9 

0.241 

0.164 

24.7 

223.4 

0.505 

0.233 

24.7 

220.8 

0.636 

0.334 

24.7 

219.8 

0.665 

0.345 

25.0 

217.8 

0.392 

0.305 

25.0 

219.6 

0.387 

0.295 

25.0 

234.1 

0.183 

0.1 11 

25.0 

236.4 

0.147 

0.074 

25.0 

238.8 

0.143 

0.070 

25.0 

248.3 

0.120 

0.048 

25.0 

257.9 

0.101 

0.031 

27.3 

211.1 

1.19 

0.715 

27.3 

216.3 

0.866 

0.477 

27.3 

220.1 

0.716 

0.370 

27.3 

222.9 

0.317 

0.225 

27.3 

226.6 

0.220 

0.138 

27.3 

232.0 

0.199 

0.114 

27.3 

236.2 

0.163 

0.080 

27.3 

238.2 

0.179 

0.094 

27.3 

241.3 

0.138 

0.058 

27.3 

243.4 

0.145 

0.063 

27.3 

246.3 

0.134 

0.052 


TABLE II — C oniimwU 


Composition 

(Wt Sr NHi) 

Temperature 

(K> 

Stokes 

viscosity 

(poises) 

Corrected 

viscosity 

(poises) 

27.3 

248.3 

0. 145 

0.0h2 

29.0 

159.9 


102,000 S 

29.0 

161.5 


29,400 S 

29.0 

162.5 


49,500 S 

29.0 

167.9 


3,570 S 

29.0 

199.9 

3.38 

2.29 

29.0 

203.5 

1.58 

1.31 

29.0 

206.7 

1.35 

1.1 1 

29.0 

209.3 

1.29 

0.784 

29.0 

209.6 

0.811 

0.658 

29.0 

210.6 

1 . 1 1 

0.652 

29.0 

211.9 

1.01 

0.813 

29.0 

212.7 

0,673 

0.538 

29.0 

213.2 

0.660 

0.527 

29.0 

215.8 

0.514 

0.400 

29.0 

216.3 

0.936 

0.530 

29.0 

217.2 

0.409 

0.309 

29.0 

217.5 

0.407 

0.300 

29.0 

217.6 

0.915 

0.497 

29.0 

219.9 

0.449 

0.344 

29.0 

220. 1 

0.420 

0.314 

29.0 

220.4 

0.416 

0.315 

o 

ON 

r i 

221.4 

0.331 

0.243 

29.0 

222.5 

0.285 

0.196 

29.0 


0.304 

0.221 

29.0 

225.7 

0.254 

0.165 

29.0 

226.8 

0.294 

0.201 

29.0 

228.3 

0.229 

0.146 

29.0 

232.0 

0.216 

0.125 

29.0 

235.9 

0.171 

0.083 

29.0 

239.2 

0.164 

0.075 

29.0 

241.5 

0.155 

0.066 

29.3 

147.7 

4.070.000 

2.300.000 S 

29.3 

154.2 

481.000 

270,000 S 

29.3 

189.3 

9.07 

7.89 

29.3 

197.3 

3.72 

3.23 

29.3 

200.1 

2.22 

1.92 

29.3 

203.2 

1.49 

(.28 

29.3 

203.8 

1.70 

1.13 

29.3 

207.7 

0.887 

0.74 

34.6 

163.9 


1,490 S 

34.6 

163.3 


2.060 S 

34.6 

162.9 


2.850 S 

34.6 

162.2 


4,300 S 

34.6 

180.6 


20.4 

34.6 

178.6 

34.5 

24.3 

34.6 

175.6 

68.3 

48.0 

34.6 

164.6 

614 

412 S 

34.6 

163.6 

893 

600 S 

34.6 

175.8 

1025 

688 17 r r crystals 

40.4 

193.9 


5.92 

40.4 

193.4 


5.93 

40.4 

192.9 


5.73 


Note. S< supercooled without crystallization. Corrected viscosities 
have been corrected for the wall effect and finite Reynolds number. The 
magnitude of the corrections can be assessed by comparison with the 
Stokes viscosities which have not been corrected. 
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FIG. 3. Temperature-composmon-log M) viscosity fit for ammo* 
nia-waier liquids. Dashed curve shows the transition temperature where 
the activation energy changes. 


through the sample. Considering the sharp temperature 
dependence of the viscosities in this thermal regime, this 
introduces a considerable uncertainty in the viscosities 
for any given temperature. 

Figure 2 shows that the viscosities of strongly super- 
cooled ammonia-water liquids exhibit far more powerful 
compositional and temperature dependencies than do the 
superliquidus mixtures. This behavior is explained in 
terms of glass formation in Section 4.1. 

3.3 , Partially Crystallized Ammonia-Water Slurries 

Complex rheologies . In contrast to the relatively sim- 
ply Newtonian viscous behavior of the crystal-free liquids 
discussed above, partially crystallized slurries are gener- 
ally characterized by an extremely complex rheology. 
Even 1 or 29c partial crystallization may introduce a rich 
variety of rheological phenomena, including hysteresis 
effects _and, most importantly, strain-rate-dependent 
(shear stress dependent) effective viscosities. The com- 
plex rheological behavior of the slurries makes it difficult 
to exactly reproduce the measured viscosities. For exam- 
ple. the strain rate dependence ensures that the same 
geometry and shear rate of a viscometric measurement 
must be reproduced exactly in order to reproduce the 
same results. 

The thermal history also affects the viscosity of ammo- 
nia-water slurries. The effective viscosities of well-stirred 
slurries maintained at constant temperature generally de- 
crease slowly, some tens of percent, over periods of about 
10 hr. Part of this phenomenon is probably thixotropy 
(shear thinning) and is common for particle suspensions 
in liquids. However, thinning also occurs to some extent 


for mixtures which are left undisturbed at some constant 
temperature below the liquidus, and for mixtures which 
are slightly warmed after cooling without otherwise dis- 
turbing them. Such "time-thinning” is probably related 
to the slow recrystallization of suspended ice, which 
may initially tend to have dendritic (snowflake-like) forms 
characteristic of quenched crystals; more compact euhe- 
dral and resorbed forms may physically interact less ex- 
tensively than quenched crystals, for a given crystal pack- 
ing fraction, therefore yielding lower effective viscosities. 
To maintain as self-consistent a data set as possible, the 
effective viscosities of slurries reported in Table I were 
measured while cooling the mixtures at more or less 
steady rates. 

Simple models of complex rheologies. The rheologi- 
cal behaviors of crystal-liquid suspensions are commonly 
approximated as "Newtonian” viscous substances, as 
"power-law" ("pseudo-plastic”) substances, or as "Bin- 
gham plastics,” the latter two comprising common types 
of non-Newtonian substances. These are sometimes use- 
ful approximations to the rheological behavior of naturally 
occurring slurries (e.g., mud and lava flows), although 
they often apply only over limited ranges of deformation 
rates and applied shear stresses. 

A Newtonian substance is characterized by a linear 
dependence of strain rate, k. on applied shear stress, cr, 
and by the lack of a yield strength. Its viscosity, /x, is 
defined simply as 

fi = cr/e. (6) 

A Bingham plastic is also characterized by a linear 
dependence of the strain rate on the applied shear stress, 
but unlike the Newtonian substance, it also displays a 
finite yield strength. A Newtonian-like effective viscosity, 
£i cfr , can still be described for a Bingham plastic using 
Eq. (6), although this effective viscosity varies with the 
applied shear stress (unlike Newtonian substances). The 
Bingham viscosity, /x B , is a constant (for constant temper- 
ature) given by 

= (cr - n/i, (7) 

where Y is the yield stress. When Y <€ cr a Bingham plastic 
appears to behave essentially as a Newtonian fluid. When 
cr < K then no permanent plastic deformation occurs and 
the substance appears to be an elastic solid. 

A pseudoplastic substance lacks a yield strength, and 
the strain rate has a power-law dependence on the applied 
shear stress 

k = Ae-W RT ir n H (8) 
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where A is a constant for a given temperature, Q is the 
activation energy for viscous flow, and R is the gas con- 
stant. The exponential term gives the temperature depen- 
dence of the strain rate, and n gives the power-law depen- 
dence of the strain rate on the applied shear stress 
(Newtonian and Bingham substances may have similar 
temperature dependencies). A value ol n = 1 signifies 
Newtonian viscous behavior. As the effective viscosity is 
simply <rlk (which is constant only in the special case 
when n = I), we may rewrite Eq. (8) in terms of the 
effective viscosity, 

^ = c r'- n *V RT lA. (9) 

Equations (7) and (9) indicate that the effective viscosities 
of Bingham and power-law substances are functions ot 
the applied shear stress, and in each case give higher 
effective viscosities as the shear stress and strain rate are 
decreased. On a graph of strain rate vs shear stress a 
Newtonian liquid plots as a straight line passing through 
the origin, a Bingham plastic plots as a straight line with 
some positive intercept on the shear stress axis, and a 
power-law substance plots as a curve passing through the 
origin. In each case the effective viscosity is given by the 
slope of a straight line passing through the origin and the 
point of interest on a diagram of stress vs strain rate. 

Newtonian character of peri tec tic slurries. Fortu- 
nately, the slurries with the greatest planetological inter- 
est (those close in composition to the peritectic) do not 
exhibit the most troublesome of rheological phenomena 
(the strain-rate dependence of effective viscosities), so we 
shall discuss these first. 

It is useful to define a relative viscosity. p r 

Mr = m7Mi- (l0) 

where is the viscosity of the slurry, and /x, is the viscos- 
ity of the intergranular liquid phase in the absence of 
suspended crystals. 

Equation (10) has broad relevance only if the slurry 
is Newtonian. This requirement is satisfied by slurries 
containing 29% or more NH-, (at least for the ranges of 
crystallization and shear rates examined in this work). 
Figure 4 shows the dependence of fi x on the percentage of 
crystals, (f>. in ammonia-water slurries. <l> was estimated 
either by the petrological "lever rule' when there w j us a 
reasonably close approach to equilibrium crystallization 
and when the temperature remained above the peritectic. 
or by the temperature rise associated with the latent heal 
released by partial crystallization of supercooled mix- 
tures. was measured, while jx, was calculated for the 
equilibrium liquidus composition at the temperature ol 
the slurry. Like terrestrial lavas, the rheology of amino- 



FIG. 4. Dependence of relative viscosity on degree of crystallization. 
Open circles are ammonia-water slurries. Silicate lava curves are from 
McBirney and Muruse (1984). Theoretical Einstein-Roscoe (E-R) for- 
mulations are for uniform solid spheres suspended in a viscous liquid 
(C = 1.35). and fora serial sphere size distribution tC = 1), and modified 
to approximately fit the ammonia-water slurry data iC = 3.7). 

nia— water liquids is extremely sensitive to its thermal 
state. Further discussion of this matter is saved for Sec- 
tion 4.2. 

Figure 5a shows a stress-strain plot for a near-peritectic 
mixture at various degrees of partial crystallization. No 
significant departures from ideal Newtonian behavior are 
discernable. The data can be fit with a power-law where 
n = I (the Newtonian condition). A yield strength was 
not evident for peritectic and near-peritectic mixtures. 
However. Murase et al. ( 1985 ) observed an apparent yield 
strength for Ml. St. Helens dacite only after 70% crystalli- 
zation. We have not examined slurries with anywhere 
near this quantity of crystals. 

Power-law behavior of water-rich slurries. As the wa- 
ter content of ammonia-water slurries increases from the 
peritectic their rheologies increasingly depart from New- 
tonian. with the apparent viscosity describable to a first 
approximation with a power-law dependence on the strain 
rate (applied shear stress). Ammonia-water compositions 
considerably more water-rich than the peritectic were dif- 
ficult to work with at high degrees of crystallization owing 
to their strongly non-Newtonian, "slushy " character. At 
lower degrees of partial crystallization these slurries had a 
"creamier'’ consistency and their rheologies were readily 
studied with the rotational method. The lever rule can be 
applied to water-rich samples to precisely determine the 
extent of crystallization as long as the temperature of the 
slurry remains above the peritectic point. 

Figure 5b shows representative stress-strain rate plots 
obtained for a relatively water-rich mixture* at various 
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Spindle Rotation Rate (rpm) 


FIG. 5. Strain rate (proportional to spindle rotation rate) vs shear 
stress (proportional to dial reading), (a) Newtonian behavior of slurries 
containing 29,4-29.6% NH 3 . (b) Power-law pseudo-plastic behavior of 
slurries containing 18.1% NH 3 . X indicates fraction of crystals, n indi- 
cates power-law index (equation 8 and Table I). 


degrees of partial crystallization. These slurries obey a 
simple power-law relation to a reasonable approximation. 
The Bingham approximation completely fails to charac- 
terize ammonia-water slurries, at least for the ranges of 
crystal contents and shear rates investigated here. 

Figures 1 and 4 show that abrupt surges in the effective 
viscosities of water-rich ammonia-water mixtures occur 
as these substances partially crystallize below their liq- 
uid!. Supercooled liquids do not display this phenomenon. 
This surge in effective viscosity corresponds to the onset 
of non-Newtonian (n > 1 ) rheological behavior. Figure 6 
shows the dependence of the exponent n (from Eqs. (8) 
and (9)) on the degree of crystallization for several compo- 
sitions. The exponent n was calculated based on viscosity 
measurements made at several different spindle rotation 
rates. 

Non-Newtonian power-law viscosities characterize the 
more water-rich slurries, with n ranging up to about 2. It 


was surprising to find that n attains its full value with 
only a few percent crystals, and either levels off or even 
decreases slightly as the degree of crystallization in- 
creases (the effective viscosity, however, certainly does 
continue to rise sharply as crystallization ensues). This 
behavior could be a hysteresis effect as described above, 
where the first few percent of quenched crystals gradually 
take on more compact, euhedral forms after several hours 
of slow cooling, thereby minimizing interactions with the 
liquid and with other crystals. 

Although this work does not support a Bingham rheol- 
ogy for ammonia-water slurries, it is possible that such a 
rheology would be found at a lower range of strain rates 
or at higher crystal contents than investigated here. Quali- 
tative observations may support this speculation; a stick- 
slip type of behavior was found for the more water-rich 
slurries where the ice crystal content exceeded about 
25%. That is, the spindle was found to stick and surge in 
an unpredictable type of rotation. Although this made 
rheological measurements unreliable, it may indicate that 
a yield strength was being attained. Alternatively, stick- 
slip behavior may have been indicating that the spindle 
was simply cutting a hole in the highly viscous (pseudo- 
plastic) medium, causing the spindle to surge through free 
space, followed by the collapse of the sticky material 
onto the spindle. Or, if the modified Einstein-Roscoe 
formulation discussed in Section 4.2 is approximately cor- 
rect, then small heterogeneities in the amount of crystals 




FIG. 6. Strain rate dependence of expressed as the power-law 
exponent n . as a function of the degree of crystallization. 
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could cause radical local variations in effective viscosity 
(lumpincss). 

3.4 . Ammonia-Water— Methanol Liquids 

Methanol (CH 3 OH) is one of the simplest and most 
abundant organic molecules in interstellar space (Irvine 
ct cd. 1985) and is readily produced in large quantities by 
ion bombardments and ultraviolet irradiations of carbon- 
bearing water ice (Khare ct cd . 1989, Roessicr and Nebe- 
ling 1988). Methanol and many other molecules may possi- 
bly constitute minor or trace components of icy satellites. 
Since the electric dipole moment of methanol (1.70 D) is 
similar to the dipole moments of water ( 1 .85 D) and ammo- 
nia ( 1 .47 D), we would expect a priori that methanol would 
interact with ammonia-water in extensive liquid solution. 
Hence', if they are present in at least minor quantities 
(even hundredths or tenths of 1%) in icy satellites, metha- 
nol and other strongly polar substances may contribute to 
the volcanological complexity and geomorphic diversity 
of icy satellites (Karge! 1990). This work explores some 
of the rheological implications when just this one of the 
many plausible minor components is added to the ammo- 
nia-water system. 

Kargel (1990) mapped out the liquidus surface of the 
ternary system H : 0-NH 3 -CH 3 0H in composition-tem- 
perature space at I atm. (Fig. 7). The detailed structure 
of this system, particularly near the methanol corner and 
near the thermal minima, are known only in approximate 
form, owing to the difficulty encountered in crystallizing 
and manipulating these highly viscous mixtures at such 
low temperatures. The ternary equivalent of the ammo- 
nia-water peritectic has been bracketed in tempera- 



FIG. 7. Liquidus surface of the system H4)-NH;-CH ? ()H (from 
Kargel 1990). Isotherms in °C. Arrows on coleclies indicate the direction 
of liquid migration upon cooling. 


TABLE III 

Viscosities of Ammonia-Water-Methanol Liquids (by 
Rotational Viscometry) 


Liquid compositions 
(weight fractions) 

Temperature 

(K) 

Viscosity 

(poises) 

A B 

(Equation (3)) 

H : 0 

NHi 

CH,OH 

0.486 

0.207 

0.307 

186.1 

26 

-34.78 

7108.2 




177.4 

266 






167.6 

1820 



0.215 

0.092 

0.693 

187.5 

1.8 

-28.44 

5390.6 




179.7 

4.3 






171.1 

[4.8 






164.1 

75 






158.0 

388 



0.241 

0.000 

0.739 

178.3 

2.2 

-21.29 

3936.9 




168.1 

8.4 



0.379 

0.090 

0.531 

198.1 

2.7 

-31.17 

6328.6 




189.9 

6.8 






184.1 

20.0 






177.9 

103 



0.442 

0.227 

0.331 

210.7 

1.0 

-29.80 

6188.8 




202.9 

2.3 






201.5 

2.9 






201.2 

-> 2 






200.2 

3.5 






198.3 

4.6 






195.8 

5.8 





* 

192.9 

7.8 






192.8 

9.0 






190.7 

12.6 






189.0 

17.1 






185.8 

23.6 






182.0 

54 






177.7 

122 






174.3 

320 






172.2 

370 






170.2 

1040 






165.1 

3260 




ture-composition space by melting experiments (giving 
the peritectic temperature) and by extrapolations of the 
liquidus temperature-composition curves (giving the peri- 
tectic composition, and less reliably, the peritectic tem- 
perature). The peritectic liquid appears to contain about 
33% CH 3 OH, 23% NH 3 , and 44% H : 0 by mass, and 
apparently freezes near 1 53 K (all values have large uncer- 
tainties at this point). 

Table III and Fig. 8 show the viscosities of ammo- 
nia-water-methanol liquid mixtures. Some of these mix- 
tures became extremely bubbly at the lower temperatures, 
and it is possible that the bubbles affected some viscosity 
measurements. Pure methanol and ammonia end mem- 
bers are less viscous than liquid water for a given tempera- 
ture. The ternary mixtures, however, are considerably 
more viscous than simple noninteractive mixing would 
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FIG. 8. Viscosities of ammonia-water-methanol liquid (near ter- 
nary peritectic? as a function of temperature, obtained by rotational 
viscometry. 


suggest. Hence, powerful molecular interactions among 
the three components are indicated by the viscosity 
data. 

This inference is supported by the molar volumes of 
ternary mixtures, which indicate highly nonideal solution 
behavior and closer intermoleculur spacings than in any 
one of the pure liquids or any of the binary liquid mixtures. 
For example, the density of a liquid close in comparison 
to the ternary peritectic (H : 0 : NH 3 : CH 3 OH = 
0.470:0.201 : 0.329 by mass) is 0.8706 gem' 1 at 298 K. 
This is 6% denser than is calculated based on the densities 
of the three pure end members and the assumption of 
noninteractive mixing. It is also 29c denser than calculated 
from the densities of pure methanol added to ammo- 
nia-water liquid assuming no interaction between the 
methanol and the ammonia-water. Croft et aL (1988) had 
earlier shown that ammonia and water interact strongly 
in solution. We infer from these density data that methanol 
interacts strongly with water and ammonia molecules, 
although less strongly (gram for gram) than the ammo- 
nia-water interaction. 

The ammonia-water-methanol peritectic mixture has 
an extraordinarily high extrapolated viscosity near the 
freezing point, about 4 x I0 4 poises, three orders of mag- 


nitude greater than the viscosity of ammonia-water peri- 
tectic liquid (which is itself 3 to 4 orders more viscous 
than either pure water or pure ammonia at their freezing 
points). This viscosity estimate is a factor of several un- 
certain because of the uncertain temperature of the peri- 
tectic and because of uncertainties in the extrapolation of 
the data (a nonlinear fit in Fig. 8 would result in a greater 
viscosity at the peritectic and reduced rms deviation). The 
extremely high viscosity of the ternary peritectic com- 
pared to the ammonia-water binary peritectic is primarily 
due to the lurge additional freezing-point depression im- 
posed by methanol. Ternary molecular interactive forces 
appear to be secondarily responsible. Just a few percent 
partial crystallization or a few degrees of supercooling 
make this mixture impossible to manipulate. 

3.5 . Comparisons of Planetologically Important Liquids 

and Ices 

Figure 9 shows the viscosities of liquid water, ammo- 
nia-water, ammonia-water-methanol, water-magne- 
sium sulfate (Kargel, 1990), and methane-nitrogen (Tou- 
loukian et aL 1975) mixtures at their respective eutectic 
and peritectic freezing points. Figure 9 also compares 
these viscosities to those of common silicate lavas at their 
liquidus temperatures. The silicate lavas are certainly not 
plausible volcanic agents on the icy satellites, but are 
included in Fig. 9 to provide a geologically familiar frame 
of reference. All of the viscosities shown in Fig. 9 are for 
crystal-free liquids. Volcanologically realistic amounts of 
supercooling and partial crystallization can increase these 
viscosities by many orders of magnitude, as shown in the 
sections above. Therefore, the viscosities in Fig. 9 are 
best utilized simply to show r the approximate relative dif- 
ferences in the viscosities of the different types of lavas 
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FIG. 9. Viscosities of cryogenic and silicate liquid mixtures. All data 

are for crystal-free mixtures. Silicates are at their respective liquidi: 
crvogenic liquids are at their respective eutectic and peritectic points. 
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FIG. 10. Lab samples of aqueous liquids and slurries of possible 
cryovolcanic importance, (al H : 0-NH? liquid near the perileetic compo- 
sition, T = 176 K, 6 r /c ice crystals, * 100 poises, lb) Supercooled 
H 2 O-NH 1 liquid also containing a few percent ice crystals. T = 160 K. 
»*I0 4 poises, (c) H;O-NHi-CH;0H peri tec tic liquid. T ~ 163 K (about 
10 K above peritectic temperature). **10* poises. 

and cryogenic substances. Assuming that superheating of 
cryolavas above the liquidus is minor or absent (this is 
true of silicate lavas), then the liquidus viscosities in Fig. 
9 correspond to lower limits on the actual eruptive viscos- 
ities. 

Figures lOa-lOc qualitatively illustrate the viscosities 
of some cryovolcanologicully plausible ammonia-water 
liquids in various thermal states. The ammonia-wa- 
ter-methanol liquid (Fig. 10c) could not be manipulated 
near its freezing point due to its highly viscous nature, so 
the experimenter had to settle for a photograph of the 
sample taken about 10 K above its peritectic temperature. 


As seen from Fig. 9 ammonia-water and ammonia-wa- 
ter-methanol peritectic liquids have about the same vis- 
cosities as basaltic and daeilic lavas, respectively, and are 
many orders of magnitude more viscous than pure water, 
salt brines, or methane-nitrogen mixtures at their respec- 
tive freezing points. 


4. THEORETICAL INTERPRETATIONS AND 
PLANETOLOGICAL APPLICATIONS 

4,L Theoretical Interpretation of Ammonia-Water 

Liquid Rheology 

Viscosity vs composition . Liquid water is a highly 
structured phase which, at temperatures well above the 
freezing point, exhibits some long-range order consistent 
with the crystal structure of hexagonal ice I (Eisenberg 
and Kauzmann 1969). Likewise, ammonia, also hydrogen- 
bonded, exhibits structure in the liquid. It is of interest to 
understand whether the mixed liquid is ordered, and what 
the nature of that ordering is. Viscosity data can contrib- 
ute in this area. In the present section, we interpret the 
viscosity data as indicating very strongly that the liquid is 
a mixture of water molecules and ammonia-hydrate and 
dihydrate molecules, rather than simply free water and 
ammonia molecules. We also show weak evidence for 
further structure in the liquid; the following section of 
viscosity versus temperature provides a more compelling 
argument for additional structure in the liquid. 

In an ideal solution, the viscosities of the two end mem- 
bers characterize the viscosity of the mixture according 
to the equation (Hummel and Arndt 1985) 

log /x mix = .v log p. Y + (I - .r) log /Z-, (11) 

where is the viscosity of component y, and .t is its 
mole fraction. In turn, the viscosity of each individual 
component at temperature T is given by Eq. (3) which can 
also be written as 



where S is the configurational entropy of component ,v. 
and A. B ^ are constants (compare with Eq. 3). In the case 
that the solution is not ideal, an additional term is added 
to Eq. (II) which expresses the contribution of the en- 
tropy of mixing of the solution. The notion that the viscos- 
ity may be expressed solely in terms of entropy corres- 
ponds to the view that material transport in liquids 
depends on the number of configurations attainable by the 
system, and that the energy required to shift from one 
configuration to another is not the rate-limiting step of the 
transport process. Such a view works well for certain 
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silicate melts such as plagioclase, even down to the glass 
temperature (Hummel and Arndt 1985). We adopt this 
view for our preliminary analysis of the viscosity-compo- 
sition data and find that it provides a useful perspective. 

We assume first that the ammonia-water liquid, above 
any of the eutectic, peritectic, or congruent melting 
points, is an ideal mixture of ammonia and water mole- 
cules. We then use the fits to the viscosity data given in 
the previous section to compare the viscosity as a function 
of composition to that predicted by the ideal mixing equa- 
tion, with the viscosities of pure water and pure ammonia 
defining the end members. We use the fits to the data 
rather than the data themselves because the trends are 
easier to see using the fits; the goodness-of-fit is thor- 
oughly described above and the reader should refer to this 
section again before drawing any conclusions. 

Figure 11a shows the result of this analysis at 280 K. 
The solid line is the viscosity; the dashed line is that 
predicted from ideal mixing. The dashed curve is up to 
65% lower in magnitude than the solid line. Discrepancies 
above the 10 to 20% level are significant, as is the disagree- 
ment in the sign of the slope for ammonia mole fractions 
less that 25%. The direction of the discrepancy is toward 
a negative entropy of mixing; i.e. , the real liquid has fewer 
configurations available to it than the ideal mixing would 
predict; hence, the real viscosity is higher than the ideal. 

The simplest approach to introducing a restriction in 
the number of available configurations is to assume that 
the liquid is composed not of water and ammonia mole- 
cules in any configuration, but instead that the molecules 
are restricted to hydrogen-bonded clusters according to 
composition. These clusters plausibly correspond to the 
stoichiometric, solid ammonia hydrate compounds. Thus 
between mole fractions of 0 and 33% the end member 
liquid components are H : 0 and 2H : O.NH 3 ; between 33 
and 50% the end members are 2H : O.NH 3 and H : O.NH 3 . 

Figures 1 lb and 1 Ic show the results of using the ideal 
mixing model with the endmember components indicated 
above. In some sense, the improved agreement is almost 
a trivial one in that we are taking smaller compositional 
slices of the viscosity curve in constructing the dashed 
line, and hence a better fit is to be expected. However, 
the great improvement in the fit, and the consistency in 
the sign of the slope, suggest that the result should be 
regarded as significant. The dashed curve fits the solid 
curve to within 6% in both the water-dihydrate, and dihy- 
drate-monohydrate regions (Figs. 1 1 b and 11c, respec- 
tively). This is much better than the scatter in the data 
and hence, to within the accuracy of the viscosity data, 
these dashed, "ideal” curves represent the viscosity data 
from 0 to 50% ammonia at 280 K. One is tempted to argue 
that the systematic, small negative residual entropy of 
mixing suggests additional structure in the liquid, but 
more accurate data are required to test this assertion. 



ommonio mole fraction 

FIG. 11. Viscosity {in poises) versus ammonia mole fraction at 280 
K, shown by the solid lines based on fits described in the text. The 
dashed line in (a) is a model viscosity based on ideal solution of individual 
ammonia and water molecules. The dashed line in (b) and (c) is a model 
viscosity based on ideal solution, in which the endmembers represent 
the compositions of the congruent melting points (or equivalently, the 
solid hydrates). 


We have performed the same test at 220 K, the lowest 
temperature at which sufficient data are available to make 
accurate comparisons, and reach the same conclusion. 
The end member water-ammonia model provides a poor 
fit to the viscosity data (discrepancies of up to 90%); using 
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the congruent compositions as end members reduces the 
errors to 17% or less at 220 K. Note that in this ease, we 
have had to extrapolate the pure water viscosity to well- 
below the freezing point which undoubtedly introduces 
additional errors at this low temperature. 

The conclusion we draw from Fig. 1 1 is that the viscos- 
ity of the ammonia-water system, well into the liquid 
regime, can be interpreted in terms of a liquid structure 
composed of clusters of molecules corresponding to the 
composition of the solid hydrates (or equivalently the 
congruently melting end members). Thus from 0 to 33% 
ammonia, free ammonia molecules do not exist; instead 
the smallest ammonia unit in the liquid is 2H : O.NH 3 . 
From 33 to 50%, the smallest units in the liquid are dihy- 
drate and monohydrate “molecules/’ Evidence for any 
additional structure in the liquid which limits accessible 
configurations is slim or nonexistent. To search for such 
additional order, higher precision in the data at low tem- 
peratures is required. However, additional evidence for 
order is provided in the following subsection, which treats 
the dependence of viscosity on temperature. 

Viscosity vs temperature . We fitted the viscosity data 
(for 33% ammonia liquid) with the expression: 

p. a Qxp(BIT) (13) 

While in the previous section we argued that the viscosity 
is given as a function of composition by a composite 
expression involving at least two terms. we are only inter- 
esting here in fixed composition and variable temperature. 
In effect, we are assigning B to be an “activation energy” 
(normalized by the gas constant) in the sense of a simple 
Newtonian fluid, and asking what is its variation with 
temperature. No more significance than this should be 
attached to the fit. 

The data show two interesting characteristics when fit- 
ted to Eq. (13). First, there is a significant change in 
activation energy at around 240 K. wherein below 240 K 
B = 4200, and above 240 K B = 2050 as described in 
Section 3.1. In both temperature regions B is nearly con- 
stant with temperature. Below 176 K. the subcooled liquid 
has a value of B which increases steeply with decreasing 
temperature below the periiectic. presaging the onset of 
glass formation at some lower temperature. The change 
in activation energy near 240 K suggests a structural 
change in the liquid involving a loss of attainable configu- 
rations at temperatures below the transition. Unfortu- 
nately. with the present data, we cannot constrain well 
the width (in temperature) of the transition: however, wc 
are confident that the data show a discrete rather than 
gradual change in the liquid structure. The effect is remi- 
niscent of the changes in physical properties which water 
undergoes as it approaches its freezing point, in which 


significant long-range order becomes rather established 
with decreasing temperature (Fisenberg and Kauzmann 
1969). 

It therefore appears that, in addition to the formation 
of single-unit hydrate clusters in the liquid over a large 
range of temperature, some long-range order involving 
many clusters becomes established at temperatures ap- 
proaching the congruent melting point of the hydrates. 
Further investigation of this effect requires higher preci- 
sion in the viscosity data. 

Likewise, the onset of the glass transition is suggested 
by the peritcctic-composition data below 176 K. In this 
regime the viscosity can be described by an equation 
which takes account of incipient glass formation at a tem- 
perature T s 

" 4I 

where D is a constant. Preliminary calculations were per- 
formed to estimate a glass temperature based on the vis- 
cosity data and are in fair agreement with the observed 
value, measured by Van Kasteren, of 120 K. However, 
the scatter in the viscosity data permits acceptable fits 
producing glass temperatures up to 140 K. Therefore, 
more precise data are required before we can assert with 
certainty that we are seeing the onset of glass formation 
in the peritectic liquid. 

42 . Volcanological Implications of Shiny Rheology 

Figures I and 4 illustrate an important rheological con- 
sequence of partial crystallization. The sharp increase in 
viscosities starting with the onset of crystallization is very 
similar to that observed for silicate lavas, particularly 
basalt. However, it is likely that relative viscosities of 
lavas in general depend not just on the crystal fraction but 
probably also on the size distribution and form of the 
crystals. Obviously, skeletal, dendritic, and rodlike crys- 
tals would interact differently (more extensively) than 
spherical or other equant crystals would. Delicate skeletal 
and rod-like crystallites form when liquid water is 
quenched ( Hobbs 1974) and when water vapor condenses 
rapidly below the frost point (forming snow and frost). 
Slow crystallization or recrystallization after rapid precip- 
itation in lavas and glacial ice causes larger, more equant 
crystal forms to develop so as to minimize the surface 
area and associated free energy. We expect that a similar 
behavior would occur in ammonia-water and other aque- 
ous slurries and that this would have rheological effects. 
Consistent with this idea, wc have observed that the vis- 
cosity of ammonia- water slurries is somewhat dependent 
on its thermal history (Section 3.3). 

Having expressed our concern that the viscosity of a 
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particular composition probably depends on much more 
than just its temperature and crystal content, we can still 
state with surety that the slurries we examined behaved 
Theologically very much like silicate lavas. This is not a 
trivial observation since it is by no means obvious that 
hydrogen-bonded aqueous mixtures should behave simi- 
larly to silicate lavas. The Fact that these substances do 
exhibit similar rheological behavior, at least under the 
conditions we examined, implies that it is reasonable to 
interpret cryovolcanic morphologies on the basis of expe- 
rience gained with terrestrial and extraterrestrial silicate 
lava flows. 

The relative viscosities of ammonia-water slurries and 
of silicate lavas are much greater than predicted based 
on Einstein-Roscoe approximations assuming spherical 
particles (Graham et al. 1964). A mere 25% crystallization 
can increase the viscosities of ammonia-water and silicate 
slurries by three orders of magnitude relative to the homo- 
geneous liquid phase. This increase can be modeled (see 
“modified E-R" curve in Fig. 4) by a modified form of 
the Einstein-Roscoe formulation 

Mr = Mid - (15) 

where the coefficient C is assigned a value of 3.7 (this 
compares to C = 1.35 for uniform spheres traveling in 
aggregates). C physically represents the increased effec- 
tiveness of highly aspherical crystals and crystal aggre- 
gates in blocking much of the liquid volume from partici- 
pating in the deformation of the slurry. 

According to Eq. (15), C = 3.7 means that a slurry 
would behave as a solid when <£ attains a value near 
30% crystals (compared to about 70% crystals for uniform 
spheres where C = 1,35, and 100% crystals for a serial 
distribution of sphere sizes where C = I). Possibly this is 
the critical crystal fraction when a basalt or ammonia-wa- 
ter slurry attains a yield strength, or at least appears to 
have one. Below this critical crystal fraction, where all of 
our measurements were made, we observed a Newtonian 
rheology for near-peritectic slurries, and a power-law rhe- 
ology for water-rich slurries (Section 3.3). 

Even in the total absence of a yield strength the effective 
cessation of motion of a cooling flow could be quite 
abrupt, occurring when the chilled flow front and margins 
crystallize the critical volume fraction of solids. When 
this happens the viscosity of the chilled crust of the flow 
increases sharply, even as the less crystallized interior 
slurry maintains a relatively low viscosity. This effect 
could be approximated with a Bingham rheology involving 
a yield strength. However, on the basis of our experiments 
and the data of others (e.g., Shaw 1969, and Murase et al. 
1985), we suspect that the commonly utilized Bingham 
model is simply a convenient way to approximate ther- 
mally and Theologically heterogenous flows having power- 


law rheologies. The final dimensions oflava flows depend 
on the dynamical problem of the characteristics of the 
eruptive event, the viscous run-out of the flow, and the 
cooling behavior of the flow margins. Any true yield 
strength probably does not appear until the final moments 
when the lava flow has already virtually attained its final 
shape. This could explain the widely observed correlation 
between lava viscosity (and composition) and apparent 
yield strength (McBirney and Murase 1984. Murase et aL 
1985, Fink and Zimbelman 1986, Moore and Ackerman 
1989). Nevertheless, the yield strength concept is widely 
used among volcanologists, and we feel that it does pro- 
vide useful results. 

4.3. Cryolava Rheology and Comparative Planetology 

Surface gravity is probably one of the more important 
variable parameters affecting the evolution and appear- 
ance of planets and satellites. Of present concern are the 
effects that intersatellite variations in surface gravity may 
impose on the dimensions and morphologies of cryovol- 
canic flows. Variations in lava viscosity affect lava flow 
morphology opposite to variations in the product oflava 
density and surface gravity. This is evident, for example, 
in problems dealing with Newtonian flow down an inclined 
plane or channel flow (Turcotte and Schubert 1982). Since 
simple rheological models of lava flows often seem to 
give incorrect results (evidently because the dynamical 
rheological structures oflava flows are not simple!), we 
find it easiest and safest to perform a first-order gravity 
scaling of cryogenic liquids so as to arrive at a gravity- 
scaled rheological analog among terrestrial lavas. Based 
on simple flow models we define a “mobility index“ = 
log m (ptf/M)* Liquids with the same mobility index should, 
to a reasonable approximation, behave in Theologically 
similar ways. 

Figure 12 shows the mobility indices of silicate and 
cryogenic liquids plotted against surface gravity. Obvi- 
ously our choice of plotting parameters ensures that there 
will be a positive correlation with a slope of unity. The 
purpose of this figure is simply to provide a visual guide 
for rheological gravity scaling. 

For example, from Fig. 12 we see that liquid water 
and salt brines on Ganymede would be about 10 times 
“runnier" than the most fluid lavas (komatiite) on Earth, 
thus probably producing submeter flow thicknesses on 
Ganymede. Mixtures of nondipolar molecules such as 
methane and nitrogen on Triton would be even more fluid. 
Ammonia-water on Enceladus should be roughly compa- 
rable to andesite on Earth, perhaps producing flows on 
the order of a few tens of meters thick on low slopes. 
Ammonia-water-methanol on Miranda or Ariel should 
behave like rhyolite on Earth, producing flows up to sev- 
eral hundred meters thick on low slopes. 
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FIG. 12. Gravity-scaling of viscosities of volcanic and cryovolcanic 
liquids. In each case the viscosity at the freezing point (liquidus) and 
in the absence of suspended crystals is represented. Data for volatile 
mixtures are the eutectic and peritectic compositions. Silicate lava vis- 
cosities are from the first author's compilation of experimental data for 
natural and synthetic rock compositions, and the temperatures for which 
these viscosities apply are the eutectic temperatures in synthetic water- 
free multicomponent mixtures approximating natural rock compositions 
(Morse 1980). 


Another important aspect illustrated by Fig. 12 is that 
while intersatellite variations in surface gravity may be an 
important variable factor controlling surface morpholo- 
gies, cryolava compositions are far more important if 
these compositions vary significantly through the outer 
Solar System. 

4.4. Observed Cryovolcanic Morphologies 

The observed morphologies of volcanic deposits on the 
terrestrial planets have been used to place crude con- 
straints on the magma rheologies (e.g., Pieri and Baloga 
1986, Dragoni et uL 1986, Fink and Bridges 1990, and 
references therein). The techniques are admittedly imper- 
fect, but seem to give reasonable results. Not all of the 
techniques can be applied to flows on the icy satellites 
because not all of the morphologic features found on ter- 
restrial flows are seen on the ice flows, probably because 
of the generally poor resolution (order I- 10 km) of avail- 
able images of icy satellite cryovolcanic deposits. A few 


techniques are useful, however. Fink (1980) presented a 
theory from which the internal viscosity of a volcanic flow 
can be estimated from the amplitude and wavelength of 
drag folds (corda or festoons) on the flow surface, when 
present. Hulme (1974) presented a method of estimating 
the yield strength of a flow from the height and curvature 
of the flow’s convex or lobate edge. Since effective yield 
strengths seem to correlate with viscosities, the latter 
gives a qualitative estimate of viscosities as well. Assum- 
ing that these techniques can be applied to deposits on 
the icy satellites, the observed morphologies imply a 
broad range of effective viscosities at the time of emplace- 
ment for deposits on various satellites. 

The coronal deposits on the Uranian satellite Miranda 
are interpreted as cryovolcanic flows (Croft 1987, 1988) 
with lobate edges hundreds of meters thick and surface 
ripples with amplitudes of hundreds of meters and wave- 
lengths of a few kilometers. Assuming the ripples are 
corda, internal viscosities of 10*-10 9 poise are obtained. 
The lobate deposit dimensions imply yield stresses of 
0.01-0. 1 bar. Both measurements are consistent with each 
other and imply rheologies comparable to terrestrial an- 
desites. On Ariel, there are long linear and lobate flows 
with convex edges on the order of a kilometer high, imply- 
ing yield strengths on the order of a bar. No corda are 
visible, but comparison with the inferred yield strength 
indicate probable viscosities in the range of I0 in poise. 
These high viscosities imply that the erupted substances 
on Miranda and Ariel were composed of some low-melt- 
ing-temperature aqueous mixture, probably involving 
several chemical components and perhaps a high degree 
of crystallinity. High as they are. these viscosity estimates 
are orders of magnitude below the 10 16 poises obtained by 
Jankowski and Squyres (1988) for the flows on Ariel. 
Straightforwardly interpreted viscosities of this order 
would suggest that the flows were erupted as solid ice (as 
the title of their paper indicates) rather than as volcanic 
slurries. However, Jankowski and Squyres pointed out 
that their viscosities are upper limits, and considered al- 
ternatives to solid ice including ice mobilized by intersti- 
tial cryogenic liquids. Jankowski and Squyres' models 
assumed a homogeneous and purely viscous rheology, 
and did not adequately consider the critically important 
rheological effects of the chilled crust of the flows on their 
mobility. Following Greenberg et al. (1990) and Schenk 
( 1989), who suggested that a 4-6 order-of-magnitude crust 
correction factor must be applied to the results of simple 
flow models, we consider (on the basis of Jankowski and 
Squyres’ uncorrected results) that eruptive viscosities on 
the order of I0 ,n -I0 i: poises are more likely for Ariel, 
and an order of magnitude less for Miranda, in rough 
agreement with the results based on corda and yield 
strengths. 

Neptune's large satellite Triton has both flat deposits 




no 


KARCjHK kt al. 


with no discernable lohatc edges and thick, lohate flows 
and ridges. The former are apparently low-viscosity de- 
posits, possibly like those on Eneeladus and Tethys, 
w hereas the latter are high viscosity deposits with rheolo- 
gies comparable to the stiff flows on Ariel and Miranda 
(Kargel and Strom 1990). 

Several Saturnian satellites show evidence of resurfac- 
ing. Most spectacular are the broad, smooth plains of 
Eneeladus. These plains arc nearly featureless at the best 
image resolutions (2 km per line-pair), except where they 
arc interrupted by arcuate and linear ridges (Squyres et 
al. 1983, Pozioand Kargel 1989). The lack of obvious flow 
margins and fronts in these plains, coupled with the very 
low surface gravity of Eneeladus, indicates that the cryo- 
volcanic liquids must have possessed much lower viscosi- 
ties than the substances which formed the Iobate deposits 
on the Uranian and Neptunian satellites. However, con- 
sidering the limitations of image resolution, substances as 
viscous as ammonia-water are allowed. 

Another large cryovolcanic plain is found on Tethys to 
the east of Ithaca Chasma. The plain is recognizable in 
contrast to the rugged cratered terrain due to its flat sur- 
face and somewhat lower crater density. The western 
boundary of the plain with the cratered terrain is very 
irregular (Moore and Ahern 1983) and suggests flooding of 
the rugged upland to a constant topographic level. There 
is no visible raised edge to the plains deposit, again sug- 
gesting a relatively low viscosity at the time of em- 
placement. 

Another style of cryovolcamsm is seen on Ganymede, 
the largest of the icy satellites. Its surface consists ot 
moderately cratered dark plains and light, grooved, 
smooth materials. The light materials were early recog- 
nized to be cryovolcanic deposits on the basis of age. 
stratigraphy, and embayment relationships (e.g.. Luc- 
chitta 1980, Header*//. 1981. Parmentier et al. 1982). The 
rougher dark plains were later recognized to consist of 
a complex of darker cryovolcanic flows based again on 
embayment patterns, age relationships, and albedo varia- 
tions (Casucchia and Strom 1984. Croft 1985, 1986). The 
edges of the light deposits are generally defined by tec- 
tonic structures. In the few cases where embayment rela- 
tions are not tectonically disturbed, there is little or no 
evidence for a raised edge where topographic variations 
on the order of 100 m would easily be seen. The bound- 
aries of albedo patches within the dark terrain (which 
presumably represent individual flows or distinct groups 
of flows) range from relatively sharp to quite diffuse. The 
highest resolution images show little or no vertical defini- 
tion to the edges of the albedo patches or to the smooth 
dark patches within the dark terrain (interpreted as local 
flooding by Casacchiaand Strom 1984). again where topo- 
graphic variations of 100 m would be visible. The thinness 
of the flows and the general lack of definable flow edges 
indicate low viscosity flows with maximum yield stresses 


of 0.01-0. 1 bar, and probably much less. The flows may 
have been emplaced as liquid water, slush, or brine with 
viscosities of 10 : poise, although ammonia-water cannot 
he ruled out based on the available data. Observations 
by Galileo will place firmer bounds on the rheology and 
composition. 

The cryovolcanic surfaces seen on the icy satellites 
vary enormously in their morphologies, and accordingly 
wc infer a great range in their rheological properties at the 
times of eruption. These properties are nonetheless within 
the range of viscosities for likely cryomagma composi- 
tions. In particular, the higher inferred viscosities on the 
icy satellites can easily be reached by partially crystallized 
cryogenic binary materials such as ammonia-water or by 
cryogenic ternary mixtures, while very water-rich sub- 
stances, salt brines, or ice-free ammonia-water can ex- 
plain the lower inferred viscosities. A direct correspon- 
dence between composition and morphology is not 
justifiable, because the inferred differences in rheology 
may be (and probably are) due to differences in eruption 
temperatures and crystallinity of the cryomagmas as well 
as differences in composition. In addition, pressure is 
an important variable affecting liquid compositions and 
cryovolcanic rheologies (Cynn et ctl. 1989: Hogenboom 
and Kargel, 1990), although we have not explicitly consid- 
ered such effects to date. Beyond these uncertainties, we 
find it intriguing that the more fluid cryolavas were erupted 
on the satellites of Jupiter and Saturn, while the more 
viscous cryolavas were erupted on the satellites of Uranus 
and Neptune. 

5. SUMMARY 

Cryovolcanic deposits have been observed on many of 
the icy satellites. We have gathered relevant viscosity 
data from the literature and have performed new viscosity 
measurements on a number of theoretically plausible cry- 
omagma compositions. The compositions included binary 
and ternary mixtures of water, ammonia, and methanol, 
both as liquids and as partially crystallized slushes. The 
viscosities of one-component liquids are orders of magni- 
tude too fluid to account for many of the flows, particularly 
those any with discernable thickness. Binary mixtures 
of liquid water and ammonia have substantially higher 
viscosities, up to 100 poise, enough to account for many 
of the observed flows. The thickest observed flows require 
the much higher viscosities that are only reached in our 
experiments by ternary mixtures or partially congealed 
binary mixtures. Collectively, these materials can account 
for all observable flows: thus, we find no justification for 
postulating solid-state extrusion. 

These experiments point out the importance of mixtures 
of simple compounds that provide a difference in tempera- 
ture between the solidus and liquidus in allowing the com- 
plex behavior of real extrusive materials, whether on the 
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terrestrial planets or on the icy satellites. Our experiments 
are still only an initial foray into the complexity of materi- 
als and material properties that almost certainly exists 
on the icy satellites, complexity that must he taken into 
account as more detailed models of their structure and 
history are constructed. 
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Cryovolcanism is among the foremost processes responsible for 
modifying the surfaces of icy satellites. Volcanic brine petrogenesis 
in ammonia-deficient satellites should mainly involve eutectic melt- 
ing in relevant salt-water systems. Carbonaceous chondrites pro- 
vide useful insights into the compositions of salts and aqueously 
altered rock in icy satellites and asteroids. Cl chondrites contain 
about one-fifth by mass of salts in various states of hydration. Many 
aspects of the petrogenesis and physical volcanology of icy satel- 
lite brines should be well described in the system H 2 0- 
MgS0 4 -Na 2 S0 4 . Minor components include sulfates of K, Ni, Mn, 
and Ca. Chondrites also contain abundant carbonates, but these 
are probably not very important in brine magmatism due to their 
low solubilities under expected conditions. Chlorides are also unim- 
portant under most circumstances because of the low cosmic abun- 
dance ratio of Cl/S. Soluble salts may have profound effects on the 
geology and structure of icy satellites and asteroids. In some models 
late episodes of water volcanism are facilitated by high buoyant 
forces due to the relatively high densities of sulfate-rich mantle and 
crustal layers. In other models early hypersaline brine volcanism 
quickly yields to plutonic magmatism due to low crustal densities. 
Europa probably has a layered crust composed of anhydrous 
Mg-Na sulfates near the base and a frozen or partially molten 
eutectic mixture of ice and hydrated Mg and Na sulfates near the 
surface. Ganymede may have a crust about 300 km thick composed 
of a 10:1 ratio of ice : mirabilite, and a mantle 500 km thick 
composed of 50% ice phases plus 50% hydrated Mg and Na 

Sulfates. © 1991 Academic Press, Inc. 


1. INTRODUCTION 

Many icy satellites have dynamic geologic histories in- 
cluding episodes of widespread cryovolcanism and exten- 
sive internal chemical differentiation. Cryovolcanism 
first predicted by Lewis (1971), is among the most im- 
portant processes documented by the Voyager project. 
Most volcanic plains on icy satellites apparently were 
extruded as aqueous solutions formed by partial melting. 
The compositions of aqueous “lavas” are likely to fall into 
two broad classes: relatively high-temperature salt-water 
solutions and truly cryogenic mixtures of ammonia, wa- 
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ter, and other substances. This work is devoted to 
sentation of the properties and probable composit 
of ammonia-free salt-water solutions, and on selee 
planetological applications of these data. 

Aqueous brines rather than ammonia-water liqu 
would have been generated wherever ammonia was , 
sent or deficient, and temperatures in the satellite's 
asteroid’s interior attained the eutectic melting poin- 
a relevant salt-water mixture. NH 3 would be abse 
regions of the Solar System where temperatures ir. 
solar or circumplanetary nebula remained too warm : 
ammonia hydrates to condense; where nebular nitrog 
was entirely in the form of N : rather than NH 3 (Prinn a 
Fegley 1988); where some NH 3 condensed buuhe amou 
was stoichiometricaily deficient with respect to high 
reactive substances such as CO : (yielding ammonium c: 
bonate) and H : CO (yielding hexamethylenetetramir. 
(Kargel 1990); or where prior melt extraction episo L 
had left the interior barren of ammonia hydrate and v. i: 
excess water ice and soluble salts. One or more of the- 
situations could have prevailed for any given satellite v 
asteroid, to the best of our knowledge. However, NH 
free brine volcanism is most likely to have occurred o 
large asteroids and the icy Galilean satellites of Jupiter 
temperatures in the nebula where these objects forme 
probably exceeded those required for ammonia hydrat 
condensation (Lewis 1974, Cameron 1978, Prinn and F ■£ 
ley 1981). 

Salt-water cryomagmatism may occur by either of tv 
paths, depending on whether salts were initially present 
If the rock component resembled type Cl or C2 carbona 
ceous chondrites then it would initially have container, 
abundant water-soluble substances. If, however, thes<. 
objects and their precursor material had never seen aque 
ous chemistry prior to or during accretion then these ob 
jects would have consisted of ice intermixed with rock 
resembling C3 or ordinary chondrites; as internal temper- 
atures attained or approached the melting point of ice 
aqueous reactions would have driven the assemblage to- 
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ward equilibrium, possibly yielding hydrated phyllosili- 
cates and soluble salts. Thus, regardless of whether salts 
were initially present, partial melting probably would 
yield brines rather than pure water. The composition and 
sal: ty of a particular brine would depend on the oxida- 
tion state of sulfur in the rock, solution pH. and other 
factors. 

This contribution will: (1) constrain the probable com- 
positions of brines generated in chondritic asteroids and 
icy satellites; (2) present an objective rationale simplifying 
these compositions to a minimum of components; (3) ex- 
amine some physical properties of simple salt-water mix- 
tures; and (4) discuss possible planetological conse- 
quences of salt-water magmatism. 

Tiis paper is restricted to physical chemistry data ap- 
propriate to I atm due to a paucity of relevant data taken 
at elevated pressures. The planetological applications at 
this juncture are exploratory and qualitative. The author 
and co-workers have recently initiated experimental stud- 
ies of high-pressure phase equilibria in salt-water systems 
necessary for more quantitatively realistic planetological 
applications (Hogenboom et al. 1991). 

Previous publications from this laboratory related to 
cr\ volcanism have dealt mainly with the densities and 
viscosities of ammonia-water liquids (Croft et al. 1988. 
Kargel et al. 1991). A review of published and original 
physical chemistry data for a broad range of plausible 
crovolcanic liquids and general applications to the geology 
of icy satellites were presented in the author’s dissertation 
(Kargel 1990). 

The compositions of volcanic brines may vary within 
broad limits. These compositions dictate the temperatures 
of partial melting and the likelihood that melting can occur 
in convecting satellites. These compositions also deter- 
mine the quantities of brine and the thicknesses of salty 
crusts and mantles that might form, the surface composi- 
tions of icy satellites, the physical properties of erupted 
brines, and the morphologic characteristics of landforms 
and landscapes produced by brine volcanism. 

2. IMPLICATIONS OF CHONDRITE SALTS FOR 

EXTRATERRESTRIAL BRINE COMPOSITIONS 

Aqueous alteration of carbonaceous chondrites. Car- 
bonaceous chondrites offer valuable insights into the na- 
ture of low-temperature water-rock chemical interactions 
*n chondritic asteroids, and the composition of water- 
soluble material produced during such alteration events. 
Aqueous alteration in carbonaceous chondrites is shown 
by a great abundance of aqueously precipitated salts, of 
secondary phyllosilicates including montmorillonite, sap- 
°n:te, chlorite, and serpentine, of certain sulfides and 
hydroxides including tochilinite and brucite, and of spher- 
“litic magnetite (DuFresne and Anders 1962, Bostrom 


and Fredriksson 1964, Kerridge el al. 1979, Zolensky and 
McSween 1988. Bunch and Chang 1980, Barber 1985, 
Keller and Buseck 1990). In a general sense, most salts 
and alteration products in chondrites are similar to the 
phases produced by chemical weathering, diagenesis, and 
low-grade metamorphism of common terrestrial rocks. 
Since primary phases in chondrites are also generally simi- 
lar to terrestrial crustal phases, including olivine, pyrox- 
enes, feldspars, and sulfides, one reasonably infers that 
many of the aqueous reactions in carbonaceous chondrite 
parent bodies were probably similar to common terrestrial 
weathering and low-grade metamorphic reactions involv- 
ing combinations of oxidation, carbonation, hydrolysis, 
and hydration. 

Formation of secondary epsomite, magnetite, and ser- 
pentine can be represented by a composite reaction in- 
volving the oxidation of troilite and hydrolysis of forsterite 
in an aqueous medium: 

12Mg ; Si0 4 + 6FeS + 54H ; 0 - 130 ; — 

6[Mg ; Si : 0 3 (0H) 4 ] + 6[MgS0 4 ■ 7H ; 0] + 2Fe 3 0 4 (1) 

Planetologists commonly assume that the rock in icy 
satellites closely resembles carbonaceous chondrites 
chemically and mineralogically. Quite incompatibly, an* 
other common assumption is that the rocky fraction of icy 
satellites does not participate in low-temperature mag- 
matic processes except as a stable inert component of the 
residuum of partial melting. While few would argue that 
the first assumption is unlikely, the second assumption, 
that rock is inert, is generally invalid. Either the rock 
contains abundant water-soluble salts if it resembles type 
1 or 2 carbonaceous chondrites, or, if it more closely 
resembles type 3 carbonaceous chondrites or unequili- 
brated ordinary chondrites, then the rock contains unsta- 
ble anhydrous minerals that normally would break down 
and yield soluble components upon reaction with water. 
As shown in Fig. la, the quantities of soluble matter 
present in extensively altered type Cl chondrites and less 
altered C2 chondrites are considerable. Commensurate 
with these abundances are the effects that the salts would 
have on the igneous and geological evolution of differenti- 
ated icy satellites and asteroids. 

Theoretical modeling of cryomagma generation and 
evolution in salty icy satellites and asteroids is much more 
complex than in the simplistic salt-free case, since the 
number of chemical components, and free parameters, 
becomes very large. However, several factors justify con- 
sideration of simplified chemical systems. 

Salts dominated by MgSO A . Salts in Cl and C2 chon- 
drites are dominated by metal sulfates (Fig. lb), particu- 
larly the hydrated magnesium sulfate, epsomite (MgS0 4 ■ 
7H 2 0). Other important salt components include sulfates 
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NajS0 4 All other sulfates 

FIG. 1. Compositions of chondritic meteorites (derived from Fredriksson and Kenidge 1988, Burgess el al. 1991, and Dodd 1981). (a) Relatri 
mass proportions of salts (dehydrated sulfates and carbonates), aqueously altered rock components (mostly phvllosilicates and hydroxides), ar. 
unaltered anhydrous minerals (including olivine, pyroxene, sulfides, and metal), (b) Relative mass proportions of sulfates, carbonates, and chloridt. 
(assuming all chlorine occurs as chlorides). Evaporated seawater for reference, (c) Relative mass proportions of magnesium sulfate, sodium sutfau 
and all other sulfates combined. The composition of all highly soluble material, ignoring the chlorides, is given by the point “without CaS0 4 " i' 
(c). 


of Ca and Na (Fig. lc), and carbonates of Ca, Mg, and 
Fe. Minor components include sulfates of Ni, Mn, and K. 
Chlorides have not been widely reported in chondrites, a 
fact attributed to a cosmic abundance ratio Cl/S — 0.01, 
and, according to one reviewer, possibly also reflecting 
the relative stability of Cl-rich chlorapatite and sodalite 
during low-temperature aqueous alteration. MgC0 3 and 
BaS0 4 have been reported from interplanetary dust parti- 
cles. (Reitmeijer 1990, Richardson 1978, Bostrom and 
Fredriksson 1964, DuFresne and Anders 1962, Nagy and 
Andersen 1964, Fredriksson and Kerridge 1988, Burgess 
etal. 1991). 


The relatively large fraction of sulfates, especially mag 
nesium sulfate, and the exceedingly low abundance c 
chlorine in chondrite salts contrast sharply with most ter 
restrial brines and evaporite deposits (Fig. lb). Terrestria 
salts are most often dominated by NaCl and are relativel; 
deficient in MgS0 4 . This distinction is easily explained 
The Earth’s continental upper crust, from which mos 
terrestrial brines ultimately derive their solutes, has bee: 
greatly enriched in Na and depleted in Mg by multistag e 
igneous fractionations, and was derived from a maniU 
that had previously seen severe chalcophile elemer. 
depletions apparently related to core formation. Na is 
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factor of 6 less abundant and magnesium a factor of 7 
m0 re abundant in Cl chondrites relative to the average 
continental upper crust (Anders and Grevesse 1989, Tay- 
lor and McLennan 1985). Sulfur. 40% of which is oxidized 
in the form of sulfates in the Cl chondrite Orgueil (Fre- 
d; ,sson and Kerridge 1988). is well over an order of 
magnitude more abundant in chondrites than in the 
Earth’s upper crust, while the condritic abundance of 
chlorine is similar to its abundance in the upper crust. 

One could usefully consider scenarios where S0 4 or Mg 
are unimportant in chondritic brines, or where Fe is as 
abundant as magnesium. For example, chemically reduc- 
ing conditions during aqueous alteration could stabilize S 
in the form of relatively insoluble metal sulfides. Reducing 
cc editions could also allow high concentrations of dis- 
solved Fe" 2 ; in view of the chondritic ratio Fe/Mg ~ 2 
(Anders and Grevesse 1989), dissolved Fe"* easily could 
be as important as Mg under these conditions. Strongly 
alkaline conditions, perhaps caused by a high NHj con- 
tent, could sharply curtail hydrolysis, thus slowing the 
rates at which major cations enter solution; furthermore, 
sufficiently alkaline conditions would stabilize Fe and Mg 
in the form of solid Fe- and Mg-hydroxides (Brookins 
1 1 preventing significant concentrations of these met- 
als in aqueous solution. Strongly acidic conditions could 
favor significant solution of carbonates. While these are 
not a priori unlikely situations, the author maintains a 
more conservative approach is assuming the prevalence 
of conditions and mineralogies represented by available 
aqueously altered cosmic material (i.e., carbonaceous 
chondrites). Thus, a substantial fraction (40%) of S is in 
the form of sulfate, particularly highly soluble MgS0 4 and 
Na,S0 4 ; carbonates and CaS0 4 are nearly insoluble; and 
the abundances of dissolved Cl and Fe * are negligible. 

Over one-fifth of the mass of the Cl chondrite Orgueil 
consists of salts in various states of hydration (Fig. 1). 
Calculated entirely on an anhydrous basis, salts constitute 
14.6% by mass of this meteorite, of which about two- 
thirds of the salts (by mass) are highly water soluble under 
normal conditions. MgS0 4 , constituting nearly three- 
fourths of the highly water-soluble fraction, is clearly the 
single most important salt in chondritic brines, followed 
by Na 2 S0 4 . 

Minor and weakly soluble salt components . MgS0 4 
and Na,S0 4 together account for 97% of the highly water- 
soluble' fraction of Orgueil. The remaining 3% consists 
primarily of the sulfates of K, Mn, and Ni. As shown in 
Fig. lb Ca-Mg-Fe carbonates are important solid compo- 
nents of chondrites, but are normally nearly insoluble, 
and could not normally constitute important components 
of magmatic brines. The same is true, to a slightly smaller 
degree, of CaS0 4 . The high abundance of carbonates and 
gypsum in carbonaceous chondrites is a result of localized 


and essentially isochemical aqueous alteration and in situ 
chemical reprecipitation. Therefore, these phases would 
remain with the solid rocky residue during anatexis and 
melt removal. Brine magmatism, as considered in this 
report, is a fundamentally different process in that it in- 
volves open-system transport of highly soluble compo- 
nents, resulting in gross chemical fractionations on large 
spatial scales. 


3. PHASE EQUILIBRIA 

Aqueous binary sulfates. Figure 2 shows solid-liquid 
phase equilibria in binary aqueous systems involving 
MgS0 4 , Na ; S0 4 , CaS0 4 , K : S0 4 , MnS0 4 , and NiS0 4 , all 
at 1 atm pressure. These binaries include all the major 
and minor sulfate components reported in carbonaceous 
chondrites, listed in order of abundance (Fredriksson and 
Kerridge 1988). The system H 2 0-MgS0 4 is by far the 
single most important of these binaries. Under unex- 
pected conditions, not examined in this work, aqueous 
FeS0 4 or (NH 4 ) 2 S0 4 could also be important components , 
and might yield very different types of petrological evolu- 
tion than considered below. 

The phase equilibria shown in Fig. 2 probably should 
remain valid, as a first approximation, to pressures up to 
2 kbars, since none of the substances are extremely vola- 
tile or unusually compressible. Higher pressures should 
progressively introduce qualitative changes, particularly 
solid-state phase changes. Kargel (1990), drawing on pub- 
lished data, showed that solid-liquid equilibria in the sys- 
tem H 2 0-NaCl remain qualitatively unchanged, except 
for the stable ice polymorph, under pressures ranging 
from 1 atm to 12 kbars. The system H : 0-MgS0 4 currently 
is being studied experimentally under pressures up to 4 
kbars; early results (Hogenboom et al. 1991) indicate no 
major surprises. However, subtle pressure-dependent as- 
pects may have significant planetological consequences. 
These results will be presented elsewhere. 

Metastable magmatic phenomena often occur in nature, 
particularly in quenched lavas. Supercooling, glass forma- 
tion, and precipitation of metastable sulfate solids occur 
regularly under laboratory conditions, and indeed are 
avoidable only with considerable patience. Metastable 
phenomena, analogous to those which occur in the labora- 
tory and in silicate volcanic flows, are especially likely 
during rapid solidification of the surfaces of aqueous sul- 
fate flows. Persistence of these metastable phases over 
geologically lengthy periods is possible, a speculation sup- 
ported by the occurrence of metastable magnesium and 
calcium sulfate minerals in terrestrial evaporite deposits. 
More often, terrestrial evaporite minerals revert to their 
stable forms when held at normal temperatures for a few 
hours or days (the surfaces of icy satellites are not "nor- 
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mal, however). Solidification of the interiors of extrater- 
resina 1 aqueous sulfate flows would likely occur over 
sufficiently long intervals (hours to days) that equilibrium 
sc J assemblages probably would be produced. More 
■r jrtantly, melting processes in asteroid and satelh e 
imenors would occur over very lengthy periods and major 
departures from equilibrium are unlikely. Thus, despite 
the likelihood that flow surfaces would be metastable^an 
aspect with possible spectroscopic significance), most 

SlTp“.‘ ’ robably wou,d involve 

The recognition that MgS0 4 is probably the most im- 
po-.ant solute in extraterrestrial volcanic brines invites 
e- Tiinat'on of multicomponent phase equilibria when 

MaSO Ch T n h d :: te , SU,fate C ° mponents are added to aqueous 

Fie ?Vnd h a nary phase equilibna are illustrated in 
Fig. 3 and are discussed below. 
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FIG. 3. Ternary aqueous phase equilibria involving third components added to the system H ; 0-MgS0 4 . (a) H ; 0-MgS0 4 -Na ; S0 4 . (b) H,0- 
MgSo' 4 -CaS0 4 . (c) H-0-MgS0 4 -K 2 S0 4 . (d) H ; 0-MgS0 4 -MnS0 4 . (e) H,0-MgSO 4 -NiSO 4 . Data sources: Binary phase diagrams and source- 
listed in Fig. 2 caption- Levi (1923), Blasdale (1927), Starrs and Storch (1930), Starrs and Clarke (1930). Benrath and Blankenstein ( 1934), Benrath 
and Neumann (1939), Soboleva (1960). 


cal valence states of Mg and Mn and their roughly similar 
ionic radii. Among the hepta- and hexahydrates, the solid 
phases have higher contents of Mg, lower concentrations 
of Mn, and higher Mg/Mn ratios than coexisting liquids. 
In silicate petrological vernacular, Mn is moderately in- 
compatible. 

System H 2 0-\igS0 2 -NiS0 4 . The chondritic ratio 
Ni/Mg ~ 0.12 (Anders and Grevesse 1989), nearly two 
orders of magnitude greater than in Earth’s upper crust 
(Taylor and McLennan 1985), allowing the possibility that 
Ni could be an important constituent of cryovolcanic 
brines. Indeed, Mg-Ni-Na sulfates in Orgueil include 
up to 28% NiS0 4 , although on average NiS0 4 is minor 
(Fredriksson and Kerridge 1988). Figure 2f shows binary 
phase relations of aqueous NiS0 4 . Comparison with the 
H 2 0-MgS0 4 system shows that sulfate solubility curves 
and the series of hydrated solids are very similar, except 
that NiS0 4 hexahydrate forms crystallographically dis- 
tinct monoclinic and tetragonal polymorphs, whereas 
MgS0 4 forms only the monoclinic variety. The apparent 


lack of a Ni counterpart to MgS0 4 • 12H : 0 might reflect 
kinetic problems in forming Ni dodecahydrate. rather than 
its lack of existence, as Mg dodecahydrate is also more 
difficult to form than heptahydrate. 

The hydrates of MgS0 4 and NiS0 4 form complete solid 
solutions, analogous to the behavior of Mg and Ni in 
silicates. This is not unexpected since the valence states 
of Ni and Mg are the same and their ionic radii are near’} 
identical. The literature sources used in the construction 
of Fig. 3e did not find Mg dodecahydrate, so the three- 
phase cotectic between ice and sulfate involves the hepta- 
hydrate. This cotectic is a simple thermal trough run- 
ning almost isothermally across the phase diagram from 
one binary to the other. The stability field of tetragonal 
NiS0 4 • 6H 2 0 terminates at a peritectic upon the addition 
of just 4% MgS0 4 to the binary H : 0-NiS0 4 eutectic liq 
uid. Therefore, except for this small tetragonal hexah; 
drate field occurring only for very Ni-rich composition: 
the aqueous MgS0 4 and NiS0 4 systems, and ternary mix 
tures of them, are isomorphic. 

Two sets of representative three-phase triangles in Fig 
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FIG. 3 — Continued 


3e connect the compositions of solid sulfates with their 
equilibrium liquids on the sulfate-ice cotectic. These data 
indicate that the partitioning behavior of Ni and Mg among 
the sulfates and liquid are very similar and of order unity, 
in accordance with these metals' chemical similarity. In 
other words, the ternary phase equilibria are insensitive 
to the Mg/Ni ratio of the solution. Hence, one may treat 
ternary mixtures in this system as pseudobinary mixtures, 
with little sacrifice in accuracy, by combining MgS0 4 + 
NiS0 4 on a molar basis into a single component. There- 
fore, the chondrite salt vein compositions reported by 
Fredriksson and Kerridge (1988), which strictly must be 
described with no less than four major components (H ; 0 
and sulfates of Mg, Ni, and Na), can be reduced to three 
components. Examined in more detail, the Ni/Mg ratios 
of solid solutions are always greater than in the coexisting 
liquid, so Ni can be considered a compatible element. 


Synthesis . The physical chemistry of carbonaceous 
chondrite-equilibrated sulfate brines may be roughly ap- 
proximated in the system H 2 0-MgS0 4 (Fig. 2a), and more 
accurately in the ternary system H 20 -MgS 0 4 -Na 2 S 0 4 
(Fig. 3a). Minor sulfate components (sulfates of Ni, Mn, 
and K) usually would not have major effects on the liquid's 
physical properties and major phase equilibria. However, 
it is possible to make a rough accounting for these compo- 
nents as long as their concentrations are low. With little 
loss in accuracy one may sum K 2 S0 4 and Na 2 S0 4 into a 
single component and MgS0 4 , NiS0 4 , and MnS0 4 into 
another component. However, as terrestrial experience 
shows very well, complex igneous fractionations can 
cause initially minor components to become major ones 
in strongly differentiated materials. In this event, quater- 
nary and higher-order chemical systems may be required 
for a more accurate accounting. 
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FIG. 4. Densities of aqueous salt solutions as functions of concentra- 
tion. Eutectics denoted by solid circles. Sources in text. 


CaS0 4 and the carbonates of Ca, Mg, and Fe are spar- 
ingly soluble under a broad range of environmental condi- 
tions. While these salts are important constituents of car- 
bonaceous chondrites, they were formed by largely in situ 
processes of solution and precipitation, and would not 
have constituted an important fraction of the ionic load in 
chondritic hypersaline brines. 

The chondritic ratio of Cl/S -0.01, sharply limiting 
the potential volumes of low-temperature chloride brines. 
Taking Ganymede as an example, eutectic melting at 
— 52 C in the system MgCL— CaCL— H 2 0, involving all the 
chlorine present, could yield 0.22 mass percent of eutectic 
brine, equivalent to 0.37 voI% of Ganymede (assuming 
Ganymede contains 60% by mass of chondritic rock 
containing 698 ppm Cl). Although this is not a large 
quantity, it is capable of forming a crust 3 km thick if 
completely extruded. If nothing else is capable of melt- 
ing, perhaps because convective adiabats are cooler 
than the melting points of aqueous sulfates, then chloride 
brine volcanism could dominate surface geology and com- 
position. 

One reviewer suggested that low-temperature petro- 
logic evolution of ice-bearing systems containing both 
chlorides and sulfates could be very complex. Alterna- 
tively, as another reviewer suggested, chlorine could re- 
main immobilized in stable, nearly insoluble phases such 
as chlorapatite and sodalite until temperatures in the range 
of hundreds of degrees Celsius are attained. In the latter 
case, late-stage eruptions of hydrothermal chloride brines 
may occur during dehydration of phyllosilicates. 

4. DENSITIES AND VISCOSITIES OF BRINES 

Densities. Published sources give the densities of 
aqueous salt solutions as functions of composition and 
temperature (CRC Handbook of chemistry and Physics, 
Fabuss et al. 1966, Korosi and Fabuss 1968, Chen et al. 
1980). The data collectively cover concentrations up to 
the aqueous salt eutectics, and extend over a wide range 
of temperatures above 293 K. Figure 4 gives the densities 


of planetologically important aqueous salts as function^ 
of concentration for 293 K and 1 bar. These data can be 
extrapolated and interpolated over data gaps to obtain the 
eutectic liquid densities at the eutectic temperatures. The 
eutectic liquid in the system H 2 0-MgS0 4 (17 
MgS0 4 ) at the eutectic temperature ( -3.9°C) has a den- 
sity of about 1.186 g cm' 3 . 

The densities of frozen eutectic mixtures of ice and 
salts are also of planetological interest, particularly in 
the systems involving H : 0, MgS0 4 , and Na : S0 4 . The 
densities of most solid sulfate phases of interest have been 
published (CRC Handbook), but the density of MgS0 4 • 
12H,0 has not. A preliminary measurement of this den- 
sity, 1.51 ± 0.01 g cm' 3 , was presented orally by Hog n- 
boom et ai. (1991). The densities of solid magnesium and 
sodium sulfates are shown graphically as a function of 
water content in Fig. 5; to the level of precision presented 
in this diagram these densities are not significantly temper- 
ature-dependent. 

Figure 6 summarizes the densities of selected liquids 
and solids of cryovolcanic interest. All liquids in Fig. 
6 are buoyant with respect to rock-ice mixtures, thus 
encouraging early aqueous volcanism. 

A frozen eutectic mixture of 47% ice and 53% 
MgS0 4 - 12H ; 0 has a density of about 1.126gcm~ 3 , about 
5% lower than the density of the coexisting liquid. A 
frozen eutectic mixture in the system H,0-MgS0 4 - 
Na 2 S0 4 has a density — 1. 133 g cm' 3 . Thus, the solid floats 
on the liquid, presenting the same buoyancy “problem” 
as occurs for water and ice. Once a thick volcanic crust 
forms subsequent volcanism is inhibited by the negative 
buoyancy of the late-arriving liquids relative to that of me 
solid crust (Croft 1985). Adapting the suggestion of Croft 
et al. (1988) for ammonia- water magmatism. early 
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FIG. 6. Densities of selected solids and liquids of volcanological 
interest. After Kargel (1990). 


salt-water volcanism may frequently have yielded to later 
plutonism. Under certain circumstances plutonism might 
generate surface tectonic features such as domical uplifts 
or patterns of parallel, radial, or concentric fractures or 
dike swarms. 

Despite the apparent “buoyancy crisis,” morphologic 
evidence exists for widespread late (i.e., post heavy bom- 
bardment) volcanism on several icy satellites. One solu- 
tion to the buoyancy crisis is apparent upon examining 
Fig. 6. Pure water, dilute brines, and eutectic sodium 
sulfate solutions are less dense than a frozen ternary eu- 
tectic assemblage; hence, these less saline solutions could 
readily extrude. Section 5.1 presents several models 
where early crust-forming ternary eutectic melting events 
exhaust the supply of magnesium sulfate in the region 
of partial melting, yielding to later production of binary 
eutectic sodium sulfate solution or nearly pure water. 
Thus, while some situations might favor an early perma- 
nent transition from hypersaline volcanism to plutonism, 
other plausible conditions would allow a late-stage re- 
sumption of volcanism by extrusions of water or sodic 
brine. 


Viscosities. The viscosity of magma is encountered 
in every problem of magma migration and lava flowage. 
Kargel et al. (1991), in an experimental study of ammo- 
nia-water liquids and crystal-liquid slurries, highlighted 
the important effects of partial crystallization on rheology. 
The relative differences in lava rheologies are well indi- 
cated by comparisons of liquid viscosities at their respec- 
tive freezing points in the crystal-free state. These data 
should not be directly utilized in rheological modeling of 
lava flows since they do not take into account the im- 
portant effects of partial crystallization. Rather, these 
data are useful for determining geologically familiar rheo- 
logical analogs (e.g., silicate lavas). The undaunted physi- 
cal modeler might choose to utilize these data after 
applying an empirical factor accounting for crystallization 
(Greenberg et al. 1991, Schenk 1990), or by employing 
an empirical function relating the relative viscosity and 
degree of crystallinity (Kargel et al. 1991). 

Available viscosity data with appropriate extrapola- 
tions are shown in Fig. 7 for aqueous magnesium sulfate. 
The eutectic H 2 0-MgS0 4 liquid has a viscosity ~ 0.06 P. 
The eutectic liquid in the system H 2 0-MgS0 4 -Na : S0 4 
has a viscosity -0.10 P. Aqueous eutectic NaCl, MgCU, 



FIG. 7. Viscosities of magnesium sulfate solutions. Sources: Korosi 
and Fabuss (1968) (solid circles), Novikova (1957) (open squares). 
Dashed sections where extrapolated. From Kargel (1990). 
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FIG. 8. Viscosities of volcanoiogically interesting liquids at their 
respective freezing points. Aqueous substances represent eutectic com- 
positions, the viscosities of which are from Kargel et at . (1991) and 
sources referenced in the text. Silicate data are based on published 
viscometry of natural and synthetic silicate liquids (Shaw 1969, Urbain 
et at. 1982, McBimey and Murase 1984, Murase et al 1985, Tauber and 
Arndt 1986. Ryerson et aL 1986. Hummel and Arndt 1985, Scarfe et at, 
1983). The melting points of silicates are based on experimental phase 
equilibria (dry) from Cox et aL 1979 and Morse 1985. The quantity 
(MgO -f Fe0)/Si02, in weight percent, is a measure of the degree of 
differentiation of silicate liquids. Strongly silica-undersaturated alkaline 
liquids do not lie on the curve. 


and CaCl 2 indicate viscosities of, respectively, 0.05, 0.2, 
and 0.4 P. Multicomponent aqueous chlorides may have 
somewhat higher viscosities, but probably not over 1 P. 

Figure 8 summarizes the crystal-free viscosities of cryo- 
volcanic and silicate liquids at their respective freezing 
points. The aqueous mixtures span a range of viscosities 
comparable to the range exhibited by silicate lavas. 
Salt-water eutectic mixtures are one to two orders of 
magnitude more viscous than water, but are two to three 
orders of magnitude less viscous than basalt and ammo- 
nia-water eutectic liquid. This physical contrast probably 
would be reflected in important differences, along compo- 
sitional lines, in volcanological behavior. 

The low viscosities of salt brines would encourage melt 
segregation from residual crystalline matrices at very low 
melt fractions (McKenzie 1985), favoring fractional partial 
melting over equilibrium batch melting, and fractional 


crystallization in preference to equilibrium crystalliza- 
tion. Fractional melting involves removal of liquid in small 
increments as melting proceeds, and may be approxi- 
mated as continuous removal of ideal eutectic composi- 
tions interrupted when one or another solid phase is ex- 
hausted in the solid residue. Melting resumes when tae 
next eutectic temperature is attained in the simpler svstem 
lacking the exhausted phase. In equilibrium batch melting 
the liquid is extracted in one or a few large increments, 
typically allowing the liquid composition to migrate up a 
cotectic after a solid phase is exhausted (see Cox et al 
1979). 

5. APPLICATIONS 
5.1. Structural Models 

Construction. This section examines possible internal 
chemical structures resulting from the differentiation of a 
series of generic satellites and asteroids in the system 
H 2 0-MgS0 4 -NaS0 4 , plus inert rock (Fig. 3a). The sup- 
ply and transfer of heat is not considered, but rather is 
assumed to be sufficient to induce complete differentia- 
tion. The deep internal pressures in large asteroids, in icy 
satellites the size of Tethys, and in the icy crust of Europa 
are on the order of a kilobar. This is a sufficiently low 
pressure that Fig. 3a and satellite structural models based 
on Fig. 3a, even though strictly valid for 1 atm, are proba- 
bly suitable approximations. However, the interior pres- 
sures in Ganymede and other large icy satellites attain up 
to 50 kbars, and qualitatively important effects on phase 
equilibria are expected. Even so, the models presented 
here are useful in advertising the probable importance of 
salts, particularly sulfates, in the low-temperature differ- 
entiation and geologic evolution of these large objects. 
But high-pressure phase equilibria in relevant systems are 
sorely needed for more accurate modeling and reliable 
insights into the histories of the larger objects (Hogen- 
boom et al. 1991 and in preparation). 

Structural models in- Fig. 9 were constructed for five 
initial rock fractions varying from 1.0 (pure chondritic 
rock) to 0.4 by mass. The percentage of sulfates in the 
rock component was calculated using the average S abun- 
dance in Orgueil, 6.25 mass% (Anders and Grevesse 
1989), and the average fraction of this S occurring as 
sulfates, 0.40 (Fredriksson and Kerridge 1988, Table 1); 
these two values yield the fraction of sulfate sulfur (S -6 ) 
in Orgueil: 2.50% by mass. This amount of oxidized S was 
apportioned among the elements Ca, Mg. and Na in the 
ratios measured by Fredriksson and Kerridge (1988, Table 
2), with the minor approximation that NiS0 4 and MnS0 4 
in Orgueil were summed along with MgS0 4 , and K,S0 4 
was summed with Na 2 S0 4 . The mass fractions of sulfates 
in the chondritic rock component, calculated on an anhv- 
drous basis, are MgS0 4 6.05%, Na 2 S0 4 2.07%, and CaS0 4 
1.79%. Only the first two components are highly water- 
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soluble, giving a total of 8.12% highly soluble salts. The 
calculated mass percentages of hydrated salts are epsom- 
ite (MgS0 4 • 7H 2 0) 8.80%, bloedite (MgS0 4 • Na ; S0 4 • 
4H.O) 4.87%, and gypsum (CaS0 4 • 2H ; 0) 2.27%. The 
quantity of water associated with the two highly soluble 
phases, epsomite and bloedite, sums to 5.55%. 

The exact meaning of “rock : ice ratio” is unclear when 
dealing with systems containing hydrated silicates and 
salts. Furthermore, the quantity of rock differs before and 
after leaching. A more precise term is needed, 

+ mj, 

wh . : e 7? , is the initial mass fraction of chondritic rock 
(including salts and water of hydration). m ri is the initial 
mass of chondritic rock, and is the mass of water 
added to the chondritic rock. The total fraction of water 
in an ice-chondrite mixture is greater than 1 - R t since 
the rock itself contains water. 7?, would correspond to the 
true rock fraction only if chondritic rock and ice are mixed 
without reequilibration or differentiation. Reequilibration 
without differentiation actually would yield more rock and 
les ice than indicated by 7? ( since epsomite and bloedite 
would take up water as they transformed to mirabilite 
and magnesium sulfate dodecahydrate. Low-temperature 
differentiation and leaching of soluble salts from the rock 
would result in less rock than 7 ? ( . The final mass fraction 
of rock, 7? f , after differentiation is 

R( = (1 - 0.1367)7?, = 0.8633 7?,. 

Mt ting and crystallization sequences can be followed in 
Fig. 3a after calculating the fractions of H 2 0, MgS0 4 , 
and Na 2 S0 4 . Fig. 2a (system H 2 0-MgS0 4 ) considerably 
assists the interpretation of petrologic evolution in the 
ternary system (Fig. 3a). The amount of water in each 
model satellite is calculated as the amount of water of 
hydration of epsomite and bloedite in the initial rock 
(5.55% of initial rock mass) plus the amount of added 
water. Phyllosilicate water is assumed to be stably bound 
to ‘.he rock. The amounts of the ternary components are 
resummed to 100%. Chondrite-water mixtures are shown 
in Fig. 3a as the straight line running from the composition 
of Orgueil to the water apex. 

Differentiation was assumed to progress by continuous 
fractional melting, resulting in a perfectly gravitationally 
stably stratified structure. Figure 9 shows the internal 
structures of the five differentiated model objects. The 
idealized evolution of each model is briefly described 
below. 

Model A. Figure 9a, with 7?, = 1.0, shows the differen- 
tiated structure of an object initially composed of pure 
chondritic rock with no added water. Epsomite and bloe- 
dite melt incongruently, yielding the peritectic liquid P2 


(Fig. 3a), containing over 30% MgS0 4 and 4% Na : S04. 
Liquid P2 is less dense than the rocky residue, so the 
liquid is readily erupted, forming a volcanic crust. Each 
successive volcanic flow would have an identical compo- 
sition until either epsomite or bloedite is exhausted in 
the object’s interior. P2 flows could solidify by fractional 
crystallization (where crystals gravitationally settle and 
segregate from the liquid almost as rapidly as they form) 
or by equilibrium crystallization (where crystals remain 
in contact and continuously reequilibrate with the cooling 
liquid). Solidification of most silicate flows approximates 
equilibrium crystallization. If liquid P2 volcanic flows cool 
analogously, the final homogenous crystalline mixture 
would be composed of 84% MgS0 4 ■ 12H 2 0 + 11% 
Na 2 S0 4 • 10H : O + 5% ice. However, as discussed in 
Section 4, the low viscosity of the hypersaline P2 brine 
might favor fractional crystallization. In this case each 
flow w ould consist of several layers of different composi- 
tions: the top and last layer to solidify would have the 
composition of ternary eutectic E. 

In the next stage of melting, the residue left by forma- 
tion of liquid P2 melts incongruently at peritectic P3 (Fig. 
3a). Liquid P3, not shown in Fig. 3a (but consult the binary 
approximation in Fig. 2a), is denser than the volcanic 
crust of composition P2, and probably would be unable 
to ascend to the surface. Rather, liquid P3 would pool 
near or intrude the base of the crust, eventually forming 
the upper mantle layer of the asteroid. 

The residue left from the second stage of melting then 
suffers a third stage of incongruent melting, yielding a 
small amount of liquid and a final solid residue composed 
of anhydrous sulfates. This stage of melting also yields a 
liquid substantially denser than the upper mantle, thus 
resulting in the formation of a lower mantle layer. The 
anhydrous sulfates probably would remain in the rocky 
core until water was released from the dehydration of 
gypsum and phyllosilicates. 

The average density of this object after differentiation 
is 2.69 g cm’ 3 , representing -2.8% volume expansion 
over the undifferentiated initial chondrite. Much of this 
expansion takes place after formation of the peritectic P 2 
crust. Since all liquids generated subsequently would tend 
to intrude beneath this early-formed crust, one expects 
extensional strains of ~-l%, sufficient to cause significant 
crustal rifting. 

Model B. Figure 9b shows the structure of a satellite 
with 7?j = 0.7. This object contains enough water that 
epsomite and bloedite contained in the initial chondritic 
rock would react with added water to form magnesium 
sulfate dodecahydrate and mirabilite, leaving excess wa- 
ter to form ice. This mixture first melts at the ternary 
eutectic E resulting in the formation of a thick volcanic 
layer (later to become the upper mantle). If each succes- 
sive volcanic flow cooled by equilibrium crystalliza- 
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FIG. 9. Structural models of differentiated objects having various initial rock fractions ranging from LO to 0 A Uft-hand colum g 
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the densities of each frozen layer. 


tion this layer would be composed of 50% ice + 44% 
MgS0 4 • 12H 2 0 + 6% Na 2 S0 4 • 10H 2 O; its density is 
1.133 g cm -3 (at 1 atm pressure). 

Fractional eutectic melting and extraction of liquid E 
(in Fig. 3a) from the interior would drive the solid residue 
toward the binary join H 2 0— Na 2 S0 4 . The subsequent 
evolution of this residue, composed of ice and mirabilite 
(plus rock), is best visualized in the binary H 2 0-Na 2 S0 4 
phase diagram (Fig. 2b). A second stage of melting could 
occur at the binary eutectic, producing a liquid containing 


4% Na 2 S0 4 and no MgS0 4 . This liquid has a density ot 
about 1.05 g cm -3 , substantially less than the density ot 
the frozen ternary eutectic layer (1.13 g cm -3 ). Therefore, 
this binary liquid could easily ascend to the surface, form 
ing a crust composed of 91% ice + 9% mirabilite, am 
leaving a residue in the rocky core composed of mirabi vS • 
Mirabilite may then melt incongruently at 32.5 S C, forming 
an anhydrous residue of thenardite and a liquid containing 
34% Na 2 S0 4 . The dense binary peritectic liquid wouU 
intrude at the base of the mantle, forming a frozen lowe 
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m antle composed of 76% mirabilite - 24% ice. The aver- 
age density of this satellite is 1.76 g cm . 

Model C. The initial rock fraction, R , = 0.6689, has 
bet chosen so that ternary eutectic melting leaves a 
rescue of ice and mirabilite in precisely the ratio of the 
binarv eutectic; thus, ternary eutectic melting is followed 
hv binary eutectic melting in the system H 2 0-Na 2 S0 4 , 
without leaving a residue of either salt or ice. R , is also 
very similar to Ganymede’s; Mueller and McKinnon 
(1 988) calculated R , = 0.64 to 0.65 when the rock 
was assumed to have a Cl chondrite composition and 
Ganvmede was' assumed to be differentiated. The differ- 
ent 7 'ion of this object is qualitatively the same as in 
mou i B except that there is no mirabilite residue and, 
thus.' no Na 2 S0 4 -rich lower mantle layers. The average 
density of this satellite is 1.68 g cm 3 . 

Model D. This model is again similar to models B and 
C except that water ice is left in the residual rocky core 
after mirabilite is exhausted during the binary phase of 
melting in the system H : 0-Na : S0 4 . The excess water ice 
melts at 273 K. yielding pure water with a density less 
tha the ternary eutectic mantle, but greater than the 
der. ,;v of the binary sodium sulfate-rich layer; hence, 
water may tend to intrude between these two salt-rich 
layers. However, when it freezes, water ice is less dense 
than the overlying layer composed of ice and mirabilite, 
so the intruding ice may tend to work its way to the surface 
overtime. In any case, the assumption that the final differ- 
entiated structures are gravitationally stably stratified 
means that the final structure indeed exhibits an ice crust. 
Exactly what really would happen is not obvious. Ice 
micct ascend slowly in small diapirs, rather like terrestrial 
salt domes, or overturn could be catastrophic. Alterna- 
tively, this water could dissolve substantial quantities of 
salt, becoming gravitationally stable at depth or mixing 
with the crust or mantle. 

The average density of this object is 1.54 g cm . it 
R = 0.6, this mixture is similar to Callisto’s composition 
as calculated by Mueller and McKinnon (1988) (although 
obviously the structure calculated here is different from 
the calculated by Mueller and McKinnon). 

Model E. This is qualitatively the same as model D, 
but the amount of water is greater, resulting in a thicker 
ice crust. The average density is 1.26 g cm' 3 , similar to 
the least dense, iciest satellites in the Solar System. 

Importance of Na#. Na is defined as Na 2 S0 4 / 
(Na 2 S0 4 + MgS0 4 ). The models in Fig. 9 assume that the 
initial rock has the composition of Orgueil, with Na# — 
0.255. However, Na# differs among various chondrites, 
anu even varies considerably within Orgueil itself, so Na# 

• should really be considered a free parameter. 

The models in Fig. 9 have Na# = 0.255, compared to 
Na# = 0.15 for the ternary eutectic liquid. Any line drawn 


from eutectic E through any of the ice-bearing mixtures 
(models B. C. D, and E) necessarily intersects the joint 
H ; 0-Na ; SO 4 (see Fig. 3a). Simply, but significantly, this 
means that ternary eutectic melting of any mixture of i ce 
and Orgueil-like rock will leave a residue containing H 2 0 
and Na 2 S0 4 . and not retaining MgS0 4 . Hence, MgS0 4 is 
always the first component to be exhausted in models 
B through E in Fig. 9, with the consequence that the 
differentiated structures always have a mirabilite-ice bi- 
nary eutectic layer overlying a ternary MgS0 4 -rich layer. 

If. however. Na# < 0.15 in the initial rock component 
then the ternary eutectic liquid has Na# greater than the 
initial rock: then ternary eutectic melting would leave a 
residue composed of H 2 0 and MgS0 4 and depleted in 
Na.S0 4 . Therefore, the petrologic evolution of objects 
having Na# < 0.15, and their resulting differentiated 
structures and geologic histories, would be very different 
from those depicted in Fig. 9. 

Fate and significance of minor components. Figure 9 
does not consider the possible petrologic roles of minor 
soluble and abundant but weakly soluble components. 
Consideration of petrologically evolved silicate planets 
such as the Earth and Moon, where minor mantle compo- 
nents often become major crustal constituents, is suffi- 
cient reason to briefly consider minor salts in the evolution 
of icy satellites and asteroids. 

Phase equilibria shown in Fig. 3 indicate that Ni is 
compatible while Mn and K can be considered incompati- 
ble with respect to Mg-rich sulfates. Therefore, early- 
generated liquids will be Ni-poor and K- and Mn-rich 
compared to the initial salts, thus enriching the solid resid- 
uum in Ni and depleting it in K and Mn. Continued melting 
will yield liquids progressively more enriched in Ni and 
depleted in K and Mn. Complex evolutionary scenarios 
could involve remelting prior melts. This could yield liq- 
uids so enriched in K and Mn that they become major 
components, and could grossly deplete Ni. It is not known 
whether any icy satellites or asteroids have such compli- 
cated geological histories, but a brief look at Europa and 
Ganymede suffices to allow this possibility. Magmatic 
redistribution of ^K-carrying salts such as K 2 S0 4 could 
be an important factor affecting the subsequent thermal 
structures and evolution of differentiated asteroids and 
icy satellites. 

K and Rb salts could eventually be useful in automated, 
in situ isotopic dating of cryovolcanic flows (Kargel 1989). 
K and Rb salts have yielded correct ages for terrestrial 
evaporites, although groundwater flow and subtle heating 
often disturbs parent-daughter systems (Brookins et al. 
1980, 1985). The frequent problems encountered in dating 
terrestrial evaporites would not be expected in deep-fro- 
zen cryovolcanic flows. 

CaS0 4 is expected to be a minor component under all 
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Na# = 0.255 and 0.083, clearly indicating that this mete- 
orite is grossly heterogeneous with respect to its sulfate 
vein compositions. Expressed as molar fractions these 
numbers are, respectively, 0.226 and 0.071. Despite first 
appearances, this heterogenity does not indicate fluid 
transport and salt/fluid fractionation on the spatial scale 
of Orgueil (decimeters). This is because each portion of 
Orgueil contains both veins and matrix, and what is pres- 
ent in the local veins is absent in the matrix so that each 
bulk sample of Orgueil (including veins and matrix) re- 
mains more or less solar in composition. 

Four individual vein analyses give Na#’s of 0.338 (Or- 
gueil vein #7), 0.336 (Orgueil vein #6), 0.335 (Ivuna vein 
#20), and 0.143 (Ivuna vein #19) (from Fredriksson and 
Kerridge 1988). In molar fractions these numbers are, 
respectively, 0.302. 0.301, 0.300, and 0. 124. Of these vein 
analyses, three (from two different meteorites) are virtu- 
ally identical in Na#, although Ni/Mg varies consider- 
ably. The other vein is similar to any of three eutectics, 
including the eutectic in equilibrium with ice, mirabilite, 
and epsomite (Fig. 3a). 

The three clustered Na#'s do not correspond to any 
known eutectic liquid in the aqueous Mg-Na sulfate sys- 
tem, nor do they correspond to any of the mixed Mg-Na 
sulfates known to occur. Fredriksson and Kerridge (1988) 
suggested a compositional similarity to nickel bloedite 
(MgS0 4 • Na ; S0 4 • 4H ; 0), but the stoichiometry is incor- 
rect for this or any other known sulfate. If bloedite is a 
major mineral then a consistent fraction of epsomite or 
another pure magnesium sulfate must be intermixed in 
these veins. The clustering of these three sulfates seems 
real, however. Since it is unlikely that three veins from 
two meteorites just happen to have the same composition, 
this clustering argues for the existence either of a pre- 
viously unknown mixed salt with the approximate stoichi- 
ometry (Mg, Ni) : Na 2 (S0 4 ) 3 • xH : 0, or of a previously 
unknown eutectic occurring near this composition. Ki- 
netic factors may have been responsible for precipitation 
of an unknown metastable mixed Mg-Na-sulfate com- 
pound (which may subsequently have reverted to bloedite 
and epsomite). Alternatively, it is possible that the pres- 
sure variable is responsible; perhaps at the modest pres- 
sures deep in the interiors of chondrite parent bodies an 
unknown high-pressure mixed sulfate precipitates with 
the stoichiometry inferred above, or perhaps high pres- 
sures result in a different invariant melting composition 
corresponding to the clustered vein compositions. The 
Ivuna #19 sulfate composition, corresponding closely to 
known eutectics, may represent melting or precipitation 
at lower pressures where phase equilibria are similar to 
the 1-atm case. In any event, an experimental study of 
polybaric phase relations in this system would be a valued 
endeavor, as would the acquisition of additional high- 
quality chondrite sulfate compositions. Unless the 


method of Burgess et al. (1991) can be significantly inv 
proved this author does not recommend the heating-oxi 
dation technique. 

Chemical data reported by Fredriksson and Kerridge 
(1988) indicate that in one Mg-Na-Ni sulfate vein, appar- 
ently the vein Ivuna #19, Na and Ni are concent; .tec 
together on one side of the vein. Such zoning pro! bly 
indicates chemical evolution of the brine from which thi> 
vein precipitated; hence, the liquid composition was not 
precisely that of a eutectic, at least not at the time the 
vein precipitated (the parent solution could be a high- 
pressure eutectic, however). Figure 3e tells us that any 
sulfates precipitated early would tend to have slightly 
higher Ni contents than the liquid, and precipitation would 
force the liquid to migrate in the direction of Ni deple on. 
The early liquidus phase assemblage was both Na- am Ni- 
rich. Candidate sodic phases are mirabilite and bloedite. 
However, even the most sodic portions of this vein con- 
tain more Mg than either of these phases, so the brine 
must have been precipitating along a cotectic involving 
coprecipitation of a magnesium sulfate in addition to a 
sodic phase. In the region of Fig. 3e there are few choices. 
Given also the Na# of this vein. 0.143, it seems auite 
probable that the early evolution of this brine invi zed 
cooling and precipitation along the cotectic betweer. mi- 
rabilite and epsomite. If cooling continued these solid 
phases would react with the residual liquid at the peritec- 
tic, forming magnesium sulfate dodecahydrate: perhaps 
eventually the liquid would attain the eutectic involving 
precipitation of water ice. This evolution would result in 
just the sort of chemical trends observed, since Na and 
Ni would be depleted in the late-forming crystal assem- 
blage. It also indicates^that the brine was not in eqc ‘ib- 
rium with ice when the' highly sodic portion of the ’• ein 
was precipitated. This conclusion is consistent with the 
calculations of Zolensky et al. (1989), who concluded that 
Cl (Cl) chondrites were probably altered from CM (C2) 
or CV3 (C3) chondrite material at temperatures ~50 to 
150°C, while CM (C2) chondrites were altered near 0°C. 

5.3. Europa 

Europa has one of the youngest, smoothest, brightest, 
and geologically most active surfaces in the Solar System 
(Smith et al. 1979). Europa’s average albedo, 0.7 (Buratti 
and Veverka 1983), is exceeded in the Solar System only 
by the surfaces of Enceladus and Triton, both of which 
also have had geologically dynamic histories. Water ice 
reflectivity and absorption features in the visible and near 
infrared spectrum of Europa require water ice to be the 
dominant component of the visible surface (Clark 19 '0. 
Clark et al. 1986). Other, more subtle spectral features 
and albedo variations, however, indicate minor impurities 
including sulfur dioxide (Ockert et al. 1987). Europa ap- 
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parently suffered global volcanic resurfacing compara- 
tively recently. Indeed, Squyres et al. (1983) suggested 
that active and very widespread resurfacing may be oc- 
curring by frost deposition associated with water volca- 
n isr Europa’s mean density, 3.04 gem" 3 , is intermediate 
bet -en the densities of aqueously differentiated Cl 
chondrites (2.69 g cm" 3 ; model A. above) and Io (3.53 g 
cm' 3 ) and is just slightly greater than the density of C2 
chondrites. Squyres et al (1983) suggested that Europa 
contains a dense, mainly anhydrous rocky core and a 
water ice crust -100 km thick formed by dehydration of 
hydrous silicates. 

Model A (Fig. 9a), though strictly applicable for the 
aqueous differentiation of Cl chondrites, qualitatively 
sug . its a different structure and evolutionary history 
than modeled by Squyres et al . (1983). Europa originally 
may have resembled C2 chondrites, roughly approxi- 
mated here as a mixture of one part Cl and two parts 
C3 chondrites supplying 6 % water (Squyres et al. 1983). 
Thus. Europa contained one-third as much water, but also 
one-third as much salt as Cl chondrites. Differentiation 
may be approximated within the ternary system H : 0- 
MgSQ 4 -Na : S0 4 (Fig. 3a). The initial stages of igneous 
eve :ion should be identical to that described in model 
A, except that the volumes of each differentiated layer 
would be one-third that shown in model A. 

Europa shows clear indications of a prolonged and con- 
tinuing igneous evolution, which must require a more 
complex scenario than outlined in model A. Indeed, Eu- 
ropa is nearly as large as Earth’s Moon, where extensive 
igneous silicate differentiation occurred. It seems ines- 
capable that radiogenic heating plus tidal deformational 
hea*:ng would have elevated internal temperatures in Eu- 
rope to the point where gypsum and hydrous silicates 
would have dehydrated (Lewis 1971, Squyres et al. 1983). 
Thus, the remainder of Europa’s water associated with 
hydrated silicates should also be in its crust. After this 
later addition of water, Europa’s crust, — l 10 km thick, 
would constitute —8.8% of the satellite’s mass (6.0% H 2 0; 
2.0% MgS0 4 ; 0.69% Na : S0 4 ; minor sulfates of K, Mn, 
and Ni; 0.09% chlorides — probably mainly MgCl 2 ; 0.01% 
CaS0 4 ; and negligible Ca-Mg carbonates assuming 
slightly acidic to basic pH). Resumming the three major 
components to 100%, the bulk composition of Europa’s 
crust would be H 2 0 69%, MgS0 4 23%, and Na 2 S0 4 8%. 

Continued heating of Europa’ s crust by tidal dissipation 
would result in further differentiation of the crust. The 
road to complete differentiation would have been very 
complex, but once complete, Europa’s crust would be 
composed of a thin lower crust —15 km thick composed 
of anhydrous magnesium and sodium sulfates, and a thick 
(- l 5 km) ternary eutectic upper crust composed of —50% 
water ice, 44% MgS0 4 • 12H : 0, and 6% Na 2 S0 4 * 10H : O. 
Any further heat input would result in repetitive congruent 


melting of the eutectic crustal layer. While trace and minor 
components could fractionate during repeated melting, 
the major components would maintain the same eutectic 
composition in the liquid phase and solid residues, irre- 
spective of the extent of melting. The eutectic liquid is 
more dense than its coexiting solid assemblage, so erup- 
tive activity would tend to be inhibited. Instead, the liquid 
would pool between the overlying frozen eutectic saline 
crust and the anhydrous lower crust, forming magma 
chambers or a saline magma ocean. 

Complete freezing of lower crustal brine would be more 
difficult than in the salt-free case since solutes would 
depress the freezing points. While sulfates would mostly 
freeze out between 260 and 269 K, the aqueous eutectic 
chloride fraction of the brine would not completely freeze 
until about 221 K (for 1 atm: pressures of —1 kilobar 95 
km deep would depress the freezing point of the chloride 
solution to about 212 K). This could greatly reduce the 
rigidity of the upper crust, thus allowing prodigious rates 
of tidal dissipative heating and vigorous geologic activity 
to be maintained into the present epoch (Cassen et al. 
1980, Squyres et al. 1983). 

Crystallization of brines in the lower crust would result 
in volume expansion and possibly tensile fracturing of 
the upper crust. Lower crustal brines could be forcibly 
extruded under pressure during progressive freezing. 
Fractional crystallization seems likely, given the low vis- 
cosity of the brine and the high density differences be- 
tween the brine and all solid phases. Complex vertical 
mineralogical zonation of the crystallizing magma could 
occur. Remelting of these layers would produce liquids of 
many compositions, some less dense and others more 
dense than the overlying ternary eutectic crust. Extrusive 
activity would selectively involve the least dense, least 
saline brines, e.g.. aqueous sodium sulfate eutectic liquid 
formed by melting ice plus mirabilite. 

Brine eruptions, whether or not resembling the model 
outlined above, should have left physical evidence on 
the Europan surface. The low relief of Europa’s volcanic 
plains (compared, say, to Triton’s) is consistent with the 
low viscosities of brines (Kargel et al. 1991). It is some- 
what surprising that spectral evidence of salts has not 
been reported, as it seems difficult to avoid melting any 
part of Europa. at any stage in its evolution, without 
involving substantial quantities of salts. However, water 
evaporated from erupted brines would tend to recondense 
widely on the cold surface of Europa (Squyres et al._ 1983) 
and might tend to hide the spectral signature of the ex- 
pected magnesium and/or sodium sulfates and other salts. 
The ultraporous microstructure of Europa’s surface as 
inferred from photometry (Domingue et al. 1991) is consis- 
tent with such frost deposition. Furthermore, cryovol- 
canic salts would be hydrated and mixed with water ice, 
so the strong spectral signature of intimately mixed ice 
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and chemically bound water may mask the spectral signa- 
ture of salts. 

5.3. Ganymede 

Geology and composition . Ganymede’s average den- 
sity, 1.93 g cm' 3 , is much lower than those of lo (3.55) 
and Europa (3.01) (Schubert et ai 1986). Ganymede's 
density falls between the densities of compressed (non- 
porous) Cl chondrites (2. 6-2. 8), and compressed water 
ice (ranging from 0.92 to 1.4 for various relevant poly- 
morphs), indicating that Ganymede, like the other icy 
satellites, contains substantial fractions of both rock and 
ice (Lupo and Lewis 1979, Schubert et al. 1986, Mueller 
and McKinnon 1988). This simple but important fact 
supports well-based theoretical predictions that ice 
should be an abundant low-temperature condensate in 
the outer Solar System and specifically in Ganymede 
(Lewis 1971, 1972). Telescopic observations confirm that 
ice is the major constituent of Ganymede’s surface, 
although there is also an important rocky component 
(e.g., Clark et al . 1986). 

The estimation of the ice and rock fractions in icy 
satellites has been an important objective of many 
workers. Schubert et al. (1986) and Mueller and McKin- 
non (1988) clearly showed that without tighter con- 
straints on the composition of the rock and on Gan- 
ymede’s state of differentiation, one cannot derive a 
unique rock fraction. This important parameter may be 
as low as 0.49 to as high as 0.65. However, the geology 
of Ganymede’s surface requires at least some degree 
of internal differentiation. Extensive if not complete 
differentiation seems probable. Modeling the case of 
extensive differentiation, Mueller and McKinnon (1988) 
calculated rock fractions ranging from 0.54 to 0.59 
depending on the type of rock composing Ganymede’s 
core. However, the authors pointed out that regardless 
of the core’s present composition it very likely evolved 
from an initial rock component resembling Cl chon- 
drites. If the rock fraction is expressed in terms of the 
initial fraction of Cl chondrite rock, a differentiated 
Ganymede would have contained 64 to 65% of its mass 
as rock, the remainder being ice. The latter set of rock 
fractions, 0.64 to 0.65, corresponds to R t in the models 
shown in Fig. 9. The quoted range of 0.54 to 0.59 would 
roughly correspond to R f , the final rock fraction after 
soluble salt extraction; however, this correspondence is 
not exact, since the models of Mueller and McKinnon 
(1988) leave the anhydrous salts in the rock while 
extracting only their water content. 

Pre-Voyager theoretical models showed the likelihood 
that Ganymede would be differentiated (Consolmagno 
and Lewis 1978). This prediction is confirmed even in 
low-resolution full-disc images showing bright and dark 
terrains on Ganymede’s surface. Higher resolution im- 


ages reveal terrains with considerably varied crater 
densities and intensely fractured and volcanicalh 
flooded terrains, indicating a complex history of differen- 
tiation (Smith et al. 1979, Head et al. 1981, Wor now 
et al. 1982, Squyres 1982, Parmentier et al. 1982, Casac 
chia and Strom 1984, Croft 1985 and 1986, McKinnon 
and Parmentier 1986, Schenk and McKinnon 1987, Croft 
and Goudreau 1987, Golombek and Banerdt 1986, Alli- 
son and Clifford 1987, Murchie and Head 1988, Mueller 
and McKinnon 1988, Tittemore 1990). Bright terrain 
typically intrudes older dark terrain (Parmentier et al. 
1982, Squyres 1982), and usually the light terra: a is 
heavily deformed by one or more sets of extenv.mal 
and probably also transtensional tectonic features (Mur- 
chie and Head 1988). 

Volcanism and extensional tectonism are closely inter- 
related on Ganymede as they are on Earth. Fracture 
zones are frequently the sources, and tectonic depres- 
sions are commonly the containment basins for erupted 
plains-forming liquids (Head et al. 1981, Parmentier ei 
ai 1982, Allison and Clifford 1987). In many -eas 
tectonic grooving was followed by cryovolcanic floe ring, 
in turn sometimes followed by reactivation of grooving 
and reflooding. Thus, near the end of a period of 
heavy impacting, Ganymede seems to have suffered late 
episodes of intense tectonic and cryovolcanic surface 
modification (Woronow et al. 1982). 

An absence of obvious flow margins in the best 
Voyager images of Ganymede's smooth plains indicates 
that flow thicknesses must be <100 m. Hence, the 
erupted liquids possessed low viscosities and rield 
strengths relative to much thicker flows on the satellites 
of Uranus and Neptune. After considering the 1/7$ 
surface gravity, the 'erupted substances may have been 
Theologically comparable to any substance as “thin" as 
or “thinner" than andesite (Kargel 1990). Of course, 
silicate lavas are disallowed in Ganymede’s situation, 
so Theologically compatible and compositionally likely 
substances would include water, salt brines, or ammo- 
nia-water liquids, the latter only as long as there are 
low contents of additional components and suspended 
crystals (Kargel et al. 1991). Since substantial quantities 
of ammonia are unlikely to have condensed in the 
relatively warm Jovian nebula, water or salt brines are 
more likely than ammonia-water liquid. 

The freezing-point depression caused by salts also 
may assist in melting in a satellite that may convectively 
regulate its internal temperatures below the me : ting 
point of pure ice. Magnesium sulfate alone would only 
lower the melting point by 4 K, but certain other salts, 
particularly the chlorides, are more effective freezing 
point depressants. Once volcanism initiates the fraction 
of rock in the solid residual core increases, causing the 
viscosity of the core to increase, resulting in higher 
convective temperatures or even a complete failure of 
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the convective mechanism, thereby promoting further 
mating. 

internal structure of a differentiated salty Gany- 
mede . Figures 9C and 9D bracket Ganymede’s compo- 
sition. Of course, these models do not consider the 
important effects of elevated pressures in Ganymede’s 
interior, nor of qualitatively different scenarios involving 
possible fractional crystallization of magma oceans. De- 
spite this outrageous neglect, it sems probable that 
Ganymede would have some type of compositionally 
layered structure; based on these models the structure 
w ild include a MgS0 4 -rich mantle —500 to 600 km 
thick, containing -15-20% by mass MgS0 4 and Na 2 S0 4 , 
plus minor quantities of other salts. This layer would 
be overlain by less salt-rich icy upper mantle and crustal 
layers an additional several hundred kilometers thick. 

Pressures near the core/mantle boundary would be 
-12-13 kbars, so the mantle and crust also would be 
layered with respect to ice polymorphs, including phases 
up to ice VI near the base of the mantle. The effects 
o*. pressures on sulfate phase stability are uncertain, 
although sparse and conflicting experimental data sug- 
gest a phase change in epsomite near 5 kbars and 
perhaps another at pressures over 10 kbars (Bridgman 
1948a, b, Livshits et aL 1963). Mantle and crustal compo- 
sitions reflect phase equilibria during partial melting 
which would have occurred under pressures ranging as 
high as 50 kbars near the center of Ganymede. Such 
high pressures would certainly impact subsolidus and 
solid-liquid phase equilibria in ways that are difficult to 
predict. Therefore, Fig. 9 is not at all a reliable guide 
to Ganymede’s interior structure, but is a valid indica- 
tion of the probable importance of salts in Ganymede’s 
and other icy satellite’s histories and interior structures. 

Taking Figs. 9c and 9d literally, Ganymede’s surface 
would be composed of volcanic flows of water or sodium 
sulfate-water solutions. However, in reality it would be 
difficult to extrude these liquids through the magnesium 
sulfate-rich mantle without substantial contamination by 
magnesium sulfate. Also* if Ganymede’s Na# < 0.15 
then sodium sulfate would be exhausted before magne- 
sium sulfate, causing a different petrologic evolution of 
late-stage liquids. Coupled with the unknown pressure- 
related effects on phase equilibria, it is not yet possible 
to predict with any certainty what Ganymede’s surface 
composition should be. Even so, it seems a priori 
highly probable that salts, particularly magnesium and/or 
sodium sulfates, would be significant components. 

6. SUMMARY 

Brine volcanism may explain resurfacing observed on 
some icy satellites, particularly Europa and Ganymede, 
and is expected to have been important in the evolution 


of chondritic asteroids. Carbonaceous chondrite mineral- 
ogies, cosmic elemental abundances, and salt solubilities 
suggest that MgS0 4 and Na : S0 4 are the most important 
solutes. Internal chemical structural models of icy satel- 
lites and asteroids suggest a wide range of possible evolu- 
tionary scenarios, some involving the formation of thick 
salt-rich crusts and mantles. 

Salt-water eutectic brines have densities lower than 
typical ice-rock mixtures, and therefore are buoyant and 
can be extruded onto the surfaces of icy satellites dur- 
ing early stages of differentiation. The ternary MgS0 4 - 
Na 2 S0 4 -H : 0 eutectic liquid has a density of -1.19 g 
cm* 3 , compared to 1.13 g cm* 3 for the equivalent frozen 
eutectic mixture; therefore, a solid eutectic crust would 
float on its own liquid. MgS0 4 -rich aqueous volcanism 
would become increasingly difficult as a eutectic-compo- 
sition crust became thicker. Early eutectic volcanism 
may generally yield to later plutonism or even formation 
of stable brine magma oceans. Late-stage water or 
Na 2 S0 4 -H 2 0 eutectic volcanism would be favored, since 
these liquids are less dense than the MgS0 4 -rich eutectic 
crust. Such late-stage water-rich volcanism could occur 
if ice and mirabilite are left in excess amounts in the core 
after exhaustion of magnesium sulfate. 

Hypersaline brines have much lower viscosities than 
familiar silicate lavas, although they are more viscous than 
pure water. Aqueous sulfate flows may produce extensive 
flows on the order of a meter thick, constructing vast 
volcanic plains of low relief such as observed on Gan- 
ymede and Europa. Chloride brines are generally more 
viscous than sulfate brines. Under low-g conditions on 
icy satellites and asteroids chloride brines would probably 
behave fluid dynamically much as basalt behaves on 
Earth, producing characteristic flow thicknesses of —2-20 
m. Chemical as well as geomorphic evidence should exist 
on the surfaces of objects affected by brine volcanism, and 
might be discemable in observations planned for Galileo. 
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